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ABSTRACT 


A FUNDAMENTAL STUDY OF NUCLEATE POOL 
BOILING UNDER MICROGRAVITY 


by 

Jamie Ervin 
Herman Merte, Jr. 


An experimental study of incipient boiling in short-term microgravity and with a/g 
= ± 1 for pool boiling was performed. Calibrated thin gold films sputtered on a smoothly 
polished quartz surface were used simultaneously for thermal-resistance measurements and 
heating of the boiling surface. The gold films were used for both transient and quasi- 
steady heating surface temperature measurements. Two test vessels were constructed for 
precise measurement and control of fluid temperature and pressure: a laboratory pool 
boiling vessel for the a/g = ± 1 experiments and a pool boiling vessel designed for the 1 3 1 
m free-fall in the NASA Lewis Research Center Microgravity Research Facility for the 
microgravity tests. Measurements included the heater surface temperature, the pressure 
near the heating surface, the bulk liquid temperatures. High speed photography (up to 
1,000 frames per second) was used in the experiments. 

With high quality microgravity and the measured initial temperature of the quiescent 
test fluid, R113, the temperature distribution in the liquid at the moment of boiling 
inception resulting from an imposed step in heat flux is known with a certainty not possible 
previously. The types of boiling propagation across the large flat heating surface, some 
observed here for the first time, are categorized; the conditions necessary for their 
occurrence are described. Explosive boiling propagation with a striking pattern of small 
scale protuberances over the entire vapor mass periphery not observed previously at low 
heat flux levels (on the order of 5 W/cm 2 ) is described. 

For the heater surface with a/g = -1, a step in the heater surface temperature of 
short duration was imposed. The resulting liquid temperature distribution at the moment of 
boiling inception was different from that obtained with a step in heat flux. 

This report constitutes a major part of a dissertation submitted by the first author in 
partial fulfillment of the requirements for the degree of Doctor of Philosophy in the 
University of Michigan, 1991. 
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CHAPTER 1 


INTRODUCTION 

The design of power generation systems and the extra-terrestrial storage of 
cryogenic fuels in microgravity for the U.S. space program require a knowledge of 
boiling and nucleation of the vapor phase. Boiling in microgravity is fundamentally 
different from boiling in earth gravity: the buoyancy force which induces liquid and 
vapor motion in boiling with earth gravity is effectively eliminated in microgravity. 
The study of boiling in microgravity can give some basic insight into the boiling 
process in the earth and space environment. 

Because of the sometimes imprecise use of terms in connection with boiling 
research, it is desirable to clearly define a number of terms to be used here. The 
following discussion is for a pure substance, one that has a homogeneous and 
invariable chemical composition throughout. 

Nucleation, pertinent to any discussion of boiling, is the formation of an 
interface between an initial parent phase and a second phase. Nucleation is classified 
as either homogeneous or heterogeneous. Homogeneous nucleation refers to 
nucleation taking place within a pure liquid bulk in the absence of any foreign 
particles or pre-existing phases at a uniform temperature as a result of local random 
density fluctuations. These density fluctuations can result in agglomerations of a 
sufficient number of molecules to form a nucleus of critical size such that subsequent 
growth can occur. (The thermodynamic requirements for this critical size will be 
given in Chapter Two.) Heterogeneous nucleation occurs at the interface between the 
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fluid and other materials, usually solid. The conditions necessary for heterogeneous 
nucleation to occur are not well understood, and may be quite different for disparate 
interfaces, such as smooth rigid surfaces (glasses for example), liquid-liquid interfaces, 
and non-smooth solid surfaces (crystalline metallic surfaces for example). 

The terms incipient boiling, initiation of boiling, and onset of boiling all refer 
to the appearance of the first observable vapor bubble. Incipient boiling can result 
from either homogeneous or heterogeneous nucleation. The first observable bubble 
may be the consequence of growth from either a pre-existing nucleus or a newly 
formed nucleus. 

Nucleation and boiling of liquids generally takes place by heating the liquid at 
a solid boundary, so that the liquid does not have a uniform superheat as under the 
circumstances of a pressure decrease. The heating of a liquid to a critical superheated 
state from a solid boundary until the onset of boiling occurs may be performed in a 
transient or a quasi-steady manner. Transient heating is inherently unsteady, while 
quasi-steady heating to the onset of boiling results from the application of a series of 
incremental steps in heat flux to a system where a steady temperature distribution is 
obtained following each of these steps until nucleation occurs. For example, a heated 
flat surface facing upward in a liquid pool in standard gravity results in natural 
convection, which permits a series of steady-state increments in heat flux before 
boiling first occurs. A quasi-steady heating process is not possible in the true absence 
of gravity, and may or may not be possible in microgravity, depending on the imposed 
heat flux level. 

This is a study of incipient boiling in which body forces of a/g = ± 1 and 
microgravity are used as a means of varying the temperature distribution in the liquid 
at the time of nucleation with all other parameters held constant. The liquid 
temperature distribution in microgravity, since buoyancy induced liquid motion is 
drastically reduced, is well defined relative to the liquid temperature distribution in a/g 
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= ± 1. The conditions in microgravity existing at the moment of boiling inception will 
be compared with those at standard earth gravity. This study is limited to short term 
microgravity conditions, because of the practical limitations of achieving long term 
microgravity at the present time. 

Incipient pool boiling was investigated in an attempt to assess the influence of 
liquid velocity under circumstances where buoyancy is anticipated to play a role. 
Short duration microgravity was available only for the pool boiling experiments. An 
imposed constant heat flux, constant surface temperature, and quasi-steady heating 
constitute the three different boundary disturbances used in the pool boiling 
experiments. 

The principal objectives of this work are: 

(1) To determine if known temperature distributions at the onset of boiling 
can be used to characterize the incipient boiling process. 

(2) To examine how the onset of boiling is influenced by the method of 
heating using an imposed heat flux, an imposed surface temperature, 
and a quasi-steadv heating procedure. 

(3) To explore the features of observed phase change dynamics, and to 
characterize the parameters which produce particular types of vapor 


motion. 
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CHAPTER 2 
LITERATURE SURVEY 

Homogeneous nucleation in a pure liquid phase is unlikely in normal 
engineering applications because heterogeneous nucleation usually takes place before 
the large superheats necessary for homogeneous nucleation are attained. A metallic 
heating surface, for example, has numerous pits which not only serve as reservoirs for 
trapped gases, but also expose more grain boundary area (grain boundaries are regions 
of higher energy) to the liquid than would otherwise be available if there were no pits 
present. 

For nucleation of a heated liquid to be homogeneous, the probability for the 
agglomeration of sufficiently energetic molecules to form a nucleus of critical size in 
the interior of the liquid must be greater than the probability for agglomeration of 
energetic molecules at its surface. A necessary condition for this to occur is that the 
heated liquid not be in contact with any vapor phase. According to Blander and Katz 
(1975), another condition which is sufficient but not necessary for homogeneous 
nucleation is that the liquid wet any other phase with which it is in contact. Jarvis et al 
(1975) have shown that if a liquid has a zero contact angle with any body that it 
contacts, then a bubble in the interior of the liquid is more stable than a similar size 
bubble at the interface. Hence, homogeneous nucleation is more likely to take place 
than heterogeneous nucleation. These conditions are met when the liquid is 
completely enclosed by a smooth surface. Examples are glass or another liquid which 
is non-volatile, as discussed in Blander and Katz (1975). 
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Homogeneous nucleation theory and its extension to heterogeneous nucleation 
is basic to any understanding of incipient boiling. A brief review of homogeneous 
nucleation theory will be presented first, followed by a discussion of heterogeneous 
nucleation. 

2.1 Homogeneous Nucleation 

Homogeneous nucleation may be considered from two perspectives. The first 
is the classical theory, attributed to Gibbs (1928), Volmer (1945), and Becker and 
Doring (1935), and refined by Frenkel (1955), and Zeldovich (1942). The second is 
the statistical mechanical point of view, associated with Reiss (1952a) (1952b), Mayer 
and Harrison (1938), Lothe and Pound (1966), and Abraham (1974), and others. The 
motivation behind the development ' f the latter perspective is that the assumption of 
the applicability of macroscopic thermodynamic properties to a small vapor bubble 
(which may be on the order of 25 molecules) was not deemed appropriate. However, 
it is difficult to obtain numerical results with the statistical approach and to relate these 
to experimentally measurable variables. As a way to circumvent these problems, some 
workers have combined concepts from both macroscopic and microscopic 
thermodynamics, bringing statistical thermodynamics into the classical expressions in 
the form of a correction factor, as in Lothe and Pound (1962). However, disagreement 
exists over the validity of this type of correction factor, as discussed in Reiss (1952), 
Lothe and Pound (1966), Reiss (1970), Reiss and Katz (1967), and Lothe and Pound 
(1968). The classical theory has been used successfully in predicting homogeneous 
nucleation temperatures for many fluids, as in Skripov (1974), Springer (1978), and 
Kotake and Glass (1981). Only the classical theory will be presented in the present 
work, as it provides insight into the basic physics of homogeneous nucleation. 
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Homogeneous nucleation of the vapor phase in the absence of solid surfaces, 
impurities, or ionizing radiation occurs within the bulk of the liquid as a result of 
molecular fluctuations in an isothermal system. These fluctuations give rise to a 
population of vapor embryos or agglomerations of molecules which may be described 
by a Boltzmann-like distribution of sizes, as discussed in Frenkel (1955): 


N(r) = A exp 


-AG(r) 

kT J 


( 2 . 1 ) 


N is the number of spherical embryos of radius r. AG is the change in the Gibbs free 
energy of an isolated isothermal system associated with the formation of a spherical 
vapor nucleus in the metastable superheated liquid. For the assumption of a spherical 
vapor nucleus in a liquid of uniform temperature, the Laplace equation provides the 
equilibrium pressure difference across a curved interface arising from the effect of 
surface tension: 


P 



With equation 2.2, an expression for AG can be derived: 


AG=r* 2 47ia 




( 2 . 2 ) 


(2.3) 


The first term on the RHS of equation 2.3 represents the surface energy contribution 
related to the formation of a spherical vapor liquid interface. The second term on the 
RHS of equation 2.3 is associated with the formation of a volume of the vapor phase, 
r* is the critical size of the embryo, above which spontaneous growth of the 



7 


agglomeration occurs. Below r*, however, the embryo collapses. The surface tension 
is assumed to be a function of T only, which may be valid for agglomerations of 
sufficient size, as explained in Springer (1978). For phase transitions, the Gibbs free 
energy of the system is a minimum for stable equilibrium, while it is a maximum for 
metastable equilibrium since the Gibbs free energy of the system decreases for either 
an increase or a decrease in the size of an embryo of critical size, as discussed in 
Kottowski (1973). The second derivative of equation 2.3 evaluated at r* is negative, 
demonstrating that a superheated liquid is in a metastable state. An embryo of critical 
size corresponds to a maximum AG, or AG* below. Hence, with r = r* in equation 
2.3: 


AG* = 4/3 it or* 2 


(2.4) 


with equation 2.2: 


AG* = 


167tg 3 

3(P V - P^ 2 


(2.5) 


The nucleation rate, or the rate of formation of embryos of critical size is given 
by Fisher (1948) as: 


J = 


nkT 


-AG*]' 

h 

exp 

L kT JJ 


( 2 . 6 ) 


This expression for the nucleation rate neglects the term for the free energy of 
activation of the motion of an individual molecule of liquid past its neighbors into or 
away from the liquid-vapor interface. AG* may be considered as a thermodynamic 
energy barrier for the liquid-to-vapor transition. In classical homogeneous nucleation 
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theory, vapor embryos are assumed to have the macroscopic properties of the liquid 
bulk phase. Hence, the surface tension in the above equations is the macroscopic flat- 
film surface tension. This assumption is often referred to as the capillary 
approximation. Different versions analogous to equation 2.6 exist but provide similar 
calculated results, as discussed by Skripov (1974). For simplicity, equation 2.6 will be 
used in the present work to represent J. The pre-exponential factor, ^ , is a large 

number, on the order of the number of molecules in the volume of liquid considered. 
T in equation 2.6 may be defined as the homogeneous nucleation temperature, denoted 
b y T hom> for a critical range of J, as discussed in Skripov (1974) and Nghiem (1980). 
Thom actually covers a very narrow range of T and vanes little even for a large change 

in the critical value of J. The maximum practical attainable superheat of a liquid is 
limited by T hom< as discussed in Blander and Katz (1975) and Sinha (1980). It is 

noted from equations 2.5 and 2.6 that the exponential term of equation 2.6 is 
temperature dependent. As the temperature of a pure liquid is raised above its boiling 
point, the surface tension decreases and the vapor pressure increases. Each of these 
two factors leads to a decrease in the activation energy barrier, resulting in an increase 
in the likelihood of vapor nucleation. 

Experimental studies of homogeneous nucleation abound and the following are 
selected as a representative sampling: 

Simpson and Walls (1965-1966) used a pulsed wire technique to investigate 
nucleation. Their homogeneous nucleation temperature measurements were within + 2 
% of the theoretical homogeneous nucleation temperature of various organic fluids. 
The experimental apparatus consisted of a thin platinum wire immersed in different 
test fluids, including benzene, acetone, n-hexane, ethyl ether, and demineralized water 
at atmospheric pressure. The wire was subjected to an intense current pulse of short 
duration, so that nucleation would occur in the superheated thermal boundary layer of 
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the liquid near the wire surface. The transient heat flux from the wire was interpreted 
to be convection-free conduction to the liquid. 

Skripov (1974) studied homogeneous nucleation in saturated hydrocarbons and 
fluorohydrocarbons by superheating floating droplets of saturated hydrocarbons and 
fluorohydrocarbons in sulfuric acid in a temperature gradient on the order of 0. 1-0.5 
°C/mm. Since the hydrocarbons were immiscible in the sulfuric acid, the nucleation 
taking place approached homogeneous nucleation. The homogeneous nucleation 
temperatures calculated by classical homogeneous nucleation theory agreed with 
measurements (differences were typically less than 1-2 °C). 

Skripov (1974) also studied homogeneous nucleation in an expanding liquid 
achieved by means of a rapid pressure drop in a bubble chamber. He observed mean 
homogeneous nucleation temperatures which were within a few degrees Celcius of 
those predicted by classical homogeneous nucleation theory for n-pentane, n-hexane, 
diethyl ether, benzene, and hexafluorobenzene for various saturation pressures. 

Derewnicki (1985) performed boiling tests using pulsed platinum wires (0.025 
mm diameter by 10 mm length) immersed in degassed, distilled water. He claimed to 
find, by photographic means, homogeneous nucleation in the thermal boundary layer 
near the wire surface which became highly superheated as a result of the fast 
temperature rise. 

Sinha (1980) measured the homogeneous nucleation temperature for Helium I 
using a transient bismuth magnetoresistive thermometer. Homogeneous nucleation 
temperatures for liquid nitrogen for reduced pressures, P/P c , ranging from 0.03 to 0.91 

were obtained by pulse heating a platinum wire (7 cm length by 0.102 mm diameter). 
The experimental temperature measurements were in good agreement with the 
classical nucleation theory. A corresponding state analysis for the case of 
homogeneous nucleation in liquid noble gases, both classical and quantum, was 
developed. 
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2.2 Heterogeneous Nucleation and Incipient Boiling 

Heterogeneous nucleation, which is more probable in most engineering 
applications of boiling than homogeneous nucleation, is influenced by the geometrical, 
chemical, and physical properties of the boiling surface. Heterogeneous nucleation in 
boiling is extremely difficult to formulate mathematically as it is different for any 
given liquid-surface combination and depends on fluid contact angles, surface cold 
working (for metallic surfaces), grain boundary area (for metallic surfaces), and 
surface chemadsorption. In addition, minute cavities which might not be entirely 
wetted by a liquid may trap air or other gas molecules as potential centers for 
nucleation. 

Several models for heterogeneous nucleation which utilize the macroscopic 
point of view based on the concept of surface tension have been proposed: Fisher 
(1948), Hirth and Pound (1963), Sigsbee (1969), Boucher (1969), Cole (1974), 
Zimmels (1976), Blander and Katz (1975), and Kotake and Glass (1981). The 
discussion in Blander and Katz (1975) will now be summarized, presenting the basic 
concepts of the macroscopic heterogeneous nucleation. 

A bubble nucleating on a smooth liquid-solid interface is shown in Figure 2.1. 
The minimum work to form a bubble of critical radius in chemical and mechanical 
equilibrium is given by: 


AG* = a. A. + 

lv lv 


' 


r 

a -a . 
sv si 

A - 

sv 

P -P. 

V 1 

v > 




(2.7) 


If the bubble is assumed to be a spherical section, the bubble volume and interfacial 
areas can be expressed in terms of the radius of curvature, R, and contact angle, <>. 
From the geometry of Figure 2. 1 : 
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V = ^(R-Z) (2R+Z) = 2-3x|H-\|/ 


A = 27ir(R-Z) = 2rtR ( 1-y) 


2 2 
A = 7TT = TtR 

sv 


1-V 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


By a force balance at the bubble edge, 


o s , = o sv +aj v cos(7t-(})) = a sv +W Iv 


( 2 . 12 ) 


Hence, 


a ,-a 

si SV 


¥ : 


(2.13) 


'lv 


Substituting equations 2.8 through 2.13 and the Laplace equation into equation 2.7 
results in: 


4 2 

AG* = TJra lv R F 


(2.14) 


where: 
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2 - 3\|/+y~ 


F = ■ 


(2.15) 


Equation 2.14 reduces to equation 2.4 with F = 1. An expression for the nucleation 
rate of the form similar to the combination of equations 2.6 and 2.4 can be obtained by 
substituting the AG* of equation 2.14 for AG* in equation 2.6. However, Blander and 
Katz (1975) suggest that J is proportional to n^ rather than to n because molecules at 
the solid-gas interface are more likely to become potential bubbles than molecules 
within the volume of a liquid droplet. Thus: 


j n^kT 

J= ~ ~ exp 


f 4 2 

| V a lv RF 


kT J 


(2.16) 


The ability of the above macroscopic model of Blander and Katz (1975) to 
predict the conditions necessary for heterogeneous incipient boiling and the resulting 
phase change dynamics is limited. For example, the microgeometry necessary for the 
determination of the F in equation 2.16 is unknown for most surfaces. In addition, the 
effects of gases adsorbed on the boiling surface are not accounted for in such a model. 

Having examined the fundamental mechanisms of homogeneous and 
heterogeneous nucleation, the conditions influencing sustained heterogeneous incipient 
boiling in which bubble nuclei are formed or pre-exist during heating of the liquid will 
now be addressed. Models developed to predict incipient boiling have demonstrated 
limited degrees of success when compared with experimental measurements as a result 
of their inability to satisfactorily account for the following elements influencing the 


process: 
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(1) Surface microgeometry 

(2) Surface energy 

(3) Liquid temperature distribution 

(4) Externally forced liquid motion 

(5) Buoyancy induced liquid motion 

The heater surface microgeometry, which is assumed to consist of microscopic 
crevices and pits, traps vapor and gases which become bubble nuclei. The heater 
surface energy affects the wettability of a surface by a liquid, and as a result, 
influences the ability of the surface microgeometry to trap vapor and gases. The liquid 
temperature distribution at the moment of incipient boiling depends on the presence of 
fluid motion, externally forced or resulting from density variations in a gravity field. 
These five elements will now be considered in some detail. 

2.2.1 Surface Microgeometry 

It is known, as indicated by the studies below, that nucleation and initial 
boiling are strongly affected by the surface microgeometry. Bankoff (1958), Griffith 
and Wallis (1960), Lorenz (1974), Winterton (1977), and Webb (1981), for example, 
purport that the role of microgeometry in incipient boiling is to trap gases and vapor, 
which then act as nuclei for bubbles. Howell and Siegel (1966) report that 
photographs of water boiling from artificial sites drilled in stainless steel surfaces 
show vapor remaining in the site after a bubble departs. The high speed photographic 
studies of surface boiling of pentane by Clark et al (1959) demonstrated that active 
nucleation sites were located at pits and scratches on the surface. Griffith and Wallis 
(1960) showed that the cavity mouth radius determines the superheat necessary for 
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conditions of isothermal incipient boiling, and that the shape of the cavity determines 
whether boiling will persist, once having begun. 

Coeling and Merte (1968) presented experimental data for natural convection 
heat transfer and for boiling inception of LH 2 and LN 2 . The experimental variables 

included heater surface material, surface roughness, and orientation. Stainless steel 
and copper surfaces were tested in both polished and lapped conditions. A surface 
coated with TFE teflon was also tested. Heater surface orientation, which included 
horizontal up, vertical, and horizontal down with respect to the direction of earth 
gravity, were test parameters. A stainless steel surface lapped with a 600 grit 
compound was found to be a more effective boiling surface for LH 2 than the same 

surface whether polished or lapped with a coarser 280 grit compound. In terms of 
vapor trapping in cavities, this experimental result indicated the existence of an 
optimum surface microgeometry for producing nucleate boiling with a minimum 
heater surface superheat. More recent works concerning the role of microgeometry in 
incipient boiling and nucleate boiling will be discussed below. 

Webb (1981) presented a survey of the evolution of special surface geometries 
that promote high performance nucleate boiling covering the period from 1931 to 
1980. The surface geometries described were designed with the intent of trapping 
vapor efficiently. These geometries included re-entrant cavities, porous coatings, and 
re-entrant grooves. 

Forest (1982) presented a general free energy analysis of the equilibrium and 
stability of a gas-vapor nucleus inside of a surface cavity for a uniform temperature in 
the liquid. He demonstrated that the spontaneous growth of a pre-existing nucleus 
occurs when the radius of curvature reaches a minimum as the gaseous volume of the 
nucleus increases. This confirmed the earlier work of Griffith and Wallis (1960). 

Tong et al (1990) described the influence of the dynamic contact angle on the 
liquid-gas trapping as related to boiling inception with highly wetting liquids. 
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Changes in the dynamic contact angle were shown to determine the trapped volume of 
gas in a heater surface crevice. 

The major problem encountered in studying the role of surface microgeometry 
on incipient boiling for a "natural" boiling surface is how to satisfactorily quantify, in 
a collective manner, the micro-valleys, pits, and scratches present on such a surface. 
Surface roughness measurements are gross measurements for the length scale of the 
cavities being considered, which are on the order of 10‘ 7 m. Surface roughness 
measurements do not necessarily detect vapor trapping re-entrant cavities or 
microscopic pits, while scanning the boiling surface at discrete locations may not 
reveal some potential boiling sites. 

Chowdhury and Winterton (1985) studied the effect of surface roughness and 
contact angle on nucleate boiling using a quenching technique with aluminum or 
copper cylinders and saturated water or methanol. They concluded that the number of 
active sites, rather than the surface roughness as measured by a profilometer, 
influenced nucleate boiling. Consistent techniques for roughening different surfaces 
of the same materials produced remarkably similar boiling curves for surface-liquid 
combinations with the same measured contact angle. 

2.2.2 Surface Energy 

The importance of the surface energy, measured in terms of the contact angle, 
to the onset of boiling has long been known. However, surface energy is not well 
understood. It is not known explicitly how surface energy affects the formation of 
nuclei, nor how surface energy affects the stability of nuclei in potential boiling 
cavities. The contact angle between a solid and a liquid characterizes the wettability 
of a certain solid surface by a specific liquid. Parameters such as the volume of 
trapped air or gas in a surface cavity, the critical bubble radius and the resulting 
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incipient superheat are influenced by surface wettability. The description of this 
wetting phenomenon is difficult, as solid surfaces are surmised to have micro crevices 
and valleys of complex forms. In addition, chemical reactions may take place at 
liquid-solid interfaces and develop impurities which concentrate at the liquid-solid 
interface. The following studies examined the relationship between the contact angle 
and boiling incipience. 

Cornwell (1982) maintained that boiling incipience was caused by the 
instability of pre-existing vapor-liquid interfaces in heater surface cavities. He 
proposed that the micro-roughness of the cavity walls results in a contact angle 
hysteresis between the advancing and receding contact angles which cause the 
meniscus within the cavity to invert from being convex to concave on the vapor side. 
This meniscus inversion provides an explanation for the trapping of vapor and gases. 
Surface heterogeneity in terms of spatially varying amounts of contaminants, including 
organic films and oxide layers, was purported to play a pan in the contact hysteresis. 

Winterton and Blake (1983) concluded that the liquid- vapor interface within a 
boiling surface cavity can adjust to external changes in liquid pressure on a time scale 
of less than a few seconds. They purported that boiling incipience on a heated surface 
induced by surface vibration can be explained in terms of the hysteresis between the 
advancing and receding contact angles. 

In a theoretical investigation, Tong et al (1990) concluded that for highly- 
wetting fluids such as R113, the nucleus for boiling incipience lies below the cavity 
mouth radius. They also concluded that for such fluids, variations in the contact angle 
during bubble embryo formation and growth resulting from changes in the liquid- 
vapor interface velocity explain the wide variation in measured incipient superheat 
values. A method for calculating the wall temperature at boiling incipience was given, 
but depended on unknown values of the static equilibrium contact angle, the dynamic 
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contact angle, and the cavity geometry. As these values could not be measured, 
assumed values were used. 

Because the heating surface microgeometry and surface energy are usually 
unknown, some studies have lumped these two effects in an energy reduction factor as 
a modification to classical expressions for homogeneous nucleation in order to predict 
boiling inception. 

Nghiem et al (1981) attempted to determine the necessary conditions for 
transient incipient boiling with a step in power input using two test fluids, saturated 
water and R1 13 at atmospheric pressure, and using three test surfaces: a platinum wire 
(0.15 mm diameter x 24 mm length), a copper wire (0.2 mm diameter x 24 mm 
length), and a thin gold film (20 mm x 2 mm x 400 angstroms thick). It was observed 
that as the steady input heat flux was increased, the delay time for the onset of boiling 
decreased. This decrease in delay time was accompanied by higher wall superheats at 
the moment of boiling inception, up to 90 K at 0.9 sec for the R113. To incorporate 
the influences of convection, fluid properties, surface energy, and microgeometry of 
the heated surface on the delay time, an empirical factor was used in combination with 
an expression resulting from homogeneous nucleation theory. This factor was defined 
as: 


F = 



I AG* J 


(2.17) 


where AG*, the free energy required to form a vapor nucleus for the process of 
homogeneous nucleation, was defined in equation 2.4. AGj, et was defined as the free 

energy required to form a vapor nucleus under conditions corresponding to 
heterogeneous nucleation. The expression of Fisher (1948) for the homogeneous 
nucleation rate of bubbles of critical size, given as equation 2.6 above, was used to 
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relate the nucleadon rate to the observed nucleadon delay time. Nucleadon is 
considered to occur when a bubble of critical size forms and grows per unit volume in 
a reasonable time, t seconds, after the onset of heating. This t was termed the 
nucleation delay time, and J in equation 2.6 was taken as the inverse of this delay time, 
as a first approximation. Since the heterogeneous nucleation appeared to take place in 
the fluid immediately adjacent to the solid surface, the fluid temperature at the instant 
of nucleation was taken to be the same as the wall temperature. For a given fluid- 
surface combination, F was a function of the fluid properties and the measured heater 
surface temperature at the moment of incipient boiling. This particular treatment of 
incipient boiling will be considered in more detail in a later section. 

A free energy reduction factor based on homogeneous nucleation theory 
paralleling that of Nghiem et al (1981) was also presented in Gerum et al (1979), 
where the value for the nucleation rate was taken as unity, using a form similar to 
equation 2.6. Calculated values of this energy reduction factor, based on measured 
quasi-steady heater surface temperatures at the onset of boiling, were correlated in 
terms of: 


B 

(2.18) 

for a number of different fluids, including R12, LN 2 , R1 14, L0 2 , R22, ethyl ether, and 

water. A and B were determined by curve fitting calculated values of F for a given 
substance as a function of the dimensionless temperature in equation 2.18. Critical 
heat flux heater surface temperature measurements were used in addition to the quasi- 
steady incipient boiling heater surface temperature measurements, as it was felt that 
the incipient boiling data alone were unreliable. However, this point is not clear. 
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2.2.3 Liquid Temperature Distribution 

The role that the liquid temperature distribution in the vicinity of the heater 
surface has on incipient boiling will now be considered. Incipient boiling models 
which include the liquid temperature distribution will be examined first, followed by a 
discussion of the effect of the mode of heating. 

Early models of initial boiling such as that of Griffith and Wallis (I960) and 
Corty and Foust (1955) assumed a liquid of uniform superheat. The predicted heater 
surface superheat at boiling inception by Griffith and Wallis agreed with the 
measurements only for a surface immersed in a uniformly heated bath. Furthermore, 
the two models referred to above placed no upper limit on the size of a potentially 
active cavity. It would appear reasonable that a cavity size must exist beyond which it 
could no longer become a nucleation site, since a cavity size large enough would be 
completely wetted. 

Models for incipient boiling which attempt to incorporate the influence of 
liquid temperature distribution on the process assume a spherically shaped pre-existing 
vapor nucleus at the entrance of a heater surface cavity, and generally relate the 
difference between the internal vapor pressure of the nucleus and interface liquid 
pressure to the radius of curvature by equation 2.2. As was noted earlier, equation 2.2 
is based on the assumption of a uniform liquid temperature and a nucleus of spherical 
geometry. This pressure difference is then related to the liquid superheat necessary by 
the Clausius-Clapeyron equation, which is based on equilibrium conditions. Analyses 
vary according to the equilibrium nucleus shapes assumed and the criteria necessary 
for the nucleus to grow, including the liquid temperature distribution. The growth 
criterion is usually that the local liquid superheat must attain a value greater than the 
saturation temperature of the vapor nucleus at some defined location. 
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The model of Hsu (1962) predicts a size range of potentially active heater 
surface cavities and will be described in some detail here since it incorporates the 
critical features of most nucleation models in varying temperature fields, and was one 
of the earliest to be presented. This model assumes the existence of a vapor nucleus of 
spherical cross-section geometry sitting at the mouth of a surface cavity. The origin of 
the nucleus is unclear. The liquid in a slab of thickness 8 adjacent to the heater surface 
is taken to be heated in a transient conduction process by the solid heater surface. It 
was assumed that the temperature of the liquid at a distance greater than 5 from the 
heater surface remained at Too. Figure 2.2a presents the coordinate system used and 

defined quantities. The model uses the following form of the energy equation: 

80 820 

< 119 > 

where 0 = T - T^. The boundary and initial conditions for the constant temperature 
case of a heater surface are: 


0(x,O) = 0 

(2.20) 

0(0, t) = 0 

(2.21) 
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The well-known solution is: 
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(2.23) 
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For a nucleus of height b above the heated surface, the nucleus is presumed to grow 
when 0 equals 0^ at x = 5-b. 0^ results from the Clausius-CIapeyron relationship and 

equation 2.2: 


2oT 


0 . = 0 + 

b s 


h r p r 
fg K v n 


(2.24) 


r n is an effective bubble radius to be defined below. In this discussion, the saturation 

pressure corresponding to a flat liquid-vapor interface will be referred to as the 
"normal saturation pressure." T s in equation 2.24 represents the saturation temperature 

in a system at a pressure corresponding to a flat interface and is similarly termed the 
"normal saturation temperature." 0 S = T s - T*, and represents the bulk liquid 

subcooling. For a saturated liquid 0 S = 0. 

A relation among the bubble height, b, the radius of the bubble nucleus, r n , and 
the cavity radius, r c , was determined assuming that the vapor bubble was in the shape 
of a truncated sphere: 
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substituting equation 2.25 for r n in equation 2.24: 
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where c was determined from the assumed geometry of the truncated sphere. 5 was 
measured experimentally or estimated, as in Wiebe and Judd (1971). Dividing 
equation 2.26 by: 


9 W — T w - Too 


(2.27) 


yielded: 
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(2.28) 


When the liquid temperature at xj, is equal to 0|j in equation 2.28, the bubble nucleus 

was said to begin to grow. For the transient conduction process within the liquid layer 
of thickness 5, this implies a tangency between the liquid temperature distribution 
curve and the equilibrium bubble temperature curve. Figures 2.2b through 2.2d 
present the sequence of events necessary for the bubble nucleus to grow as described 
above. Figure 2.2b presents the liquid temperature distribution before heating. Figure 
2.2c shows the liquid temperature distribution during heating, but the liquid 
temperature distribution depicted does not yet pass through the equilibrium bubble 
temperature. Figure 2.2d presents the liquid temperature distribution passing through 
the equilibrium temperature, and the nucleus of equivalent radius r n is now assumed to 

grow. 

Cavities which required an infinite time to become active were considered to 
be inactive. At t = «>, the temperature distribution in the liquid, from equation 2.23, is 
linear: 
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When equations 2.28 and 2.29 are set equal to each other at x = 5-b, the upper and 
lower limits for a range of cavities which can be active may be found: 




(2.31) 


where: 


c, = (l+cos<|))/sin<|) 


(2.32) 


and 


c 3 = l+cos<t> 


(2.33) 
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Hsu used the values of 1.6 and 2 for the constants c 3 and c,, respectively. The model 
of Hsu described above gives a range of cavities on a constant temperature boiling 
surface which are potential sites for boiling. For a saturated liquid: 


e s =o 



For saturation conditions, it is required that: 
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(2.34) 


(2.35) 
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(2.37) 


for a physically reasonable cavity size range. With this in mind, Hsu gave a criterion 
for the boiling inception in terms of a critical 0 WO for the case with an imposed heater 

surface temperature: 
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(2.38) 


Boiling incipience occurred if the imposed heater surface temperature was greater than 
0 WO . If 8 were known by some measure, the surface temperature for incipient boiling 

could be predicted for a given subcooling from equation 2.38. 

The case for a constant heater surface heat flux was also analyzed. The range 

of potentially active cavities for this situation was expressed by equations 2.30 and 
2.31 by replacing 0 W with q"8/k. A criterion for the boiling inception in terms of a 
critical q" Q was also given for the case with an imposed heat flux: 
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(2.39) 


For boiling incipience q" > q" 0 , where 8 again must be determined in some manner in 
order to calculate q" 0 . 

A disadvantage associated with the Hsu model for incipient boiling is that the 
constants used in the model, including 8, cannot be measured, and must be assumed 
empirically to provide the best fit of the experimental data. Wiebe and Judd (1971) 
found that measured values of the heating surface temperature at boiling incipience 
agreed well with the 0 WO of equation 2.38. Shoukri and Judd (1975) found that 

although equation 2.38 predicted a surface superheat of 0.25 °C for the onset of 
boiling, no boiling was observed until a surface superheat of 5 °C was attained. 
Shoukri and Judd felt that boiling inception did not occur precisely at 0 WO because of 
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the nonexistence of sufficiently large cavities with the prerequisite trapped gas. They 
also found that measured boiling cavity radii data fell in the lower portion of the range 
of cavity radii predicted by the Hsu model. 

A conduction model for boiling inception similar to that of Hsu is the model of 
Han and Griffith (1965). In this model the physical situation was formulated as a 
transient conduction problem in a semi-infinite liquid, where the temperature 
distribution has been solved by the Laplace transform method. The spatial derivative 
of the liquid temperature evaluated at the heating surface was set equal to the slope of 
an assumed linear temperature distribution in the liquid. The transient thermal 
boundary layer thickness was defined as: 



(2.40) 


The time from the onset of heating of the liquid to boiling initiation was defined as the 
waiting period, t w . During the waiting period, the bubble nucleus was assumed to be 
hemispherical. The liquid temperature at a distance of 1.5 r c from the heating surface, 

located on the interface of the vapor nucleus, was equated to the vapor bubble 
equilibrium temperature, yielding: 
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The distance 1.5 r c , rather than a presumed cavity radius of r c , comes from an analogy 
with potential flow, but this analogy is not clear. The effective cavity range was found 



27 


in a manner similar to that of Hsu (1962) and defined in terms of a maximum thermal 
boundary layer: 
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c max 
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c min 



(2.42) 


(2.43) 


For T w = T s , r cmin and r cmax are both zero, which means that boiling does not occur. 

The maximum thermal boundary layer thickness was considered to be capable 

of being determined from convection correlations. A disadvantage of the Han and 
Griffith model for the onset of boiling is the use of empirical constants, including r c or 

5. The heater wall temperatures computed from the Han and Griffith model were 
compared to the measured maximum transient wall temperatures in the work of Hall 
and Hamson (1966). The measured maximum wall temperatures were 30-40 °C 
higher than the computed wall temperatures at boiling incipience. Sakurai and Shiotsu 
(1977) claimed that Hall and Harrison should not have compared the measured 
maximum transient surface temperature with the incipient surface temperature as 
calculated by the Han and Griffith model, but rather felt that the measured incipient 
boiling surface temperature, which is not necessarily the maximum transient surface 
temperature, would more closely follow the Han and Griffith model. Ibrahim and 
Judd (1985) found that the Han and Griffith model does not correctly predict the 
bubble waiting time at high degrees of subcooling. 



28 


The Bergles and Rohsenow (1963) model of boiling initiation assumes a linear 
liquid temperature distribution near the heating surface, and relies on convection 
correlations for the determination of the active cavity range. Rohsenow (1970) later 
indicated that the trend predicted by this model could not be observed clearly in a large 
number of experimental data. 

Sakurai and Shiotsu (1977) used a conduction model similar to those of Hsu, 
Han and Griffith, and Rohsenow, with different empirical constants to describe the 
geometry of the pre-existing nucleus. Experiments using an imposed exponential rise 
in heat flux with a horizontal platinum wire verified the predictions of their model. 

The above conduction models neglect any turbulence effects which may be 
present in a flow field near a heating surface, and may explain the large errors 
observed when predicting the heater surface temperature at the onset of boiling. 

2.2.4 Externally Forced Liquid Motion 

The temperature distribution in the liquid adjacent to the heated surface is 
known to affect boiling inception. For example, the liquid temperature distribution 
produced by rapid transient heating is quite different from that as a consequence of 
quasi-steady heating. Homogeneous nucleation of various organic liquids initially at 
atmospheric pressure has been produced by heating a platinum wire immersed in the 
test fluids with an intense current pulse of short duration, in the studies of Simpson 
and Walls (1965-1966). However, what appears to be heterogeneous nucleation only 
has been obtained by quasi-steady heating from a solid surface immersed in organic 
fluids. 

Externally forced liquid motion, as in forced convection boiling systems, 
influences the liquid temperature distribution. In a gravity field, the degree of this 
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influence is determined by the magnitude of the ratio of the buoyancy forces on the 
liquid to the inertial forces, and is expressed as the Richardson number: 

Ri = -^ (2.44) 

Re 

Ri = 1 : corresponds to "mixed" convection 
Ri » 1: corresponds to "natural" convection 
Ri « 1: corresponds to "forced" convection 

These parameters are defined in the List of Nomenclature. Figure 2.3 shows a solid 
surface with an imposed step in heat flux. Two mean liquid temperature distributions 
are presented: one mean temperature distribution in a flowing liquid and the other in a 
quiescent liquid. The effect of fluid velocity, regardless of the specific forces inducing 
the fluid motion, produces a temperature distribution which is lower than that of the 
liquid with zero velocity. The heater surface temperature at incipient boiling and the 
ensuing phase change dynamics are surmised to be different for the two temperature 
distributions. 

Incipient boiling with forced convection and with transient heating of the 
surface is defined here for operational purposes as the appearance of the first visible 
bubble anywhere on the surface. Incipient boiling with forced convection as used here 
is distinctly different from what is generally termed incipient flow boiling, which 
usually refers to the quasi-steady heating of a flowing liquid until a vapor bubble 
appears either as the initial one on the heated surface or as an additional one at a 
location where no previous boiling had taken place. The latter will be referred to here 
as quasi-steady incipient boiling. Quasi-steady incipient flow boiling studies are 
reviewed in Yin and Abdelmessih (1977), Hino and Ueda (1985), and Lee et al (1988) 
and are not considered in the present work. 
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Few works have been reported on transient forced convection boiling, with 
even fewer reported for transient forced convection incipient boiling. Fabic (1964) 
used 4 mil thick nickel-iron (Deltamax) ribbons in flowing water to measure transient 
mean surface temperatures and measured the time from when the ribbon was energized 
to the appearance of the first bubble. The power input to the ribbon increased 
exponentially with time. Measurements of the mean surface temperatures by means of 
reading ribbon voltages from Polaroids of an oscilloscope had large uncertainties, on 
the order of + 12 to 19 °C. Fabic presented a model based on a uniform temperature 
for incipient boiling which relied on assumed values of the advancing and receding 
contact angle. He also found that the heat flux predicted by the Bergles and Rohsenow 
(1963) correlation for incipient boiling with quasi-steady heating; 


f u 0 . 0234 . . 
(2. 30 /(p | )) 

q"j = 15.6 P| ' 56 (Tw-Ts) 


(2.45) 


did not agree with the transient measured heat flux at the onset of boiling. 

2.2.5 Buoyancy Induced Liquid Motion 

When the local liquid temperature near the heating surface exceeds the 
saturation temperature by some amount, determined by the microgeometry of the solid 
surface, the surface-liquid properties, the surface temperature of the solid, and the 
temperature distribution in the liquid, incipient boiling occurs. The heater surface 
temperature and the liquid temperature distribution, in turn, depend on the heat flux 
and velocity distribution of the liquid as affected by gravity induced buoyancy. For 
forced convection boiling with Ri « 1, then the buoyancy forces are negligible. If Ri 
= 1, then the buoyancy forces and the inertial forces acting are of similar importance. 



31 


Forced convection studies involving boiling inception were described above in the 
previous section. 

For pool boiling in earth gravity, the liquid is initially stagnant. Buoyancy is 
the only force to initiate bulk liquid motion during the heating of the liquid with Ri » 
1. The onset of single phase natural convection is an instability phenomenon induced 
by buoyancy forces. Studies with the onset of natural convection occurring before 
boiling inception in a heated liquid pool at earth gravity or greater acceleration are 
reviewed in Ulucakli (1987). Additional works concerning the transient heating of a 
liquid pool in earth gravity include that of Oker and Merte (1978) and Nghiem (1980). 
In Oker and Merte (1978), the departure of the measured transient gold film 
temperature from the one dimensional conduction prediction was considered to be the 
onset of natural convection. For a flat rectangular heater surface horizontal up in earth 
gravity, a semi-empirical correlation for the onset of natural convection was found and 
is referred to in Chapter 5. In Nghiem (1980), transient boiling in water and Rl 13 
with various heater surfaces was investigated experimentally with an imposed step in 
heat flux. Prior to boiling, the onset of natural convection was observed with 
holographic interferometry for heated wires and a narrow (20 mm x 2 mm) thin gold 
film heater. Nghiem proposed an empirical heterogeneous nucleation factor which 
also included natural convection effects, and was described above in section 2.2.2. 

2.3 Relationship to Proposed Research 

The influences of surface microgeometry, surface energy, liquid temperature 
distribution, externally forced liquid motion, and buoyancy induced liquid motion in 
previous studies have been presented above in section 2.2. Information lacking in 
previous studies but pertinent to the objectives of this work is now discussed. 
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In Oker (1973), it was shown that fluid velocity as a consequence of gravity 
induced buoyancy can affect initial boiling and the early motion of the newly formed 
vapor phase. In Figure 2.4, incipient boiling with transient heating occurred at a 
considerably larger surface temperature for the horizontal up orientation than with 
microgravity, with all other test conditions held constant. This difference in surface 
temperature at initial boiling demonstrates that the temperature distribution in the 
liquid, which is effected by the fluid velocity, influences incipient boiling. The effect 
of a known liquid temperature distribution at the moment of boiling inception and the 
resulting phase change dynamics has not been previously examined and is an objective 
of the present study. 

In transient pool boiling in the absence of buoyancy, an initially quiescent fluid 
remains stagnant upon heating until the onset of boiling. The transient temperature 
distribution at the moment of initial boiling in a motionless fluid can be determined 
from a conduction heat transfer analysis with a high degree of certainty. A heated 
liquid in microgravity does not have buoyancy induced fluid velocities, and thus, 
microgravity provides a known liquid temperature distribution at boiling inception. 

Microgravity conditions may be obtained with drop towers, aircraft in 
parabolic flight trajectories, and sounding rockets. Although the available time in drop 
towers is short, up to a maximum of five seconds, a high vacuum drop environment 
can provide high quality microgravity, down to a/g = 10* 5 . The time available for 
reduced gravity in aircraft flying a parabolic trajectory is typically greater than the 
time of microgravity in drop towers, up to 25 seconds. The quality of the reduced 
gravity when compared to microgravity conditions of a drop tower is considerably 
less, on the order of a/g = 0.02. In addition, the initial quiescent states required to 
obtain a well defined temperature distribution for studies of incipient boiling are not 
possible since operation requires large gravity fields preceding each trajectory, in 
addition to the normal disturbances associated with aircraft in flight. Sounding rockets 
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have been used to achieve high quality microgravity, claimed to be on the order of a/g 
= 10" 4 for up to 140 seconds. The small size inherent in such vehicles requires 
limitations in the size of the heating surface and in the ability to obtain precise 
measurements. Repeated rocket launches are costly, and the possibility for non- 
reproducibility of the experimental conditions between tests is increased, particularly 
if new test equipment must be fabricated because of damage during retrieval. The 
present study is concerned with incipient boiling, and Merte and Litdes (1971) 
presented a summary of incipient boiling work in microgravity. 

Table 2.1 presents a survey of transient incipient boiling experimental studies. 
Oker (1973) examined the results of initial boiling of saturated R113 and LN2 in a 
drop tower (a/g = 0.004 for 1.4 seconds) using a semi-transparent gold film on pyrex 
glass. Oker was the first to measure the transient heater surface superheat required for 
boiling inception in microgravity and in earth gravity for various orientations. The 
degree of initial temperature uniformity was unknown, and special consideration was 
not taken to absorb the infrared radiation from the 200 watt mercury arc lamp used for 
the high speed filming. In addition, the system pressure was not controlled. 

Straub et al (1990) summarized the results of a German microgravity boiling 
test program from the mid 1970s to the present which used sounding rockets and 
parabolic aircraft flights. Many of the tests were performed using heated wires. The 
diameters of these wires were of the same size as the thermal boundary layer 
surrounding the heated wires and comparable to the size of the vapor bubbles which 
nucleated on them. Surface tension effects on vapor motion were quite different for 
wires as compared to flat surfaces. Often the wires were completely enveloped by 
vapor at boiling inception, and subsequent surface tension effects maintained a 
"pseudo film” boiling only because of the wire geometry. Straub et al concluded that 
the effect of buoyancy on pool boiling was very small. This conclusion, however, was 
based largely on the results of the small diameter wire experiments (Weinzierl (1982), 
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Zell et al (1984)). In addition to the small heated wires, Straub et al also used thin 
film resistance surface heaters for boiling tests. 

In transient heating, different boundary disturbances result in different 
temperature distributions. An imposed step increase in heat flux has a constant 
temperature gradient at the heater surface and provides an elementary boundary 
disturbance* all other boundary disturbances can be constructed from a series of steps 
in heat flux. An imposed heater surface temperature is another elementary boundary 
disturbance important in heat transfer and produces a temperature distribuuon in the 
liquid quite different from that of a step increase in heat flux. In Table 2. 1, it may be 
noted that no previous study used an imposed surface temperature to study transient 
incipient boiling, as was used in the work described here. 

The present work is the first experimental study of incipient boiling with 
transient heating in microgravity where the quiescent fluid had an initial known 
uniform temperature, and the pressure of the system was precisely controlled for 
various subcoolings. The microgravity tests involved a single step in heat flux 
imposed on a large flat heating surface for a five second duration with a/g = lfr 5 . As a 
result of the well defined initial conditions of the tests and the single step in heat flux, 
the temperature distribution in the liquid near the heating surface was computed with 
certainty not previously possible. The present study is also the first to utilize an 
imposed surface temperature to specifically examine incipient boiling. 
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Figure 2.1 


Gas bubble at a solid-liquid boundary 
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Figure 2.2a Pre-existing vapor nucleus protruding from a cavity 



Figure 2.2b Initial liquid temperature distribution 



Figure 2.2c Liquid temperature distribution after beating, but below 
the equilibrium bubble temperature, equation 2.28 



Figure 2. 2d Liquid temperature distribution tangent to the 

equilibrium bubble temperature, and nucleus now 
begins to grow 
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Figure 2.3 Influence of convection on mean liquid temperature 
distributions 
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CHAPTER 3 

EXPERIMENTAL APPARATUS 

This is an experimental study of incipient boiling with variable subcooling for 
pool boiling. Pool boiling experiments were conducted in both earth gravity and 
microgravity, Two different test vessels were used: a pool boiling vessel for earth 
gravity studies (designated as the laboratory pool boiling vessel) and a pool boiling 
vessel for microgravity tests in a drop tower (designated as the drop vessel). The heat 
flux, heater surface temperature, pressure near the heating surface, and liquid bulk 
temperature were measured in both test vessels. High speed video recordings for the 
earth gravity experiments and high speed black and white movies for the microgravity 
tests were obtained. The videos and movies were recorded using two orthogonal 
views, across and through each test surface. 

3.1 Common Experimental Elements 

The heat transfer test surfaces, heater surface power controllers, test fluid, data 
acquisition systems, thermocouples, optics, pressure control system, camera 
synchronization, and vessel power supply were common to both of the test vessels and 
will be discussed before the unique details of each of the individual vessels are given. 
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3.1.1 Heat Transfer Test Surfaces 

Two different kinds of heating surfaces were used in the experiments: a thin 
gold film deposited on a quartz substrate and a massive gold coated copper surface. 
Two types of boiling tests were performed in this study, transient and quasi-steady. 
The gold film heaters were used in both the transient and the quasi-steady tests, while 
the large thermal mass of the copper heater restricted it to the quasi-steady heating 
tests. 

Figure 3.1 presents the concept of the gold film coated quartz heating surface. 
A 400 angstrom thick gold film is sputtered on a quartz substrate which was pre- 
sputtered with a 30 angstrom tantalum layer for improved gold adhesion. Before the 
tantalum layer was applied, the surface was cleaned by ion bombardment to remove 
about 10 angstroms of the quartz substrate. The thin film served simultaneously as a 
heater and a resistance thermometer, and had a negligible time constant associated 
with the transient temperature measurement. The temperature of the thin film thus is 
identical to that of the surface of the quartz substrate. The quartz substrate was highly 
polished by the manufacturer, using a 1.4 micron pitch polish. The semi-transparent 
gold surface consisted of a rectangular shape 1.91 cm by 3.81 cm (0.75 inches by 1.5 
inches), and was expected to be considerably larger than the size of the bubbles to be 
formed above it. This semi-transparent gold surface had a higher resistance than the 
10,000 angstrom thick gold surface near the power taps, which acted as the current 
connection. The voltage measuring leads also consisted of 10,000 angstrom thick gold 
films. Kovar screws passing through holes drilled in the quartz were used as power 
and voltage taps. The low linear thermal expansion coefficient of Kovar closely 
matched the linear thermal expansion coefficient of the quartz substrate, which was 
mounted in a Teflon holder. 
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As in the work of Oker (1973, 1981), it was found that the electrical resistance 
of the gold surfaces varied with time, but that the slope of the linear relationship of 
surface resistance and temperature, dR/dT, for each particular surface remained 
constant. This resistance change was attributed to an aging effect. Following the 
suggestion of Oker (1973), the surfaces were tempered at a temperature of about 275 
°C for 20 minutes, which accelerated the aging to a stable level. It was also found in 
early testing that sputtered gold surfaces were more durable than those deposited by 
vacuum evaporation. The temperature-electrical resistance relation for each thin gold 
film surface was determined by calibration, and the details are presented in Appendix 
B.l. On the day an experiment was to be performed, the thin film heater was 
calibrated in place within the test vessel before heating the vessel and the test fluid to 
the test temperature level. This calibration is referred to here as the single point 
calibration. The single point calibration consisted of passing a small current, on the 
order of 20 milli-amps, across the thin film heater to produce negligible ohmic heating 
while simultaneously measuring the differential voltage across the heater surface and 
the differential voltage across a heater shunt, which was in series. With the 
differential voltage across the heater surface and the differential voltage across the 
shunt known, together with the calibrated shunt resistance, the heater surface 
resistance was determined from Ohm's relationship. Since the test vessel and fluid 
remained at the ambient temperature, the temperature of the heater surface during the 
resistance measurement was known from the calibrated test vessel thermocouples. 
Because the slope of the linear temperature-resistance relationship remained constant 
from previous calibration, this single measured heater surface resistance and 
corresponding measured fluid temperature determined the constant required for the 
temperature-resistance relationship. Additional details of the single point calibration 
are given in section 3.1.4. 
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Figure 3.2 shows the design of the typical metal heater used for the study of 
initial boiling under quasi-steady heating. The heater consisted of a copper block with 
a film heater imbedded in a Stycast electrical encapsulating potting compound. A 
copper foil, 0.00254 cm (0.001 inch) thick was soldered to the copper block and to the 
stainless steel holder to expose a continuous surface to the heated fluid, thereby 
preventing preferential boiling at the crevice that would otherwise be present. The 
copper foil was sputtered with gold to provide the identical heater surface-liquid 
energy combination as the gold deposited on the quartz. A chromel-constantan 
thermocouple with a bead of 0.025 inches was inserted into a 0.030 inch diameter hole 
in the copper heater in order to determine the temperature of the copper foil. The hole 
for the thermocouple was filled with a thermally conducting copper compound. The 
tip of the thermocouple was about 0.030 inches from the gold sputtered copper foil 
and was essentially at the same temperature as the heating surface. To prevent any air 
trapped between the stycast and the foil from expanding and deforming the foil during 
the period when a vacuum was imposed on the upper surface (as required for gold 
deposition and for the charging of the test vessel with the test fluid), it was necessary 
to embed an evacuation tube passing through the stycast to a point just under the 
copper foil. The other end terminated in a "quick-disconnect" fitting which could be 
attached to a vacuum pump when required. 

3.1.2 Heater Surface Power Controllers 

Two modes of heating were used in this study, transient and quasi-steady. For 
the transient heating studies two different conditions were imposed on the heater 
surface: a step increase in heat flux and a step increase in temperature. A heater 
surface controller was designed for use with both the step increase in heat flux and the 
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quasi-steady heating, and a different heater surface controller was designed to provide 
the step increase in surface temperature. 

The power supply for the step increase in heat flux and the quasi-steady 
heating maintained a nominally constant voltage across the thin film surface for the 
boiling tests. Since the resistance of the thin film test surfaces during heaung changed 
on the order of a few percent, the product of the current and voltage remained 
essentially constant. The constant voltage power supply could accept either a 
microprocessor input or a manual voltage pot signal as an input, and provided an 
output signal to the base of the power transistors, which in turn powered the thin film 
surface heater or copper block heater. The heater surface voltage was used as a 
feedback signal, compared with the set point voltage, and the difference amplified to 
furnish the appropriate signal to the power transistors. Details of the heater surface 
controller circuit for the step in heat flux and quasi-steady heating are given in 
Appendix E.l. 

The imposed step increase in surface temperature tests, performed only in the 
laboratory pool boiling vessel, required a power supply different than that for the 
constant voltage heater surface controller. With the heater surface inverted in earth 
gravity for heating periods on the order of milliseconds, conduction heat transfer from 
the heater surface to the liquid may be used to determine the transient heat flux 
required to impose a constant heater surface temperature. The solution of equation 
2.19, the one dimensional transient conduction equation, together with the boundary 
condition of an imposed step in heater surface temperature for two semi-infinite media 
with heat generation at their common plane, is well known. The result is a two 
domain solution, a temperature distribution for each of the two media. The partial 
derivative of the solution for each domain was evaluated at the common plane and 
multiplied by the respective thermal conductivity to obtain the heat flux from the plane 
of energy generation to each of the semi-infinite media. The transient heat flux to 
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each of the two domains was summed to obtain the required total energy input to the 
thin gold film heater: 




T -T 



(3.1) 


The details of the above analysis resulting in equation 3.1 are given in Appendix F. In 
the present experiment, a series of 24 imposed steps in heat flux were used as an 
approximation to equation 3.1. An illustration of this approximation to the exact 
expression for the total imposed heat flux is shown in Figures 3.3 and 3.4. The total 
heat flux-time curve initiated at a pre-set maximum total heat flux for 400 
microseconds, then decreased in 23 discrete steps ranging in time intervals from 22 
microseconds to 1.4 seconds, so that equation 3.1 was approximately followed. 

A set of nichrome wire dropping resistors was fabricated to vary the voltage 
across the heating surface. The initial connection with the thin film heater was a direct 
connection to the power supply, providing the maximum power to the heater. The 
subsequent resistance steps reduced the power to the heating surface according to the 
time schedule given in Table 3.1, which apply for an initial maximum voltage. Each 
successive power level was produced by a different resistor and switch in series with 
the heater. One side of each resistor was connected to the supply voltage, and the 
other side was connected to a switch which was in series with the heater. Figure 3.5 
illustrates the arrangement of the switches, resistors, electronic controls, and power 
generator. 

The times for switching each resistor to approximate equation 3.1 were used in 
a microprocessor program written in assembly language to provide the switching 
signals, at the appropriate times, to the switching circuits. The switching currents 
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ranged from 37.5 amps to about one amp. Twenty-four MOSFETs were used as the 
solid state switching devices with the capability to perform the actual switching 
process in less than one microsecond. The MOSFETs were actuated by the output port 
controlled by the microprocessor, using timing pulses. (A MOSFET is a metal oxide 
semi-conductor field effect transistor.) Details of the circuitry for the imposed step in 
surface temperature are given in Appendix E.2. 

3.1.3 Test Fluid 

Commercial grade R113 (Refrigerant 113 or trichlorotrifluoroethane, 
CC1 2 FCC1F 2 ) was used for all of the experiments. This fluid, with a normal boiling 

point of 47.6 °C, was chosen as the test fluid because of its well established 
thermodynamic properties. The low normal boiling point permitted variation of 
subcooling by the control of system pressure. R113 is also electrically non- 
conductive, which prevented "deplating" of the thin gold film heater as would occur if 
water were used. Because of its vapor pressure and surface tension characteristics, 
R1 13 is representative of "space fluids" such as cryogens and other fluorocarbons. 

A degassing and purification unit was necessary because of the high solubility 
of air in R113. For example at 25 °C and 1 atm, R113 saturated with air has a mole 
fraction of air dissolved nearly equal to 1.15 x 10 3 or a volume fraction of air 
dissolved of about 0.23 Dupont (1990). The degassing unit was fabricated of 304 SS, 
and a sketch of the system is shown in Figure 3.6. The concept of the degassing 
system involved a combination of distillation at room temperature, leaving behind the 
high boiling point components such as oils and solids, and freezing the R113 on the 
fins by LN 2 within the inner vessel. A sufficiently low pressure was maintained such 

that the air components, except for water vapor, were removed by a vacuum pump. 
Davison Chemical molecular sieves made from synthetic silico aluminate zeolite, with 
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a four angstrom effective pore size, were used to remove the water vapor, and also 
served as a mechanical filter. 

3.1.4 Data Acquisition Systems 

Because of the large differences required for the rate of heating in the imposed 
constant heat flux experiments and the imposed surface temperature experiments, two 
different voltage measuring systems were used for data acquisition: the CR7X and the 
Codas system. 

The Campbell Scientific Co. CR7X (Serial No. 1448) is a programmable 
digital measuring and controlling device used to measure differential voltages from the 
thermocouples and the pressure transducer and across the heater surface and the heater 
shunt, in addition to sequencing events involving the drop vessel. The CR7X could 
record data within eight voltage intervals from +1.5 mV (with a resolution of 50 
nanovolts) to ± 5 V (with a resolution of 166 microvolts) at two different speeds, 
designated as fast or slow. The CR7X sampled the voltage input channels sequentially 
in time. 

Four tests were performed to determine the performance characteristics of the 
CR7X under various measuring conditions. The results of these tests are summarized 
here, with the details given in Appendix C. The maximum uncertainty in voltage on 
the fast sampling rate mode, the mode in which all tests were run, was ± 0.6 
microvolts. This maximum uncertainty in voltage readings corresponded to an 
uncertainty in temperature as indicated by the thermocouples of about + 0.01 °C for 
the 15, 50 and 150 mv ranges at the bulk liquid temperature level of 49 °C. The 15 
mV and 50 mV voltage measurement ranges were the most desirable in terms of 
stability and accuracy in the 20 to 60 °C ambient temperature range, and had the 
smallest errors when compared to a standard reading from a potentiometer at the slow 
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sampling rate. It was determined that an ice reference junction for thermocouples 
provided more accurate temperature measurements than when the internal reference 
junction compensation of the CR7X was used. The dynamic measurements of the 
heater shunt and the heater surface voltages as measured by both the CR7X and the 
7700 series 8 channel Sanborn recorder, which was used as a standard, were 
compared. The CR7X dynamic voltage measurements agreed with the Sanborn 
recorder measurements. 

A schematic of the voltage measurement circuit for the determination of the 
surface temperature of the thin gold film is shown in Figure 3.7. The heater power 
supply set the voltage level for the circuit. Voltage differentials across the shunt and 
surface were measured with the CR7X in all of the tests except for the experiments 
involving the imposed step in surface temperature, where the Codas system was used. 
A calibrated shunt, 0.016314 + 0.000002 12, was used for the determination of the 
current in tests in the laboratory pool boiling vessel. For the single point calibration 
in all of the laboratory tests, the National Bureau of Standards calibrated standard 
resistor (Leeds and Northrup Model No. 4221-B), with a resistance of 0.1000008 ± 
0.000001 £2 was used. In the drop vessel, a relay was used to switch voltage 
difference measurement from a calibrated shunt, 0.014788 ± 0.000002 £2, used for the 
drop tests to a calibrated resistor, 3.770 + 0.001 f2, used for the single point surface 
calibration. Figure 3.8 presents a schematic of the drop vessel heater surface 
experimental run and calibration circuit. 

In order for the CR7X to sample in the fast mode, it was desirable to have 
voltages from the heater surface divided proportionately so as to be in the same 
voltage range as the shunt (150 mV). The heater surface voltages in all of the tests 
were divided by stable calibrated resistors. 

For the imposed surface temperature experiments, the Codas system was used 
to measure differential heater and shunt voltages, in addition to the CR7X, because the 
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initial portion of the power transient was too rapid for the CR7X to sample. The 
Codas system consists of the software on an IBM PS/2-30 286 for acquiring data, an 
analog and digital interface card, and 16 single-ended (8 differential) analog input 
channels. The analog-to-digital converter has a 12-bit resolution, 1 part in 4096, with 
a conversion time of 20 microseconds (including acquisition) and a maximum 
throughput rate of 25 kHz for two channels. The full scale voltage input range of the 
Codas was -2.5 to +2.5 volts. The Codas software permitted programmable data 
acquisition and waveform analysis (integration, differentiation, rectification, and 
filtering) of the acquired data. The uncertainty in the Codas voltage measurements is 
± 1.0 millivolts, presented in Appendix C. 

3.1.5 Thermocouples 

Chromel-constantan teflon sheathed thermocouples with ungrounded junctions 
were used for the measurement of the R 1 1 3 temperature. These thermocouples were 
fabricated from a 30 gage thermocouple spool from Omega Engineering Co. The 
measuring junctions of the thermocouples were carefully welded to produce a 0.051 + 
0.013 cm (0.020 ± 0.005 inch) bead. The thermocouple reference junction was 
maintained in a dewar containing crushed ice and water, and each of the thermocouple 
wires and extension wires were numbered for identification purposes. The dewar for 
the drop vessel had an o-ring seal to avoid leakage under the external vacuum 
conditions. Copper wires were used as extension wires from the ice junction to the 
CR7X terminals. The thermocouple spool was calibrated with a platinum resistance 
thermometer, and the details of this calibration are given in Appendix B.2 The final 
accuracy of the temperature measurement as determined by a calibrated thermocouple 
was on the order of ± 0.01 °C, as discussed in Appendix D.2. 
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3.1.6 Optics 

Two views of the heating surface were recorded in all of the transient 
experiments for purposes of determining the incipient boiling time and for observation 
of phase change dynamics. No images of the heating surface were recorded for the 
quasi-steady heating tests. 

The optical arrangement for the laboratory pool boiling vessel and the drop 
vessel were essentially the same. Figure 3.9 shows the optical arrangement for the two 
orthogonal views of the heating surface. Light passed across the heater surface and 
simultaneously through the surface for the two simultaneous views. The light 
sources, labelled LI and L2, were 360 watt bulbs powered by a 110 volt a.c. 
transformer to set the desired light intensity in the case of the laboratory pool boiling 
vessel. For the drop vessel, LI and L2 were 80 watt lights. It was found that neither 
the 360 watt lights, set at their highest intensity in combination with the water cooled 
infrared filter (which contained water as an infrared absorbing media) for use with the 
laboratory pool boiling vessel, nor the 80 watt lights, used in combination with the 
cold mirrors of the drop vessel, produced any detectable increase in R113 temperature 
as measured by the vessel thermocouples over a 3 minute period (a much longer time 
period than any of the transient tests). Milk glass or ground glass was used to diffuse 
the light to permit improved imaging. Three mirrors, M2, M3, and M5, and a prism, 
PI, directed the light to the high speed camera. A light block reduced the amount of 
light entering the camera across the heater so as to produce a picture of the same 
brightness as the image formed from the light which passed through the surface. The 
mirrors and prisms, Fixed to a vessel-mounted frame, were adjusted so that the two 
light paths were nearly equal. 

The laboratory pool boiling vessel and the drop vessel had different imaging 
cameras. The high speed camera for the laboratory pool boiling vessel was a Kodak 
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Ektapro 1000 Image System, which was operated at speeds up to 1000 frames per 
second. This high speed camera system made video recordings which were stored on 
magnetic tape cartridges. A D.B. Milliken high speed camera (Model DBM 3A RFI, 
28 volts d.c.) recorded boiling phenomena in the drop vessel on black and white film 
at 400 pictures per second with timing marks at 100 Hz for reference. 

3.1.7 Pressure Control System 

A schematic of the pressure control system is shown in Figure 3.10, and was 
essentially the same for both test vessels. The pressure system actively maintained and 
measured the pressure level of the liquid within each of the test vessels, and consisted 
of four main parts: a bellows, valves, electronic control circuits, and a pressure 

transducer. 

A seven inch diameter stainless steel bellows isolated the test fluid, located on 
the inner side of the bellows, from the air (or nitrogen in the case of the drop vessel) 
on the external side of the bellows. The valves, which opened or closed as directed by 
signals from the electronic control circuits, admitted high pressure air to the bellows to 
increase the system pressure or released air from the bellows to lower the system 
pressure. The valves included two (four in the case of the laboratory pool boiling 
vessel) pneumatic solenoid valves (Kip Inc. No. 0141013, 12 volts d.c.), a proportional 
valve (Skinner valve No. BP2EV0035, 24 volts d.c. supply, 0-5 volts d.c. control), and 
a needle valve. 

A description of the electronic control circuits for controlling the pressure is 
presented here, while the circuit diagrams are given in Appendix E.3. The pressure 
control card received a set point either from a microprocessor, in the case of the drop 
vessel, or from a manual input, for the laboratory pool boiling vessel and. The output 
signal from the pressure transducer is also transmitted to the control card. The set 
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point voltage is compared with the actual pressure transducer signal, and their 
difference, suitably conditioned, used as inputs for actuating the vent and fill solenoid 
valves and the proportional valve. While the proportional valve operated 
continuously, the vent and fill solenoid valves were normally closed. They were 
actuated only when the pressure error exceeded a pre-determined value. The vent and 
fill solenoid valves controlled the large pressure transients while the proportional valve 
controlled small pressure fluctuations. The vent and fill solenoid valve controllers 
opened the valves at a given error level and caused them to remain open until the error 
was reduced to half the original error. The proportional valve controller used a 
combination of proportional and integral controller modes. The proportional mode 
responded to pressure transients while the integral controller mode attempted to drive 
the steady state error to zero. The proportional valve was operated within its linear 
range of operation by adjusting the needle valve on the "fill side" so that as voltage to 
the proportional valve was increased or decreased, the flowrate from the bellows 
changed in direct proportion to the voltage. 

A Heise Model 623 pressure transducer was used for the measurement and the 
pressure control in the laboratory pool boiling vessel. The Heise pressure transducer 
had a maximum pressure rating of 344.74 kPa (50 psia) and a sensitivity of 0.01724 
kPa (0.0025 psi). The Heise pressure transducer provided a pressure measurement 
with an uncertainty of + 0.028 kPa (+ 0.004 psi). The pressure level was controlled to 
within +0.172 kPa ( + 0.025 psi) of the set point pressure, corresponding to ± 0.06 °C 
(±0.1 °F) of the saturation temperature The pressure transducers were mounted 
horizontally at the center of rotation on the laboratory pool boiling vessel so that 
compensation for the change in hydrostatic pressure at the heater surface on test 
package rotation took place requiring only a sign change, provided by the software 
which translated the transducer signal into a pressure measurement. 
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The pressure system of the drop vessel was similar to the laboratory drop 
vessel, except for a few minor differences. The gas side of the bellows in the drop 
vessel was pressurized by an on board high-pressure nitrogen tank supply rather than 
by a building air supply as for the laboratory pool boiling vessel . Instead of a 
calibrated Heise 623 pressure transducer, a calibrated Setra model 270 pressure 
transducer which could withstand a 50 g impact was mounted on the top vessel wall 
and provided pressure determination within the desired uncertainty range of ± 0.172 
kPa (± 0.025 psi). In addition, as a guarantee against the possible detrimental effects 
of the drop tower vacuum, the CR7X was mounted in a pressurized box. A complete 
pneumatic diagram for the drop vessel is shown in Figure 3.1 1. 

A discussion of the uncertainty estimates in the pressure transducer 
measurements is presented in Appendix D.3. For pressure calibration of the Heise and 
Setra pressure transducers, a Ruska Pressure Dead Weight Calibration System (2400 
series dead weight gage) was employed, and the details of this calibration are given in 
Appendix B.3 The Ruska system was also used for calibrating a Heise Bourdon tube 
gage which in turn was used for in-place periodic check calibrations of the pressure 
transducers. 

3.1.8 Camera Synchronization 

The camera was synchronized with the powering of the heater surfaces to 
establish a time reference for the observation of incipient boiling phenomena. The 
description of the synchronization apparatus for each of the vessels will now be 
presented, beginning with the laboratory pool boiling vessel. 

To synchronize the high speed video system with the switching on of the heater 
power for the constant heat flux tests in the laboratory pool boiling vessel, a light 
emitting diode was placed in the field of view of the camera. The light emitting diode 
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was in parallel with a copper knife switch to the 24 volt power source. When this 
knife switch was closed, the quartz heater and diode were simultaneously powered. A 
five millisecond delay exists until the light emitting diode provides sufficient 
illumination for viewing by the camera. As a result of this delay, there was an 
uncertainty of + 6 milliseconds in the time interval between heater powering and 
visible light emitting diode illumination. 

For the imposed step in surface temperature experiments in the laboratory pool 
boiling vessel, a light emitting diode was placed in the field of view of the camera to 
synchronize the camera with the powering of the thin gold film heater. The light 
emitting diode, which was powered by a nine volt battery, was switched on by the 
same signal that turned on the first MOSFET switch of the heater driver. The 
uncertainty in the time interval between the powering of the heater and recording of 
the light emitting diode at 1000 fps is + 1.5 ms. This uncertainty in the time interval is 
smaller than that associated with the imposed step in surface heat flux because the 
voltage switching for the diode was done by electronic components rather than by a 
copper knife switch. 

As in the case of the laboratory pool boiling experiments, a light emitting diode 
was placed in the field of view of the camera for the synchronization of the high speed 
camera with the powering of the thin film heater for the drop vessel experiments. The 
CR7X powered the light emitting diode with a 5 volt signal 2.9 milliseconds after the 
thin film heater was energized. The diode took 5 milliseconds to reach maximum 
intensity. The uncertainty in the time interval between powering of the diode and of 
the heater was ± 7 milliseconds. 
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3.1.9 Vessel Power Supply 

For the imposed step in heat flux tests, with the lab pool boiling vessel an 
essentially noise-free 24 volt d.c. power supply consisting of two 12 volt batteries 
connected in series was used. A Sorenson three phase 240 volt a.c. power converter 
supplied up to 150 d.c. volts to the power heater supply for the step in surface 
temperature experiments. The ± 5 and ± 15 d.c. volts required for the control circuits 
were supplied by a 110 volt a.c. convener. An additional 110 volt a.c. converter 
supplied 12 volts for the operation of the solenoid valves. 

The 24 volt d.c. power supply for the drop vessel heater and pressure controls 
were provided by batteries mounted on the drop bus. The ± 5 and ± 15 volt d.c. power 
supply consisted of a 24 volt d.c. power converter mounted on the heater surface 
controller cards. The 12 volt d.c. supply needed for the operation of the solenoid 
valves was provided by batteries mounted on the drop bus. 

3.2 Laboratory Pool Boiling Vessel 

The laboratory pool boiling experiments were conducted at a/g = ± 1. This 
variation in body forces was achieved by inverting the test vessel from horizontal up 
(a/g = + 1) to horizontal down (a/g = -1). Heat flux, liquid subcooling, surface type, 
and orientation of heating surfaces were the test parameters. 

The laboratory vessel shown in Figure 3.12 contained two heater test surfaces, 
both located on the same surface of the vessel. Trunions located on opposite sides of 
the vessel permitted rotation between the horizontal facing up and the horizontal 
facing down orientations with respect to the earth body force direction. The laboratory 
vessel was vacuum helium leak tested and was found to be leak free for the system 
pressurized with R1 13. 
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The inner part of the laboratory vessel containing the R1 13 was surrounded on 
all sides, except where the heating surface and windows were located, by a circulated 
mixture of water and ethylene glycol (50% ethylene glycol by volume). This mixture 
was used as it had desirable viscosity and heat transfer characteristics. A precision 
constant temperature bath (Neslab Endocal RTE-210) pumped the thermostatically 
controlled water-ethylene glycol mixture through the laboratory vessel. The 
temperature stability of the bath was ± 0.02 °C. The temperature non-uniformity was 
0.08-0.11 °C (0.15-0.20 °F) at the operating level of 49 °C, as determined by twelve 
boiling chamber thermocouples. The circulation system of the laboratory pool boiling 
vessel is shown in Figure 3.13. 

The heating surface and viewing windows were surrounded by thick aluminum 
walls to help provide uniform temperature conditions at the start of a test. The 
viewing windows were made of optical quality quartz, and were enclosed by a glass 
cover to minimize convection losses. Fins were located in the annular space between 
the inner and outer sections of the laboratory vessel to promote heat transfer and to 
guide the flow of the circulating water-ethylene glycol mixture. Fins were also welded 
to the inner walls of the inner vessel to increase the heat transfer to the test fluid. 

3.3 Drop Vessel 

The drop vessel experiments were performed in the evacuated 132 meter drop 
tower at the NASA Lewis Research Center, which provided 5.18 seconds of free fall at 
accelerations of 10' 5 g. All of the experimental hardware associated with the drop 
vessel were loaded on a one meter diameter "bus" for the drop. The drop chamber was 
evacuated by a series of pumps to a final pressure near 10‘2 torr. The bus was 
decelerated in a 6. 1 meter deep cylinder filled with polystyrene pellets. 
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The drop vessel and associated electronics were designed to take a 50 g impact 
at the bottom of the drop tower under the conditions of a vacuum. The drop vessel 
was made of aluminum and of essentially the same design as that of the laboratory 
pool boiling vessel, except that the drop vessel was not double-walled and, 
consequently, did not rely on external ethylene glycol and water circulation to 
maintain a uniform R113 temperature at the start of a test. Temperature uniformity 
was achieved through the use of film heaters mounted externally to the drop vessel. 
The drop vessel was vacuum leak tested with helium and found to be leak free for a 
pressurized system. An overview of the drop vessel is shown in Figure 3.15. 

The desired uniform temperature of the test fluid before the drop was achieved 
by film heaters attached to the external surfaces of the test vessel. A laboratory test to 
determine the amount of temperature stratification in the liquid when using the 
external wall film heaters was performed. After a time period of 5 hours following the 
switching off of the external wall film heaters, the temperature stratification was 0.47 
°C from the bottom surface to the top surface of the insulated vessel. A smaller 
temperature stratification (about 0.2 °C) was expected for the vessel in the drop tower 
before the drop, during the drop tower evacuation period for the following two 
reasons: (1) the dominant mode of heat loss in the evacuated drop tower was 

radiation, which was less than the convective heat loss during the laboratory 
stratification test, (2) the time from the switching off of the drop vessel external film 
heaters was less than 2 hours, as compared to five hours for the laboratory 
stratification test. 

The temperature of the six sides of the vessel were controlled automatically 
and independently: the temperature of the bottom surface was set with respect to the 
ice thermocouple reference junction, while the temperatures of the other five surfaces 
were controlled relative to the bottom surface. A thermopile summed the voltages of 
the individual chromel constantan thermocouples distributed on the vessel wall. Mean 
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values of temperatures for a wall were obtained by the average of the thermocouple 
voltages on the wall. Since each of the vessel walls was heated with respect to the 
bottom wall, a differential thermopile measurement was made using the bottom 
surface as the reference junction. Power was supplied to the wall film heaters by an 
a.c. current heater drive. A schematic of the external film heating system is presented 
in Figure 3.15. Diagrams of the associated circuits appear in Appendix E.5. 
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Figure 3.2 Copper block heater for the quasi-steady tests 
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Table 3.1 Power Schedule for Imposed Surface Temperature 


% Of Initial Power 

Step Number 

Resistance (Ohms) 

Time Step (Sec 

100 

1 

0 

4.47E-04 

95 

2 

0.103 

1.96E-05 

90 

3 

0.216 

2.29E-05 

85 

4 

0.338 

2.71E-05 

80 

5 

0.471 

3.23E-05 

75 

6 

0.618 

3.90E-05 

70 

7 

0.78 

4.76E-05 

65 

8 

0.96 

5.90E-05 

60 

9 

1.163 

7.44E-05 

55 

10 

1.392 

9.57E-05 

50 

11 

1.656 

1.26E-04 

45 

12 

1.961 

1.70E-04 

40 

13 

2.323 

2.38E-04 

35 

14 

2.76 

3.47E-04 

30 

15 

3.301 

5.34E-04 

25 

16 

3.998 

8.86E-04 

20 

17 

4.942 

1.63E-03 

15 

18 

6.325 

3.52E-03 

10 

19 

8.646 

1.01E-02 

5 

20 

13.884 

5.43E-02 

4 

21 

16 

4.07E-02 

3 

22 

19.088 

8.80E-02 

2 

23 

24.28 

2.52E-01 

1 

24 

36 

1.36E+00 



Power 

Supply 


Heater 

Surface 



Figure 3.5 Schematic of switches and resistors for imposed beater 
surface temperature 
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Figure 3.6 


Schematic; of R 1 1 3 degassing unit 
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Figure 3.7 Schematic of the voltage measurement circuit in 
the determination of the gold film temperature 
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Figure 3.9 Schematic of the optical system for orthogonal views in the 
laboratory pool boiling vessel and the drop vessel 
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Figure 3.11 Drop vessel pneumatic diagram 
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Figure 3.12 Laboratory pool boiling vessel 
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Figure j.13 Laboratory pool boiling vessel temperature control 
system 






73 


Test Vessel 


Heaters 



Thermopile Measur .merit 
Junctions 


H*»ter Pow*r 


1 IOy A.C. 


Reference / 
Junctions on Bottom 


Analog Controller 
and Heater Drive 


Thermopile 

Voltage 


Figure 3.14 Drop vessel temperature control system 




Men ler Surface ( 2 ) 



CHAPTER 4 


EXPERIMENTAL PROCEDURE 


The experimental procedures that were followed are described below. Several 
of the test procedures are common to both test vessels used and are presented first: the 
charging of the vessels and the calibration of the thin film heater surfaces. The filling 
of the test vessels will be described first. 

4.1 Common Procedures 

Because of the high solubility of air in R113, precautions were taken in the 
filling process to insure that no contact was permitted between the degassed R1 13 and 
air. The description of the degassing system was given previously in Section 3.1.3. A 
filling cart fabricated from angle iron held a vacuum pump, a valve system and 
associated piping, and a pressure gage. The vacuum pump was used for the evacuation 
of the test vessel and the connecting line to the test vessel. The piping connected the 
test vessel to the tank holding the degassed R1 13, and the valve system permitted the 
evacuation of the test vessel and piping without opening the valve of the R113 
container. The vessel and associated piping and valves were evacuated with a vacuum 
pump for eight hours to remove air and water vapor. The pressure gage provided an 
approximate vacuum indication before the R113 was allowed to flow from the storage 
container. In addition, this filling cart supported the elevated storage tank of the 
degassed R1 13 which provided the hydrostatic pressure for filling the test vessel. 
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The storage container was irradiated with heat lamps to maintain a positive 
pressure relative to the atmosphere. When the pressure gage mounted on the storage 
container indicated a gage pressure of 34.47 kPa (5 psig), the valve between the 
storage container and the stainless steel lines was opened and R1 13 flowed into the test 
vessel. The system was flushed with R113 vapor several times before filling with the 
liquid, to further remove any possible residual gases. To avoid contamination of the 
degassed R113 by atmospheric air, the pressure in the boiling vessels was always 
maintained above the ambient, either by the beilows actuated pressure system or by 
heating. 

For every test, a single point calibration of the thin film heater surface was 
performed before heating the test fluid to the operating temperature of the experiment. 
As discussed in Section 3.1.1, the linear temperature-resistance relationship changed 
with time, with observable changes occurring after about a week. However, the slope 
of the temperature-resistance relationship remained constant, and hence, only a single 
calibration point was necessary for the determination of the new temperature- 
resistance relationship. Prior to the heating of the test fluid, the surface temperature 
was known with certainty since the entire test package was at the accurately known 
ambient temperature, determined with the calibrated thermocouples of the vessel. This 
known surface temperature was determined simultaneously with the heater surface 
resistance for the calibration point data. The procedure for the single point calibration 
of the gold film heating surface was described in detail in Section 3.1.1. 

The experimental procedures required for conducting the experiments with the 
laboratory pool boiling vessel and the drop vessel were slightly different from one 
another due to the variation in vessel instrumentation. Observation of the boiling 
inception times, phase change dynamics, and natural convection phenomena were 
obtained with the high speed camera recordings in the transient experiments. 
Subcooling in all of the tests refers to the difference between the saturation 
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temperature corresponding to the system pressure and the bulk liquid temperature. In 
the following experiments, the total heat flux, q" T , refers to the total thermal power 

generated in the heating surface divided by the area of the heater surface. 

4.2 Laboratory Pool Boiling Vessel 

Transient and quasi-steady experiments were performed in the laboratory pool 
boiling vessel following the calibration of the thin film heating surface. The 
laboratory pool boiling vessel required two hours to reach a steady-state uniform R1 13 
temperature distribution within the vessel, nominally 48.5 °C (with a spatial variation 
of + 0.06 °C as indicated by the internal calibrated thermocouples), as ethylene glycol 
and water circulated in the outer vessel jacket from the constant temperature bath. The 
nominal operating temperature of 48.5 °C was selected since this is slightly above the 
R 1 1 3 saturation temperature (46.7 °C) corresponding to atmospheric pressure. 

Subcooling can be obtained by one of two methods: (1) with a constant liquid 
temperature by pressurizing the liquid above the saturation pressure corresponding to 
the constant temperature or (2) with a constant pressure by cooling the liquid below 
the saturation temperature corresponding to the constant liquid pressure. At the 
beginning of a run, the desired liquid subcooling level of either 0, 2.8, or 11.1 °C, was 
set by the former method with the pressure controller, with the pressure maintained 
constant for each test. 

After the desired subcooling was attained, the camera lights were switched on 
with water circulating through the double-walled infrared water filter, necessary to 
eliminate the heating of the R 1 1 3. A cool-down period of 45 minutes between tests 
provided the same reproducibility of the boiling inception surface temperature (within 
5% of the previous value) as that obtained by waiting 36 hours. This cool-down 
period allowed the heating surface to cool to the original liquid temperature before 
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energization of the heater, and permitted the observable vapor bubbles to collapse. 
Tire procedures for the transient heating measurements will be presented first, 
followed by that for quasi-steady heating. 

4.2.1 Transient Heating Measurements 

The transient experiments performed in the laboratory pool boiling vessel 
imposed either a step in the heater surface heat flux or a step in the heater surface 
temperature. The purpose of using these two different types of transients was to 
examine the effect of the heating rate on the incipient boiling process, as manifested 
by the heater surface temperature at the appearance of the first vapor bubble, and by 
the associated delay time from the onset of heating. 

4.2. 1.1 Imposed Constant Heat Flux 

For the transient imposed constant heat flux tests, a light emitting diode was 
placed in the high speed video camera field of view to synchronize the time when the 
thin film heater was energized with the operation of the camera. A copper knife 
switch connected the heater controller to the heater surface. The 1.5 volt diode was in 
parallel with the knife switch such that when the switch was closed, the diode lighted. 

A command given to the CR7X data acquisition unit via a portable computer 
initiated the data acquisition, and the camera was then switched on. At this point the 
knife switch was closed, and the thin gold film heater and diode were energized. 
Shortly after boiling was initiated and the data taken, the heater circuit was opened at 
the knife switch, and the video recording and data acquisition terminated. 

The independent variables in the imposed heat flux tests were: surface 

orientation, horizontal up and horizontal down with respect to earth gravity; surface 
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type, four different quartz surfaces; bulk liquid subcooling, nominally 0, 2.8, or 11.1 
°C; and total heat flux, nominal levels of 2, 4, 6, and 8 W/cm 2 . The total heat flux is 
the power input to the heater divided by the heater surface area. The dependent test 
variables were the voltage across the heater surface and the voltage across the heater 
shunt, which provided the instantaneous mean heater surface temperature and heat 
flux. The bulk liquid temperature near the heater surface and the system pressure were 
also recorded. 

4.2. 1.2 Imposed Constant Surface Temperature 

All tests conducted with the step increase in surface temperature were 
performed with the thin film heating surface in the horizontal down orientation. With 
heating at a given heat flux in the horizontal down position, the onset of natural 
convection is delayed for a considerably longer period of time than that for heating 
with the identical heat flux in the horizontal up position. The temperature distribution 
in both the liquid and the quartz substrate can thus be determined from a conduction 
analysis with a reasonable certainty, since each experiment is completed before the 
onset of natural convection occurs. This is discussed in detail in Oker (1973). 

As stated in Section 3.1.2, the software for acquiring data in the fast portion of 
the power transient, together with the analog and digital interface card and analog 
input channels, are referred to as the Codas system. With the maximum sampling rate 
limited to 25 kHz per each of two channels, the Codas system was programmed to 
sample the differential voltages across the heater surface and shunt sequentially. The 
CR7X, as previously described in Section 3.1.2, sampled the shunt and heater surface 
voltages only for the slow part of the power transient, but with much greater accuracy 
than was possible with the Codas system. The measurement accuracy was thus 
compromised with the use of the Codas system for the high sampling rate, while the 
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high sampling rate was compromised by using the CR7X for the measurement 
accuracy. Details of the measurement accuracy are given in Appendices C and Dl. 
The pressure transducer voltage and bulk fluid temperature thermocouple voltages 
were also recorded by the CR7X at the start and end of a given test. The onset of 
boiling and initial vapor dynamics were recorded by the high speed video camera. The 
imposed steps in heater surface temperature were in the range of 0.5-25 °C. 

An experiment was initiated by turning the camera on, followed by execution 
of the switching program on the first computer, which caused a light emitting diode on 
one of the circuit boards to flash, indicating that the heater surface would be powered 
in one second. At the time of the one second warning, the Codas data acquisition 
process began on the second computer. The test was terminated after heating for 1.8 
seconds by the switching software, which opened the heating circuit automatically. 

The independent variables in the imposed heater surface temperature tests 
were: surface type, two different quartz surfaces; bulk liquid subcooling, nominally 

0, 2.8, or 11.1 °C; and the initial heater surface voltage level, which governed the level 
of the step increase in heater surface temperature. The dependent test variables were 
the voltage across the heater surface and the voltage across the heater shunt, which 
provided the instantaneous mean heater surface temperature and heat flux. The bulk 
liquid temperature and system pressure were also recorded. 

4.2.2 Quasi-Steady Heating 

Both the thin gold film heater on quartz and the solid metal heaters were used 
for the quasi-steady measurements. The desired power setting was set, and the system 
was allowed to reach a steady-state surface temperature. If boiling inception did not 
occur, the heater input voltage was increased by 10% of the previous level, and 
another steady-state heating surface temperature was obtained. Following the point of 
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boiling initiation, the power was switched off and the surface allowed to cool for more 
than 45 minutes. Power at a slightly lower level than that which produced the 
previous boiling inception point was then applied to the surface to obtain a non-boiling 
steady-state surface temperature. The power was then increased slightly again in 
steps, until boiling occurred. This incremental approach was iterated to determine 
more precisely the temperature and heat flux at which incipient boiling occurs. 

Because of the large uncertainties as to when incipient boiling would take place 
with quasi-steady heating, high speed photographs were not taken. The CR7X 
recorded the shunt and heater surface voltages, for the calculation of the heat flux (and 
mean heater surface temperature in the case of the gold film heater on quartz), the 
pressure transducer voltage, for the determination of the pressure level, and the voltage 
output of the thermocouple imbedded in the copper block near the heater surface for 
the determination of the metal heater surface temperature. The bulk fluid temperature 
increased by no more than 0.25 °C during any of the quasi-steady heating tests. 

The independent test variables in the quasi-steady heating tests were: surface 
orientation, horizontal up or horizontal down with respect to earth gravity; surface 
type, quartz substrate or metallic surface; and bulk liquid subcooling, nominally 0, 2.8, 
or 11.1 °C. The dependent test variables were the voltage across the heater surface, 
the voltage across the heater shunt, and the metal heater thermocouple voltage, which 
provided the mean heat flux and heater surface temperature. The bulk liquid 
temperature and system pressure were also recorded. 

4.3 Drop Vessel 

Following charging with degassed R 1 1 3, the drop vessel, cameras, lights, and 
associated electronics were placed on the drop bus, which was then balanced with lead 
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weights. Only one of the two thin film heaters installed in the drop vessel was used 
for a particular test. 

The day before a drop test, the single point calibration of the heater surface 
resistance-temperature was performed, since the entire drop vessel was at the 
accurately known and uniform ambient temperature. The measured shunt, surface, 
and thermocouple voltages acquired by the CR7X during the surface calibration were 
stored in a portable computer. The a.c. temperature controller regulating the power to 
the thin film heaters mounted on the outer drop vessel walls was switched on, and the 
vessel was then covered with an insulating blanket. The controller was adjusted to 
provide a temperature of approximately 58 °C when the vessel was to be prepared for 
installation in the drop tower on the following day. 

Preparation for the drop test the next day then required that the a.c. heaters be 
disconnected from the power source. The heat loss to the ambient at about 20 °C 
during the 20 minute period required to transfer the vessel to the actual drop location 
and during the two hour period required for the evacuation of air from the drop tower 
resulted in a decrease in temperature of the test vessel and its contents from the initial 
58 °C to the desired level of about 47.7 °C. 

Three tests were performed sequentially in connection with each drop 
experiment: a pre-drop test with the high speed photography, a drop test with high 
speed photography, and a post drop test with no photography. The pre-drop test was 
identical to the drop test, using the drop tower control room circuitry, but was 
conducted at earth gravity. This test provided a direct comparison between the boiling 
process in earth gravity and boiling in microgravity, and included the transient heater 
surface temperature, heat flux, pressure, and photographic measurements. The 
exposed film was removed for developing after completion of the pre-drop test, and 
the data downloaded from the CR7X to a portable computer. The camera was then 
reloaded for the drop test. After the test package was dropped and retrieved from the 
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bottom of the drop tower, the film was again removed for processing, and the drop 
data downloaded from the CR7X to a portable computer. The post drop test was then 
performed to verify the proper operating condition of the circuitry and mechanical 
hardware following a drop test, and was identical to the pre-drop test except that no 
movies were filmed. 

The actual drop test procedure will now be described in greater detail. 
Following the pre-drop test, the camera was reloaded, the pre-drop data was 
downloaded to the portable computer, and the batteries were recharged. Data 
acquisition and control instructions were transferred from the portable computer to the 
CR7X mounted on the drop bus. The software instructed the CR7X to measure the 
bulk liquid temperature and then set the specified subcooling level, based on this 
measured bulk liquid temperature, by adjusting the pressure level inside the vessel 
through actuation of the bellows pressure control system described previously. The 
CR7X also set the voltage level to be applied to the heater surface to produce the 
desired heat flux. The CR7X control instructions also checked that appropriate 
voltage signals were being received from the drop tower control room by means of the 
drop bus umbilical. The protective outer skins were installed on the drop bus, and the 
bus placed into position at the top of the drop tower, at which time evacuation began. 

At the desired tower vacuum level and R113 temperature, the control room 
sent a signal to the drop vessel (called the energize signal) which caused the CR7X to 
set the pressure level based on the current vessel temperature in order to achieve the 
desired subcooling level. The CR7X sampled the pressure transducer signal and 
calculated the pressure based on the pressure transducer calibration. The pressure set 
point (command voltage) was adjusted until the measured pressure was near the 
desired pressure. At this point, the CR7X transmitted a five volt signal to the control 
room, indicating that the R1 13 in the vessel was at the desired subcooling. The control 
room then switched on the camera lights and started the high speed camera, which was 
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followed by the release of the drop bus in the drop tower. The onset of the free fall 
then disconnected the umbilical from the test package, which in turn produced an open 
circuit signal called the drop signal. The drop signal was an indication to the CR7X to 
(1) power the heater at the pre-determined voltage level by sending a calibrated set 
point signal between zero and five volts to the heater controller and to (2) begin data 
acquisition. A diode powered by a three volt digital signal from the CR7X in the field 
of view of the camera was used to synchronize the application of power to the heater 
surface with the camera operation. After five seconds, the CR7X stopped measuring 
the shunt, heater surface and pressure transducer voltages, and terminated power to the 
heater surface. In some of the drops, the voltage of a thermocouple in the liquid near 
and over the center of the heater surface was sampled. The camera continued to film 
during impact. An outline of the free fall drop sequence is presented in Figure 4.1. 

The independent variables in the drop vessel tests were: surface type, two 
different quartz surfaces; bulk liquid subcooling, nominally 0, 2.8, or 11.1 °C; and 
total heat flux, nominal levels of 4, 6, and 8 W/cm 2 . The low heat flux level of 2 
W/cm 2 used in the laboratory pool boiling experiments was not used in the drop 
experiments because it was observed that the time required for the onset of boiling at a 
power input of 4 W/cm 2 in microgravity was apparently greater than the five second 
microgravity period available in the drop tower. The dependent test variables were the 
voltage across the heater surface, the voltage across the heater shunt, and the pressure 
transducer voltage, which provided the mean instantaneous heater surface temperature, 
heat flux, and system pressure. As will be noted in the results, short pressure 
variations occurred during some of the more dynamic nucleation processes, although 
attempts were made to maintain a constant pressure level. At times, a liquid 
temperature measurement near the heater surface was also included. 
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4.4 Data Reduction 

The procedures for the calculation of the thin film heater surface temperature, 
the metal heater surface temperature, the bulk liquid temperature, the liquid pressure, 
and the heat flux from the recorded data are presented below, along with the 
uncertainties in the determination of these quantities. 

4.4. 1 Spatial Mean Film Heater Surface Temperature 

The spatial mean thin gold film heater surface temperature, T w , is determined 
from the measured test surface resistance, R w , the single point calibration resistance, 
R^, with its corresponding calibration temperature, Tq and the slope of the calibration 
curve, in the following equation: 


T 

W 


m dT 
T c + dR 


R -R r 

w C 


(4.1) 


For the constant heat flux tests, R w , R^, T c , and ^ were determined with 


representative values and uncertainties of 2.9985 ± 0.0008 Cl, 2.6266 + 0.0002 £2 , 
20.00 ± 0.01 °C, and 214.83 + 0.23 °C/Q, respectively, as given in Appendix D.l. As 
a result of the uncertainties in these quantities, the heater surface temperatures were 
determined with an uncertainty of + 1.0 °C, as described in Appendix D.l. 

For the imposed step change in surface temperature tests using the Codas 
system during the fast portion of the power generation curve discussed earlier, Ry, R^, 
T(-, and ^ were determined with representative values and uncertainties of 2.7989 ± 


0.020 Cl, 2.6266± 0.0002 Cl , 20.00 ± 0.06 °C, and 214.83 ± 0.23 °C/Q, respectively. 
The calculation of these uncertainties is also described in Appendix D.l. Asa result of 
these uncertainties, the heater surface temperatures were determined with an 
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uncertainty of + 6 °C. When the CR7X was used during the slow portion of the power 
generation curve, the heater surface temperature was determined with an uncertainty of 
± 1 °C, as shown in Appendix D.l. 

4.4.2 Metal Heater Spatial Mean Surface Temperature 

The metal heater surface temperature was determined from the chromel- 
constantan thermocouple imbedded in the copper heating block near the heater surface, 
as shown in Figure 3.2. The error in assuming that the temperature at the location of 
the thermocouple was equal to the surface temperature was estimated to be less than 
0.002 °C, by assuming one dimensional conduction in the copper block and negligible 
contact resistance between the soldered copper foil and the copper block. An error of 
this magnitude is less than the uncertainty of ± 0.01 °C associated with the 
thermocouple temperature measurement and, as a result, was neglected. 

4.4.3 Bulk Liquid Temperature 

The liquid bulk temperature was measured by calibrated chromel-constantan 
thermocouples, described in Chapter 3. The CR7X was used for all measurements of 
bulk liquid temperature in all three vessels, and the uncertainty of + 0.6 microvolts for 
the thermocouple voltage determination translates to a + 0.01 °C uncertainty in the 
liquid bulk temperature. A detailed discussion of the estimate of this uncertainty in the 
temperature determination from a thermocouple voltage measurement is in Appendix 
D.2. 
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4.4.4 System Pressure 


The liquid pressure near the heater surfaces in all three vessels was measured 
by calibrated transducers. The calibration and resulting calibration equations of the 
pressure transducers are described in Appendix B.3 The uncertainty of the pressure 
measurement in the laboratory pool boiling vessel, the drop vessel, and the forced 
convecdon boiling loop was less than the desired uncertainty of + 0.172 kPa (± 0.025 
psi). The details of the estimate of the uncertainty in the pressure measurements are 
given in Appendix D.3. 

4.4.5 Heat Flux 


The total power input per unit area, q" T , for both the quasi-steady heating tests 
and the constant heat flux tests was calculated from: 
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T sh 


Q T 


A R , 

sh 


(4.2) 


The voltages were in volts, the resistance in ohms, the heater surface area in cm 2 , and 
q” T in W/cm 2 . The maximum uncertainty in q" T was ±0.18 W/cm 2 . The details of the 
estimate in the uncertainty associated with q" T are presented in Appendix D.4. 

The heat transfer to the liquid, q"|, could be calculated from the measured total 
heat flux, q" T , under circumstances which permitted reasonable analytical assumptions. 
For the transient constant heat flux tests in the pool boiling vessels q", can be 
determined by a one dimensional, semi-infinite media conduction analysis provided 
two conditions were met. The first condition for a valid semi-infinite media 
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conduction analysis is the absence of significant natural convection while the fluid is 
being heated. This condition was met for all of the microgravity boiling tests up to the 
point of the onset of boiling. The conduction analysis was valid for experiments 
conducted in earth gravity where boiling occurred in times less than that required for 
the onset of natural convection. The time for the onset of natural convection for a 
heated surface horizontal up in earth gravity was given by Oker (1973) in a semi- 

empirical equation for the combination of a pyrex substrate and R-113 at atmospheric 
pressure, where q" T is in W/cm 2 and t^ is in seconds: 



0.77 


(4.3) 


For illustrative purposes, a q" T of 8 W/cm 2 corresponds to a time, t^, of 0.27 seconds 
required for the onset of natural convection. Equation 4.3 is based on the fluid near 
the heating surface reaching a certain critical Rayleigh number. This condition was 
not met by most of the pool boiling experiments conducted here in earth gravity and, 
as a result, q"! was not computed for most of these earth gravity tests. The second 

condition necessary for a valid semi-infinite media analysis is that the thermal 
penetration depth be less than the thickness of the heated substrate material. An 
acceptable thermal penetration depth in this study corresponds to periods of heating 
that are less than six seconds. All of the microgravity tests were completed in less 
than six seconds, as were many of the laboratory pool boiling experiments with the 
heating surface horizontal down. 

Using the one dimensional semi-infinite media solutions presented in Appendix 
G, for a constant plane heat source, the ratio of the heat flux transferred to the liquid 
and to the quartz is a constant, given by: 
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Evaluating the properties of quartz and R 1 1 3 at a mean temperature of 70 c C, the 
fraction of the total input heat flux to the liquid becomes, from equations 4.4 and 4.5: 


— = 0.175 


(4.6) 


The heat flux to the test fluid is then evaluated from equation 4.6 using the 
measurement of q" T . The uncertainty in the calculation of heat transfer to the liquid as 

given by equation 4.6 comes from the uncertainties in the tabulated property values 
and the uncertainty involved in using a representative mean temperature for the 
selection of the tabulated properties. The total uncertainty in q", calculated from 
equation 4.4 is ± 4%. In the non-boiling portion of the quasi-steady tests q”| is 
computed from subtraction of the losses from the total power input to the heater. The 
losses are estimated from the following natural convection correlations (Incropera and 
Dewitt (1984)) for the test surfaces heated in air: 


\ 

Nu = 0.15 Ra 3 (horizontal up) 


(4.7) 
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l 

Nu = 0.27 Ra“* (horizontal down) (4.8) 

The uncertainty in q", results primarily from the uncertainty in the calculation of the 

losses computed from the natural convection heat transfer correlations, which are on 
the order of + 20%. 
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Drop Sequence 


Initial State: 

Program downloaded to CR7 datalogger 
Command voltage to heater = 0 V 
Power (+24 )V to timer card 

t=- 5 Minutes: E nergise «igna l (open contact) 


1 

Control 

\ 


U Room 

[ 

Datalogger 

r 


t=0(Drop): 



—? l . s*t pressure based on temperature 
to get subcooiing by fitting command 
voltage 

2. Set relay *3 to turnon+24V lor 
pressure controller /instrumentation 

| ^ 3. Open shutoff valye# 

— ^ 4. tfeasure pressure, adjust command 
voltage until measured pres sure is 
close to set point 


Drop signal 
generated by 

umbilical release 


-) 5. Switch relays 1 and 2 to switch Irom 
calibration to test circuit 





Datalogger 




-) 1 . Turn on surface heater and diode 
L — } 2. Acquire voltage data 


t=5.0 sec (Impact): 


Initiate activities 
alter timer stops 


Datalogger 


1. Stop taking data 

— )2. Turnoff heater and diode 

^3. Switch relays 1 and 2 

— ^ 4. Close shutoff valves 
— )5. Switch relay *3 turning 
off instrumentation 
power 


Figure 4.1 Sequence of events during a free fall microgravity test in 
the NASA Lewis RC drop tower 
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CHAPTER 5 

EXPERIMENTAL RESULTS 

The experimental results obtained during the current study are presented 
below. The transient experimental results are presented first, followed by the quasi- 
steady heating results. 

5.1 Transient Boiling 

The transient boiling experiments were performed in both the laboratory pool 
boiling vessel and the drop vessel. A step in heater surface heat flux was used in the 
drop vessel. Both step increases in heat flux and in surface temperature were used in 
the laboratory pool boiling vessel, as well as the quasi-steady heating process which 
will be presented in Section 5.2 below. A step change in heat flux is the most 
elementary form of imposed heat flux variation possible, as all other time varying 
imposed surface heat fluxes can be constructed from a series of steps in heat flux. 
These two basic and elementary thermal boundary disturbances, step changes in heat 
flux and in temperature, produce distinctly different liquid temperature distributions, 
which determine the heater surface temperature at the moment of boiling inception and 
the resulting phase change dynamics. The transient pool boiling results are presented 
first, followed by the transient forced convection boiling results. 
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5.1.1 Pool Boiling-Imposed Constant Heat Flux 

5. 1 . 1 . 1 Transient Heater Surface Temperature 

Since the resistance of the gold film heater surface is computed from the 
measured current and voltage drop across the entire heater surface, the resulting 
calculated surface temperature is a spatially averaged temperature. The non- 
uniformity of the surface temperature was estimated by solving the three dimensional 
transient conduction model for the quartz substrate-Rl 13 system using a finite 
difference solution technique. The conduction assumption with R1 13 is valid for short 
times in earth gravity and in microgravity up to the point of incipient boiling. The 
Crank-Nicholson finite difference method was used, with a uniform spatial grid size of 
0.79375 mm (0.03125 inches). The procedure is described in Appendix H. Figure 5.1 
shows the computed surface temperature rise for a total heat flux input of 8 W/cm 2 at 
one second, the approximate time from the onset of heating for boiling inception. 
Because of symmetry, only one quarter of the film surface is shown, and the computed 
surface temperature is noted to be quite uniform except near the heater surface edges. 

In the remainder of this section, selected representative heater surface 
temperature-time plots are presented in order to demonstrate the distinguishing 
features in common for the same body force magnitude and orientation. The identical 
heater surface, designated as Q-5, is used in all of the results presented here. 

The transient heater surface temperature for Run 79 is shown in Figure 5.2, 
along with a tabulation of the test conditions, and applies to the case where the heater 
surface is facing horizontal up in earth gravity. The onset of natural convection, 
appears as an irregularity in the temperature-time plot. The time from the energization 
of the heater to the onset of natural convection is designated as t^ . Experimental 
values for this t nc were obtained by the appearance of two independent phenomena. 
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which coincide with the irregularity: the observed onset of fluid motion above the 
heater was characterized by a wave-like disturbance recorded by the high speed video 
camera, and by the departure of the heater surface temperature from the one 
dimensional semi-infinite media transient conduction solution. 

Oker (1973) found a semi-empirical relationship for based on the critical 

Rayleigh number predicted by instability theory for the case of a semi-infinite surface 
in the horizontal up position in earth gravity with an imposed heat flux: 
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a, Pgq" 
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(5.1) 


Ra c is the critical value of the Rayleigh number for the onset of fluid motion, where 
the Rayleigh number is defined as: 
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5 is the approximation to the penetration depth of the thermal disturbance, and was 
taken as v/6at in the work of Oker, from the integral solution of the one dimensional 
transient conduction equation with a step in heat flux for a media of semi-infinite 
extent. With the properties evaluated at a mean temperature of 66 °C and Ra c taken as 

800, as suggested by Sparrow et al (1964) for a semi-infinite surface with an imposed 
step in heat flux without an upper rigid boundary, equation 5.1 becomes: 
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(5.3) 


In Figure 5.3, the prediction of Oker for t nc , equation 5.3, is less than the 
experimentally measured t nc of this study. When 8 is taken as 2v/at from the exact 
one dimensional solution to the conduction equation for an imposed surface heat flux 
with the properties also evaluated at a mean temperature of 66 °C, the following 
equation is obtained: 
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0.463 
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(5.4) 


Equation 5.4 predicts times for the onset of natural convection which are also lower 
than the corresponding experimentally measured times. Using a larger Ra c of 3600, 

with 5 taken as 2y[at , gives: 


t 

nc 


0.982 



(5.5) 


and describes the experimental measurements better than the previous expressions. 
Since the heating surface is not semi-infinite (w/1 = 0.5), the edge effects of the film 
heater during heating can be expected to affect the onset of natural convection, 
requiring a Ra c different than that predicted by the earlier analyses based on a semi- 
infinite heating surface. 

In an earlier study, Ngheim (1980) observed the onset of natural convection 
from a narrow gold film heater (w/1 = 0.1) on a large quartz substrate in R113 with 
high speed holographic interferometry simultaneously with the measurement of the 
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film heater resistance. The resulting temperature-time plots, similar to that in Figure 
5.2, also showed the same type of disturbance at t„ c , and were also correlated well by 
equation 5.5. 

The next significant event in Figure 5.2 following the onset of natural 
convection is incipient boiling, as indicated. Incipient boiling, or the onset of boiling, 
is defined as the appearance of the first vapor bubble on the heating surface. In some 
tests, this incipient boiling takes place at the point when the mean heater surface 
temperature reaches a maximum, while in other tests, the onset of boiling occurs prior 
to this point. For tests in which the latter took place, the level of the ensuing 
maximum heater surface temperature depended on the manner in which the boiling 
propagated across the heated surface. If incipient boiling occurred as an almost 
explosive event over the entire heater surface, then the temperature associated with 
boiling incipience was the maximum surface temperature. If incipient boiling 
occurred at a location near a corner of the heater, for example, the heater would 
become cooled locally while the remainder of the heater surface continued to rise in 
temperature. The heater surface temperature measurement is a spatially averaged 
quantity, as described earlier, and the measurement would continue to rise until the 
boiling spread sufficiently to produce a subsequent decrease in the mean value. This 
boiling propagation will be defined in more precise terms later. 

Following the onset of boiling, a quasi-steady boiling region is noted in Figure 
5.2. In this domain of the temperature-time plot, boiling has spread across the entire 
heating surface, and the quasi-steady boiling temperature level is less than the 
maximum heater surface temperature, but above the saturation temperature for the 
liquid at the system pressure, as expected. 

The heater surface temperature-time plot of test Run 75 is presented in Figure 
5.4 with the experimental conditions as indicated. These experimental conditions are 
virtually identical to those of Figure 5.2, except that the heater surface is inverted, or 
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horizontal down (a/g = -1). As can be seen, the temperature-time curve is quite 
different from that for heating in the horizontal up position. The departure from the 
one dimensional conduction prediction marks the onset of natural convection, t^, as 

was also observed in Figure 5.2 for the horizontal up surface orientation. In this case 
the fluid flow that eventually takes place with heating results in a lateral outward flow 
for the heated surface in the horizontal down position, and depends on the dimensions 
of the heater surface and on the geometry of the enclosure. The buoyancy force tends 
to hold the liquid against the heated surface. A uniformly heated semi-infinite surface 
with this orientation produces a stable configuration with no liquid motion, and heat 
transfer by conduction alone exists. Heaters of finite size, however, produce lateral 
temperature gradients in the liquid, which induce subsequent lateral liquid motion. In 
the present study, the edge effects on the liquid temperature resulting from the finite 
extent of the heater surface were such that the t Tlc for a heater horizontal down, with 

the same imposed heat flux and initial conditions as a horizontal up heater, was larger 
than t nc for a horizontal up heater (1.9 seconds versus 1.2 seconds in Figures 5.4 and 

5.2). Similar observations were reported in the transient experimental work ol Oker 
(1973). 

The incipient boiling point is indicated in Figure 5.4, and appears as a 
disturbance in the temperature-time plot. The initial phase change and propagation of 
this phase change across the heated surface causes a drop in the spatial average 
temperature. Once boiling has spread across the entire heating surface, the vapor 
bubble covers the heater, since buoyancy holds the R113 vapor at the heater surface 
with a/g = -1, unlike heating with a/g = +1 where the buoyancy force lifts the vapor 
from the heated surface. Since the thermal conductivity of R113 vapor is much less 
than that of the liquid, the heater becomes, in effect, insulated from the liquid. As a 
consequence, the temperature of the heater will rise, if permitted, until the heater 
reaches an unduly high temperature or is destroyed. This process is sometimes 
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referred to as dryout, or heater burnout. Boiling with a surface horizontal down (a/g = 
-1) in the absence of forced liquid motion will result only in an unsteady heater surface 
temperature. Periodicities superimposed on the mean surface temperature are 
observed in Figure 5.4 because the vapor slug continuously collapses and grows near 
the edges of the heater. 

Transient heating in microgravity is different from that in earth gravity, 
regardless of whether heating horizontal up or horizontal down. Figure 5.5 presents 
the spatial mean temperature-time plot for microgravity drop vessel experiment 
PBMT1 1 1 5.605, which used a heated thin gold film surface. The measured 
temperature follows the one dimensional semi-infinite media conduction solution 
reasonably well until incipient boiling takes place. At the level of microgravity 
present in these initially quiescent experiments (a/g = 10' 5 ), there was effectively no 
liquid modon resulting from the heating of the liquid. Following the onset of boiling, 
the surface temperature decreases and then immediately rises, as was expected after 
the drop tower movie films were viewed. Following the initial boiling, a single large 
vapor bubble covered the heater surface, effectively insulating the heater surface. 
Similar to boiling in earth gravity in the horizontal down orientation, boiling on a flat 
surface in microgravity is unsteady. 

Figure 5.6 presents the predicted and measured thin gold film surface 
temperature rise for a variety of tests with surface Q-5, a nominal q" T of 6 W/cm 2 and 

nominal subcoolings of nominally either 0 or 11 °C. The final two digits of the test 
names listed in Figure 5.6 give the subcooling in degrees Fahrenheit. The results are 
scaled by the total heat flux to both the R 1 1 3 and the quartz substrate, and plotted as a 
function of the square root of time. The predicted heater surface temperature rise is 
from the solution to the one dimensional transient heat conduction problem for two 
semi-infinite media with constant energy generation at the plane separating the two 
media given in Appendix G and evaluated at the common plane: 
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(5.6) 


Measurements from tests with heater surface orientations of both horizontal up and 
horizontal down in earth gravity are included along with data from two microgravity 
tests. Departure from the linear prediction results either from the onset of liquid 
motion arising from the buoyancy effects during heating or from the onset of boiling. 
Although heating of the film surface is not truly one dimensional, the departure from 
the one dimensional linear prediction shown in Figure 5.6 provides a conservative 
estimate for the onset of natural convection with a surface in the horizontal up or 
horizontal down orientation in earth gravity, as was demonstrated in Figure 5.2. For 
the microgravity experiments PBMT1 101.600 and PBMT1 121.620, the departure from 
the conduction predicted surface temperature was the consequence of incipient boiling. 
The lack of departure from the linear prediction earlier than at incipient boiling gave 
credence to the technique used for heater surface temperature measurement. 

Figure 5.7 presents the results of constant imposed heat flux tests in the 
laboratory pool boiling vessel, with a surface orientation of horizontal up. The tests 
were all conducted with the same heat flux but with three different levels of 
subcooling. For the early times the curves are superimposed on one another, as 
anticipated. However, at later times they were expected to separate because of the 
differences in subcooling and delay times to incipient boiling. Figure 5.7 confirms 
this expectation, and provides an indirect check on the heater surface temperature 
measurement. Because of the difference in subcooling, the quasi-steady surface 
boiling temperature is different for each of the surfaces by roughly the amount of the 
liquid subcooling, as indicated. 
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5. 1 . 1 .2 Transient Liquid Temperature 

"The measured heater surface temperatures provide data with which the 
calculated temperatures can be compared. In the absence of liquid motion, it is also 
possible to compute the transient temperature distribution in the liquid. A calibrated 
thermocouple was included in the drop vessel, and transient measurements obtained 
which can be compared with the computations. Figure 5.8 shows the temperature rise 
determined by the thermocouple located at a nominal distance of 1.6 mm above the 
center of the heater surface. The uncertainty in the location of the thermocouple 
relative to the heating surface was + 0.2 mm, and the calculated temperature variations 
are presented for the possible locations of the thermocouple at 1.4 mm, 1.6 mm, and 
1.8 mm above the heater surface. The calculations for the liquid temperature at these 
locations are from the solution to the one dimensional transient conduction equation 
for two semi-infinite media given in Appendix G. The measured thermocouple 
temperature variation is closest to the liquid temperature variation at 1.6 mm, the 
nominal location of the thermocouple. 

5. 1.1.3 Transient Pressure Measurements 

The laboratory pool boiling vessel and the drop vessel pressure transducers 
measured the transient pressure of the liquid near the heating surface during the 
experiments. This section presents representative results of the transient pressure 
measurements, which give an indication of the quality of the pressure control possible 
with set pressure levels in the laboratory and drop vessel pool boiling experiments. 

Figure 5.9 shows representative pressure measurements for boiling tests 
performed in the laboratory pool boiling vessel at the three pressure levels used to 
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obtain three different subcooling levels with an initial liquid temperature of 48 °C. 
The pressure level was maintained within ± 0.34 kPa (± 0.050 psi) of the set point 
except at boiling incipience, where the pressure increased by as much as 15 kPa ( 2.2 
psi) as a result of the sudden generation of a considerable amount of vapor. This spike 
in the pressure measurement, when it was large enough to be observed, always 
correlated with the incipient boiling observed on the recorded videos. When the 
boiling incipience was particularly violent, as seen on the recorded videos, a pressure 
system oscillation was observed. The response time of the pressure control system 
required to restore the pressure to the original level before boiling was less than 30 
milliseconds, except for the more violent incipient boiling events. 

Figure 5.10 shows representative pressure measurements obtained with the 
drop vessel experiments in microgravity. The three pressure levels necessary to 
produce the three different subcooling levels used in microgravity are presented, and 
one pressure level from a pre-drop test measured with the heater surface horizontal up 
in earth gravity is included for comparison. In the drop vessel pressure measurements, 
as in the laboratory pool boiling measurements, a spike of as much as 6.8 kPa (1 psi) 
sometimes occurred at incipient boiling. The pressure level in the drop vessel 
otherwise was maintained to within + 0.6 kPa (± 0.09 psi) of the set point. On 
occasion, however, an initial small decrease in pressure occurred in changing from 
earth gravity to microgravity on release of the test package before the pressure reached 
a nearly constant level, because of the dynamics of the pressure control system. 
Incipient boiling occurred after this initial pressure disturbance, so this anomaly in 
pressure was not detrimental to the results obtained. 



102 


5.1. 1.4 Transient Heat Flux Measurements 

Figure 5.11 illustrates representative examples of transient total heat flux data 
for boiling tests performed in the laboratory pool boiling vessel and in the drop vessel. 
The heater controller maintained the heat flux constant to within ± 3% of the set point 
for most of the laboratory pool boiling vessel experiments, and to within ± 5 % of the 
set point for most of the drop vessel experiments. 

5. 1.1. 5 Heater Surface Superheat at Boiling Inception 

Four different gold thin film surfaces were used in the tests conducted with an 
imposed constant heat flux and were designated by: Q-l, Q-2, Q-5, and Q-6. Surfaces 
Q-5 and Q-6 were used in both the laboratory pool boiling vessel and the drop vessel, 
while surfaces Q- 1 and Q-2 were used only in the laboratory pool boiling vessel. The 
heat flux imposed on these surfaces, q T " , the delay time to incipient boiling, t*, and 
the heater surface superheat temperature at the time of boiling inception, AT* sup , are 
related to one another. The following figures will attempt to demonstrate these 
relationships in a variety of ways. It will become obvious, of course, that as t* is 
reduced for whatever reason, in the absence of natural convection effects AT* sup will 
also be reduced for a given level of q" T . Certain data points shown are flagged, and 
represent the heater surface superheat at the time of boiling inception. The flagged 
data point is connected by a dashed line to an unflagged data point, which in turn 
represents the heater surface superheat measured at the time when the boiling begins to 
propagate across the surface. Data points with no flags or connecting lines represent 
data where the onset of boiling resulted in virtually the immediate spread of the 
boiling process. This amounted to a delay of less than 1 millisecond for tests 
conducted in the laboratory pool boiling vessel, and of less than 3 milliseconds in the 
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drop tower. This latter figure is greater because of the uncertainty associated with the 
lower camera speed. Nominal subcoolings of 0, 2.8, and 11.1 °C were used in all of 
the tests. With the heater surface horizontal up, a q T " of 2 W/cm 2 resulted in a boiling 

inception delay time of more than five minutes. As a result, this combination of 
orientation and heat flux was not used in any of the laboratory vessel boiling tests. 
The lower bound of the heat flux range for the microgravity tests was set by the 
limited fall time of 5.18 seconds in the drop tower. In two drop vessel microgravity 
experiments with a nominal q" T of 4 W/cm 2 , boiling inception did not occur. As a 

result of these two preliminary tests, the heat flux range for the microgravity tests was 
limited to 6 to 8 W/cm 2 . 

The heater surface superheat results for surface Q-l shown in Figure 5.12 were 
obtained in the laboratory pool boiling vessel with the heater surface in either the 
horizontal up or the horizontal down position. For the initial test conditions as 
indicated, Figure 5.12 implies a significant time interval between boiling inception and 
boiling spread at the lower heat flux of q" T = 4 W/cm 2 which was not observed for any 

of the other surfaces. Figure 5.12 also indicates a trend of a decrease in the heater 
surface superheat associated with the spreading of the boiling process as the heat flux 
is increased. This decrease was not observed with the other heating surfaces, where 
the other corresponding heater surface superheats tended to remain constant as the heat 
flux increased. The onset of boiling with Q-l always occurred at a large visible 
scratch on the thin gold film. From the presence of this scratch and the inconsistency 
between these data and that from the other surfaces, the results in Figure 5.12 are 
considered to be an anomaly. Since the thin gold film was undamaged prior to 
installation, it is surmised that the heater surface was accidentally scratched as it was 
placed in the vessel. 

Figure 5.13 show's the results of wall superheat at boiling inception from the 
laboratory pool boiling experiments for surface Q-2. The initial conditions are given. 
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Larger heater surface superheats at boiling incipience resulted for a/g = - 1 than for a/g 
= + 1 with q T " near 4 W/cm 2 , but not near 8 W/cm 2 , for tests with the same nominal 

subcooling. The wall superheat for most tests fell into a band of data between 40 and 
55 °C, except for wall superheats corresponding to q T " near 8 W/cm 2 which produced 

lower heater surface superheats. 

Figure 5.14 presents the heater surface superheat incipient boiling results for 
surface Q-5 for a/g = ± 1 and microgravity, with different subcoolings. A nominal 
initial fluid temperature of 49 °C was used for the a/g = ± 1 tests and a nominal initial 
fluid temperature of 47 °C was used for the microgravity tests. For total heat flux 
levels of 2, 4, 6, and 7.5 W/cm 2 , the maximum superheats at boiling incipience 
occurred at a/g = -1 or microgravity, and at a nominal subcooling of either 0 or 2.8 °C. 
Boiling inception followed by a measurable time delay to boiling spread occurred for 
q x " greater than 7 W/cm 2 . It appears from the surface Q-2 and Q-5 data in Figures 

5.13 and 5.14, respectively, that the time delay to boiling spread is most pronounced at 
the higher heat flux levels. Except for the initial boiling phenomena near the heat flux 
of 8 W/cm 2 , all of the heater surface superheat at boiling inception fall within a range 
of 50 to 75 °C. 

For all tests conducted with surface Q-5 after 10/13/89, the superheats for 
boiling incipience consistendy fell in the range of 50 to 75 °C. Before this date, the 
superheat level increased continually from about 10 to 38 °C. The boiling incipience 
character of the surface was somehow changing during this interval, after which the 
heater surface superheat at the onset of boiling became consistent. Similar behavior 
was observed by Tong et al (1990) in an experimental study of boiling of R113 with 
thin platinum thin films. Tong et al speculated that the R113 had somehow broken 
down above 100 °C and formed deposits on the heater surface. However, the 
decomposition temperature of R 1 1 3 is greater than the observed heater surface 
temperatures of the present tests: a trace decomposition of R1 13 exposed to quartz for 
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30 seconds was found to first occur at a temperature of 299 °C in a testing described 
by E.I. DuPont de Nemours (1969). It is believed that the surface became more 
"regular" in incipient boiling character after the first few tests and the resulting boiling 
removed some of the adsorbed gases and impurities. 

Figure 5.15 shows the heater surface superheat at incipient boiling for surface 
Q-6 as a function of total heat flux with a/g = + 1 and microgravity with different 
subcoolings. The tests at a/g = -1 were conducted in the laboratory boiling vessel, 
while those in microgravity were conducted following the tests for surface Q-5. 
Surfaces Q-5 and Q-6 were installed in the drop vessel at the same time. The nominal 
initial fluid temperature for the earth gravity tests was 49 °C and 47 °C for the 
microgravity experiments. Except near 8 W/cm 2 , the heater surface superheats lie in a 
range between 48 and 62 °C. The lower level of this heater surface superheat range is 
close to the lower range of heater surface superheats for surface Q-5. The heater 
surface superheat at incipient boiling is lower than that for the undamaged surfaces Q- 
2 and Q-5, at the high heat flux levels. 

The incipient boiling heater surface superheat levels for surfaces Q-5 and Q-6 
were similar, but greater than those of Q-l and Q-2. The similarity of behavior for 
surfaces Q-5 and Q-6, compared to surfaces Q-l and Q-2, is presumed to be due to 
their having been fabricated together by sputtering, while surfaces Q-l and Q-2 were 
vacuum deposited with gold at the same time. Surface Q-2 shows a more constant 
superheat range than does Q-l, most likely because of the presence of the scratch 
observed on Q-l. In all of the plots of AT* sup , the delay from incipient boiling to the 

onset of boiling spread is most pronounced near the heat flux of 8 W/cm 2 . 

In the present pool boiling experiments with an imposed heat flux at the heater 
surface, the subcooling level, from near saturation to 11 °C, had little influence on 
AT* As is shown below, the present subcooling range also had little effect on the 

incipient boiling delay time, t*. 
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Figures 5.16 through 5.19 show the incipient boiling delay time, t*, as a 
function of the total heat flux. Figure 5.16 is for the damaged surface Q-l and Figure 
5.17 for heater surface Q-2. It may be noted that t* decreases with increasing q" T . It 

is also noted that the onset of a significant time interval between incipient boiling and 
boiling propagation is first observed at q x " of 4 W/cm 2 on Q-l and q T " of 6 W/cm 2 on 

Q-2. Further, the incipient boiling delay times with the surface orientation that is 
horizontal up are greater than those with a horizontal down orientation, for the same 
conditions. This is to be anticipated since the fluid motion induced by buoyancy 
transfers energy from the horizontal up heater surface at a rate greater than from the 
horizontal down heater surface. As a result, for a given total heat flux more time is 
required for the horizontal up heater surface to reach the heater surface temperature 
range in which incipient boiling is expected. On the other hand, no particular 
influence of bulk liquid subcooling on the incipient boiling delay time is to be 
observed. 

The incipient boiling delay times for three subcooling levels and with a/g = ± 1 
and microgravity are presented in Figures 5.18 and 5.19 for surfaces Q-5 and Q-6, 
respectively. Similar to the data with surfaces Q-l and Q-2, t* decreases with 
increasing q" T , and the experiments associated with a/g = + 1 have greater delay times 
for incipient boiling for the same heat flux than either a/g = -1 or microgravity for q" T 
< 7 W/cm 2 . It is also observed in Figures 5.18 and 5.19 that for q" T > 7 W/cm 2 , the 

values of t* for the heater in the horizontal up orientation approach the values of t* for 
the heater in the horizontal down position or in microgravity. Also, for q" T > 7 

W/cm 2 , boiling inception occurs before the onset of significant natural convection. As 
in the tests performed with surfaces Q-l and Q-2, the time interval between the onset 
of boiling and boiling propagation is greatest near 8 W/cm 2 . t* measured for surfaces 
Q- 1 and Q-2 in the horizontal down orientation for q" T > 4 W/cm 2 were similar to the 

t* determined from the drop tower experiments at the same heat flux. This 
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reproducibility of t* between the horizontal down tests and the microgravity tests 
indicates that the determination of the incipient boiling delay time of heater surfaces 
may be done with the heater surface in the inverted position rather than in the more 
expensive microgravity conditions. The slopes of the curves of the incipient boiling 
delay times for surfaces Q-5 and Q-6 were similar to each other, but different than that 
for surfaces Q-l and Q-2. Since surfaces Q-5 and Q-6 were sputtered at the same 
time, they should have physical characteristics similar to each other but different from 
those of surfaces Q-l and Q-2, which were vacuum deposited at the same time, as was 
pointed out earlier. 

The quartz surfaces are representative of "smooth" boiling surfaces, and the 
metal surfaces are representative of "rough" boiling surfaces. The smoothness of the 
quartz surfaces is characterized below using the measured surface superheats. The thin 
film surfaces were polished before sputtering with a 1.4 micron pitch polish by the 
manufacturer. Heater surface superheat levels ranged from 75 °C to 12 °C in the 
microgravity tests. The Laplace equation, equation 2.2, taken together with the 
Clausius-Clapeyron relationship for equilibrium conditions give an estimated cavity 
radii range of 0.02 microns to 0.8 microns. Scanning electron microscopy of the thin 
film on quartz surface Q-9 revealed minute scratches which had widths of 
approximately 0.2 microns. This width is in the above cavity radii range. 

5. 1.1. 6 Incipient Boiling and Boiling Propagation 

Of the earlier transient boiling studies described in Chapter 2, none had 
photographed simultaneous orthogonal view's of the heating surface, across and from 
beneath the heating surface, using high speed cameras. A single view of the heating 
surface is inadequate for distinguishing the characteristics of the vapor dynamics since 
much of the interfacial motion is obscured by other vapor motion. In the current 
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experiments, the boiling inception and the resulting vapor motion was observed from 
two orthogonal views: directly across the heater and through the bottom of the quartz 
substrate. Distinctly different categories of vapor dynamics were observed in both the 
laboratory and the drop vessel pool boiling experiments. This section will present 
types of behavior, some previously unobserved, which will be referred to here as 
boiling spread or boiling propagation. Boiling spread may be classified into two 
distinct categories: (a) The first is the formation or nucleation of vapor ahead of a 
moving liquid-vapor front, which is sometimes pushed ahead by this front and is at 
other times absorbed into the front, (b) The second is the growth of an interface 
similar to that of vapor bubble growth due to evaporation at the interface, regular in 
some cases and in other cases chaotic, but without the formation of vapor bubbles 
ahead of the interface. The parameters under which these differences in behavior 
occurred will be described, but their detailed origins are as yet unknown. 

The following discussion includes photographs which show the chronological 
sequences of the behavior observed. The measured initial conditions, including the 
applied power input, system pressure, bulk liquid temperature, a/g level, and 
orientation are given in the first figure of each chronological sequence. Six different 
categories of boiling spread are listed and then described below: 

A. Advancement of vapor by irregular protuberances 

B. Formation of mushroom-like vapor mass followed by spreading 

C. Single pancake-like bubble with a "smooth" interface 

D. Orderly growth followed by onset of interface disturbances 

E. Explosive growth of a vapor mass 

F. Marangoni motion of bubbles toward region of higher temperature 
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A. Advancement of vapor by irregular protuberances 

The first category of boiling propagation is the advance of the vapor mass by 
irregular protuberances from bubbles, shown in the sequence of Figures 5.20a through 
5.20e. The upper photo is through and from beneath the heating surface, and the lower 
portion presents the view across the heating surface. In the first figure of this series, a 
vapor bubble has already appeared. The onset of natural convection may be seen in 
Figure 5.20a, in both the view through the bottom of the heating surface and in the 
view across the heating surface. In the view across the heating surface, the onset of 
natural convection is characterized by the newly formed thermal plumes, parallel to 
and slightly above the heater surface. In the view from beneath the gold film heater, a 
disturbance around the periphery of the heater surface is visible. In the next figure, 
5.20b, three bubbles are now present, and the peripheral disturbance, a result of natural 
convection, may be seen in the view from beneath the gold film heater. In Figure 
5.20c, the boiling has spread over a large area of the heating surface, and bubbles 
depart from the heating surface as a result of buoyancy as the boiling propagates over 
the heating surface. The last figures of this boiling sequence are Figures 5.20d and 
5.20e. In Figure 5.20e, boiling covers most of the heater surface, and bubbles coalesce 
near the edge of the outwardly growing bubble front. This type of boiling spread was 
observed only with the heating surface in the horizontal up position at a/g = +1, and 
for the total power input to the gold film heating surface > 7 W/cm 2 . 

B. Formation of mushroom-like vapor mass followed by spreading 

The second category of boiling spread presented in Figures 5.21a through 
5.2 Id shows the formation of a mushroom-like vapor mass in the center of a thin 
vapor layer, which rises above the thin vapor layer due to buoyancy, while the thin 
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vapor layer moves outward beneath the large mushroom-like vapor mass. The first 
figure of this series, Figure 5.21a, shows the appearance of the first bubble in both 
views. The upper photo is through and from beneath the heating surface, and the 
lower photo of the figure looks across the heating surface. Convective stagnation lines 
are clearly observable in the lower portion of the photo. Figure 5.21b shows a thin 
vapor layer surrounding the large mushroom-like bubble. The bubble in the center of 
the thin outward growing vapor layer becomes larger in Figure 5.21c. The vapor layer 
has vapor bubbles on its surface. Figure 5.21d shows that boiling has covered most of 
the heater, and in the view across the heating surface, the large bubble is departing 
from the heater surface. The upward bubble motion entrains liquid and causes more 
disturbances over the outward growing vapor layer. This kind of spreading was 
observed only for a/g = +1 with 2 < q" x < 7 W/cm 2 . 

Okuyama et al (1988) observed this type of vapor motion using a large step in 
heat flux (q"| on the order of 20 W/cm 2 to 40 W/cm 2 ) applied to a 7 pm thick copper 

foil (0.4 cm wide by 4 cm long) that was laminated to an epoxy plate, in R 1 1 3 at P < 
0.5 MPa. A horizontal up heater orientation was used. In the present study, by 
contrast, the heat flux to the liquid is considerably smaller, with the largest q"| on the 

order of 1 .4 W/cm 2 . 

C. Single pancake-like bubble with a "smooth" interface 

A third category of boiling spread consists of a single pancake-like bubble with 
a "smooth" interface. In Figures 5.22a through 5.22d, using the laboratory pool 
boiling vessel, the heater surface orientation is horizontal down. The upper portion of 
each figure is a view through the heater surface from above, and the lower portion of 
each figure is a view across the heating surface. The initial bubble grew outward in an 
orderly manner from one location near the edge of the heating surface. This initial 
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bubble eventually moved away from the location at which it first originated, and a 
second bubble formed at this location, as shown in Figures 5.22a and 5.22b. In the 
underside view of Figure 5.22b, the smaller second bubble is in front of the larger 
bubble (in the darker region of the heater surface), but in the side view of Figure 
5.22b, the smaller bubble is behind the larger one. Figure 5.22c shows only one 
bubble in both views. The larger bubble of the previous two figures expended the 
available thermal energy for both bubbles while growing, and this large bubble 
eventually coalesced with the smaller bubble. The rate of growth of the larger bubble 
varied linearly with the square root of time, until the bubble reached the outer edge of 
the heater, an indication that the growth of the bubble was thermally controlled. After 
the bubble had grown past the edges of the film heater as shown in Figure 5.22d, it 
collapsed at some peripheral locations, and the liquid rewetted the gold film heater. 
This category of boiling spread was observed only in earth gravity with the heating 
surface oriented horizontal down, such that the upward buoyancy force flattened the 
bubble and held the vapor against the heating surface. 

D. Orderly growth followed by onset of interface disturbances 

The fourth category of boiling propagation consists of initially orderly growth 
of either a hemispherical bubble, in the case of microgravity, or a pancake-like one for 
a/g = -1. This orderly growth then gives way to violent growth, for some as yet 
unknown reason. The initially smooth bubble interface was observed to rupture either 
on the surface or at the base of the bubble at the beginning of the violent growth. 
Figures 5.23a through 5.23d show the history of violent boiling spread associated with 
the rupturing of the bubble occurring at the periphery of the bubble base. The bubble 
was initially smooth, but this is not shown. The heater surface shown is in the 
laboratory pool boiling vessel in the horizontal down position. The upper photo 
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presents the view across the heater surface, and the lower photo presents the view 
through the heater surface from above. Figure 5.23a shows a ruptured bubble one 
frame after boiling inception. In Figure 5.23b, the violently growing vapor mass 
develops large-scale satellite interface bubbles, which sometimes break away from the 
vapor mass. The view across the heater surface (oriented horizontal down) shows that 
the vapor is being pushed down into the subcooled liquid away from the heating 
surface. This motion was simultaneously accompanied by a spike increase in the 
pressure transducer reading. These two observations suggest that this later vapor 
motion is driven by pressure forces arising from the phase change. The heating 
surface is eventually covered by the vapor mass, as shown in the subsequent two 
photos. Figures 5.23c and 5.23d. 

The type of violent boiling spread by the rupturing of an initially smooth and 
symmetrical bubble surface is shown in the photo sequence of Figures 5.24a through 
5.24d. In these photos, the heater surface is again in the horizontal down position with 
the view through the heater surface from above in the upper portion of the photo and 
the view across the heater surface in the lower portion of the photo. Rupturing of the 
surface of the bubble is beginning in Figure 5.24a. The breaking of the liquid-vapor 
interface occurs on the bubble periphery near the warmest part of the heater, the heater 
center. The following two Figures, Figures 5.24b and 5.24c, show smaller bubbles on 
the outwardly moving interface of the larger bubble. In the last photo of this 
sequence. Figure 5.24d, the vapor mass with the bubbles at the liquid-vapor interface 
has nearly covered the heater surface. 

Category D phase change dynamics were observed to occur either in 
microgravity or in earth gravity with the heater surface in the horizontal down 
position. The heat flux producing category C boiling spread was greater than that 
associated with the category D boiling propagation, for the heater surface in earth 
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gravity in the horizontal down orientation, with the same subcooling and heater 
surface. 

E. Explosive growth of a vapor mass 

The fifth category of boiling spread is referred to as the explosive or energetic 
boiling propagation, characterized by initially uniform small scale protuberances over 
the entire liquid vapor interface. The sequence of photos in Figures 5.25a through 
5.25g are from the drop vessel experiment PBMT1 101.600 and illustrate the phase 
change dynamics associated with this kind of boiling propagation. In these figures, the 
upper photo of the figure is the view across the heating surface, and the lower photo of 
the figure is the view through and from beneath the heating surface. In the first of 
these figures, the upper part in Figure 5.25a, the conduction thermal boundary layer is 
visible as the line parallel to and slightly above the heater surface. This photo shows 
the onset of incipient boiling. In the view from beneath the heater surface, two slight 
shadows are seen on the heater surface. In the upper corresponding view at the same 
location, a shadow is present within or at the edge of the thermal boundary layer. The 
next photo. Figure 5.25b, shows boiling at five locations in the view from beneath the 
heating surface, and a large vapor mass above the thermal boundary layer in the view 
across the heater surface. Small scale protuberances are visible, but are blurred 
because of the high initial growth rate. The protuberances are more distinguishable in 
Figures 5.25c and 5.25d, and they are still present in Figure 5.25e. Figure 5.25f shows 
that the vapor has coalesced into a single large vapor mass, and some small bubbles 
have separated from this vapor mass. The protuberances are also larger. The vapor 
mass is observed to lift off the heater surface in the last photo. Figure 5.25g. The 
pressure disturbance associated with the initial rapid growth of the vapor mass imparts 
momentum to the liquid bulk which then lifts the vapor off the heating surface. 



114 


It is believed that the small scale protuberances, which may be the result of 
local pressure disturbances related to instabilities, greatly increase the liquid-vapor 
interface surface area, which then results in a large mass evaporation rate. The 
pressure disturbances resulting from this large mass flux create more protuberances 
which drive the unstable growth of the vapor mass across the heating surface. This 
type of boiling spread was observed only in microgravity, and only near saturation 
conditions with q j < 6.5 W/cm 2 . A lower limit for this boiling propagation category 

could not be established at this time because the time available in microgravity was too 
short for boiling incipience to take place at the low heat flux level. In the drop tower 
tests, a q T of 4 W/cm 2 did not result in incipient boiling. A space experiment is 

required to obtain information concerning the lower limit of total heat flux for this 
explosive type of boiling spread. This category of boiling spread will be discussed in 
more detail in Chapter 6. 

F . Marangoni motion of bubbles toward region of higher temperature 

The final type of boiling propagation observed in the pool boiling experiments 
consisted of small bubbles forming and then moving toward the center of the thin film 
heater, the warmest part, and coalescing there as shown sequentially in time in Figures 
5.26a through 5.26e. In these figures, the upper part of the photo is the view across 
the heater surface, and the lower part of the photo presents the view through the heater 
surface from beneath. The following photos resulted from the drop vessel pool boiling 
experiment PBMT013 1.820, with the test conditions as given in the first figure. In 
Figure 5.26a, the first bubble appears at the edge of the heater, and the thermal 
boundary layer is not yet visible as in the case of Figure 5.26b. In Figure 5.26b, the 
first bubble has moved away from the heater edge, and two more bubbles have formed 
at the original incipient boiling location. Another location on the left side of the heater 
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surface has become activated as a boiling site with the thermal boundary layer now 
visible. In the next photo. Figures 5.26c and 5.26d, a number of additional bubbles 
have emanated from the two boiling locations. Some of the smaller bubbles have 
coalesced with larger ones, and new bubbles form at the edges of the vapor mass. In 
the last figure. Figure 5.26e, one large vapor mass finally results, which covers the 
heater surface. Some vapor has separated from the vapor mass and formed individual 
bubbles floating in the R1 13. The growth of these bubbles was orderly and slow. This 
type of boiling spread was observed for microgravity and for a/g = -1. q" T for these 

instances was > 7.2 W/cm 2 . 

Table 5.1 summarizes the total heat flux ranges, subcooling levels, and body 
forces which characterize the conditions necessary for the aforementioned types of 
boiling propagation as known to date. 

5.1.2 Imposed Heater Surface Temperature 

5. 1.2.1 Comparison of Measurements and Conduction Prediction 

An imposed step in heater surface temperature, like an imposed step in surface 
heat flux, provides a simple and known temperature distribution in the test fluid, in the 
absence of buoyancy induced liquid motion. In order to observe the effects of 
different liquid temperature distributions on the resulting incipient boiling, these two 
different elementary liquid temperature distributions were used. The results with the 
imposed surface heat flux have been presented, and the results with a step increase in 
heater surface temperature will now be presented. The theory behind the design of the 
equipment for imposing a step in the heater surface temperature is presented in 
Appendix F, and a description of this equipment was given previously in Chapter 3. 
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Two different gold film surfaces, Q-6 and Q-14 were used in the tests 
conducted with an imposed increase in surface temperature. Both surfaces were 
heated in the horizontal down orientation, and the results from the two different 
surfaces were similar. Figure 5.27 presents the heater surface temperatures and 
uncertainties, as measured by the Codas system during the initial, rapid, portion of the 
heater power transient, and for the later part of the heater power transient as measured 
through the use of the CR7X. As discussed in Chapter 3, this later portion of the 
heater power transient was longer than the earlier portion. Because the initial series of 
steps in heat flux comprising an approximation to the heat flux necessary to impose a 
constant heater surface temperature were of short duration while the later steps were of 
much longer duration, two different measuring devices were used to verify that the 
system imposed a step in heater surface temperature. The less accurate Codas system 
measurements indicated that the surface resistance was constant, with an uncertainty in 
the resistance measurements of + 2%. This uncertainty in heater resistance translates 
as a large uncertainty in the measured heater surface temperature, on the order of + 7 
°C, resulting from the small slope of the resistance-temperature linear relationship. 
The CR7X provided more precision in the measured heater surface temperatures for 
the slower portion of the heater power transient, on the order of + 2.5 °C. The 
discussion concerning the estimation of the uncertainties in temperature measurement 
using the Codas and the CR7X is presented in Appendices C and D.l. The initial 
liquid temperature, T j( and the computed heater surface temperature are shown for 

reference. 

Figure 5.28 shows the comparison of the measured heater surface temperature 
with the computed one dimensional conduction heater surface temperature prediction, 
as presented in Appendix F. The measured heater surface temperatures for surface Q- 
6 shown in this figure were determined with the CR7X, and follow the one 
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dimensional prediction well until near the value of E~ npul /100 = 15. At this input 

voltage level, the larger disagreement between the measured and predicted heater 
surface temperature may be a result of multi-dimensional heating effects. The one 
dimensional predicted heater surface temperatures are greater than the measured heater 
surface temperatures, as would be expected, since the edge effects are neglected in the 
one-dimensional prediction. 

Figure 5.29 shows the incipient boiling delay times for the step increase in 
heater surface temperature experiments performed in the laboratory pool boiling vessel 
for surface Q-14. The initial system temperature was 49 °C, and the liquid was 
subcooled by increasing the system pressure above the saturation pressure for R 1 1 3 
corresponding to 49 °C. As the initial voltage E input , which corresponds to the 

imposed step increase in heater surface temperature, was increased, the delay time 
from the energization of the heater to the onset of boiling decreased somewhat. This 
decrease in the delay time to incipient boiling is expected intuitively since more 
energy goes into the liquid with a larger step in temperature. The delay times 
associated with incipient boiling are similar in magnitude for both surfaces tested here. 

Figure 5.30 presents a comparison between the incipient boiling delay times 
and heater surface temperatures at incipience for both imposed surface temperatures 
and imposed surface heat fluxes for surface Q-14. The initial system temperature was 
49 °C, and the liquid was subcooled by increasing the system pressure to a value 
greater than the saturation pressure corresponding to 49 °C. As demonstrated in 
Figure 5.30, the constant surface temperature experiments always resulted in boiling 
inception occurring at significantly lower surface temperatures and delay times when 
compared to the corresponding values of the constant heat flux tests. It thus may be 
concluded from these results that the manner of heating, and not merely the heater 
surface temperature, as viewed by some, affects boiling inception. Although Figure 
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5.30 only shows surface temperatures rather than surface superheat temperatures, the 
temperature differences are greater than those represented by different subcoolings. 

The type of boiling propagation observed with the imposed heater surface 
temperature experiments always consisted of a single pancake type bubble, category C, 
described in section 5. 1.1. 6 for boiling propagation observed in the constant heat flux 
experiments. This observation is related to the consistently small superheated thermal 
boundary layer formed in the imposed surface temperature tests and will be discussed 
in Chapter 6. 

5.2 Quasi-Steady Heating 

Quasi-steady heating was performed in the laboratory pool boiling vessel using 
both the gold film heater, Q-2 and the copper block heater, SH-2. The liquid 
temperature in the laboratory pool boiling vessel was 49 °C and the liquid was 
subcooled by increasing the liquid pressure above the saturation pressure 
corresponding to 49 °C. The intent of these tests was to find the minimum heater 
surface temperature and heat flux tor which the onset of boiling would just occur. The 
following figures show the heat flux to the liquid as a function of the difference 
between the heater surface temperature and the bulk liquid temperature. q" t was 

determined, as discussed in Chapter 4, through natural convection correlations. Figure 

5.31 is for surface Q-2 in the horizontal up orientation. Figures 5.32 and 5.33 are for 
the quasi-steady heating tests in the horizontal up and the horizontal down orientations 
involving the copper block heater, SH-2. Measurements at subcoolings of 1.23, 3.19, 
and 11.03 were obtained for the horizontal up orientation, and measurements at 
subcoolings of 4.08 and 1 1.87 were obtained for the horizontal down position. 

For surface SH-2 which had the gold-sputtered copper foil exposed to the test 
liquid, the heater surface superheat temperatures associated with the minimum power 
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input to the heater at which incipient boiling would just occur were in the narrow 
range of 39.0 to 44.8 °C for both the horizontal up and horizontal down orientations 
regardless of the nominal subcooling levels of 1, 3, and 11 °C. The heater surface 
superheat temperatures for surface Q-2 corresponding to the minimum power input to 
the heater at the point where incipient boiling would occur were in the range of 32.0 to 
39.6 °C irrespective of the nominal subcooling levels of 1, 3, 11, and 16 °C. The 
minimum heater surface superheat range required for incipient boiling was slightly 
greater for the metal foil surfaces rather than for the "smoother" quartz surfaces. This 
might at first appear contradictory to previous observations as described in Chapter 2. 
However, the measured heater surface temperatures of the present experiments are 
spatially averaged temperatures, and the local heater surface temperatures are 
estimated to vary by as much as 25 °C from each other, which can account for this 
discrepancy. (This variation in local heater surface temperature in earth gravity with 
the presence of buoyancy can be contrasted with the variation in local heater surface 
temperature in the microgravity tests with essentially no buoyancy. The local film 
temperature variation in microgravity was estimated to be less than 2 °C for locations 
1mm inward from the film heater edge.) The tests using the quartz surface and the 
metal foil surface were conducted with the same initial conditions. Prior to the 
experiments, both test surfaces were placed together in the laboratory pool boiling 
vessel as described in Chapter 3. 

For heater surface Q-2, the transient tests using an imposed step in heat flux 
resulted in larger mean heater surface temperatures at boiling inception for heat flux 
levels less than 7 W/cm 2 than in the case of quasi-steady heating for the same initial 
conditions. A step in heat flux of less than 7 W/cm 2 resulted in a mean gold film 
temperature of surface Q-2 at incipient boiling in the range of 42.0 to 62.9 °C. 

Figure 5.34 presents natural convection data obtained before boiling inception 
during quasi-steady heating. The Nu calculated from the heat flux using the thin film 
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surface had larger uncertainties than those computed for the metal heater. This larger 
uncertainty in the Nu of the quartz surface is a consequence of the relatively greater 
uncertainties in the corresponding heat loss estimates. The natural convection 
correlation of Merte and Clark (1961) describe the measured data better than equation 
4.7. Natural convection data is also shown for the metal heater surface in the 
horizontal down orientation. The Nu for the inverted metal surface are quite different 
from the Nu of the horizontal up orientation, as anticipated, because of the lateral fluid 
motion which results from buoyancy effects. 
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Figures 5.1 - 5.34 follow. 
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Figure 5.2 


Thin gold film temperature with onset of natural 
convection, a/g = +1 
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Figure 5.3 Predicted and measured time for the onset of natural 
convection, a/g = 10' 5 
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Figure 5.4 Thin gold film temperature wiih onset of natural 
convection, a/g = -1 
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Figure 5.5 Thin gold film temperature with heater surface 
dryout, a/g = 10* 5 
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Figure 5.6 


Departure of measured surface temperature from the ID 
semi-infinite media prediction 
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Figure 5.7 Effect of subcooling on the thin film surface 
temperature with a/g = +1 
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Figure 5.10 Representative pressure transducer measurements from the 
drop vessel at a/g = 10* 5 and a/g = +1 
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Figure 5.11 Representative measurements of the constant imposed total heat 
flux from the laboratory and drop vessels 
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Figure 5.12 Heater surface superheat at incipient boiling with constant 
imposed heat flux, surface Q-l 
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Figure 5.13 Heater surface superheat at incipient boiling with constant 
imposed heat flux, surface Q-2 
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Figure 5.14 Heater surface superheat at incipient boiling with constant 
imposed heat flux, surface Q-5 
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Figure 5.15 Heater surface superheat at incipient boiling with constant 
imposed heat flux, surface Q-6 
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Figure 5.17 Incipient boiling delay time with constant imposed heat 
flux, surface Q-2 
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Figure 5.19 Incipient boiling delay time with constant imposed heat 
flux, surface Q-6 
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Figure 5.20a Category A: appearance of the initial vapor bubble, time = 0.894 sec 




Figure 5.20b Category A: vapor bubbles immediately after the boiling 

propagation has begun, time = 0.934 sec 





Figure 5.20c Category A: convection cells with propagating vapor, 

lime = 1.126 sec 
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Figure 5.20d Category A: vapor propagating outward, time = 1.306 sec 
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Figure 5.20e Category A: vapor covering nearly the entire heating surface. 

Large bubbles visible from coalescing of smaller ones, time = 1.498 sec 



Thermocouple Above 
Heater Surface 


146 



Figure 5.21a Category B: appearance of the initial vapor bubble with 

convection cells, time = 0.237 sec 
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Figure 5.21c Category B: the thin vapor layer has advanced outward radially, 

time - 0.257 sec 
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Figure 5.2 Id Category B: heater surface is nearly covered by vapor mass and 

the large bubble departs upward from the surface, time = 0.289 sec 
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Figure 5.22b Category C: the two bubbles growing and moving from their 

site of origin, time = 0.419 sec 
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Figure 5.22c Category C: the two bubbles have coalesced into one, time = 0.505 sec 
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Figure 5.22d Category C: the single vapor bubble with a "wavy” interface 

remains at the heated surface, time = 1.151 sec 
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Figure 5.23a Category D: the single bubble has ruptured at the periphery of 

the bubble base and boiling propagation has started, time = 10.214 sec 
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Figure 5.23b Category D: large-scale satellite interface bubbles with the 

vapor pushed downward away from the heating surface, time = 10.226 sec 
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Figure 5.23c Category D: the heater surface is nearly covered by vapor, 

time = 10.238 sec 
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Figure 5.23d Category D: the single bubble remains at the heated surface 

with bubbles on parts of the larger bubble interface, time = 10.270 sec 
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Figure 5.24a Category D: rupturing of the bubble interface is beginning, 

time = 28.890 sec 
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Figure 5.24b Category D: vapor moves outward from bubble interface 

disturbance, time = 28.894 sec 






Figure 5.24c Category D: smaller bubbles are on the outwardly moving 

vapor mass, time = 28.910 sec 





Figure 5.24d Category D: the vapor mass with interface bubbles has nearly 

covered the heater surface, time = 28.946 sec 
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Figure 5.25a Category E: incipient boiling in the presence of the conduction 

thermal boundary layer, time = 1.130 sec 
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Figure 5.25b Category E: boiling at five different locations, time - 

4.813 sec 



Figure 5.25c Category E: the vapor mass covers nearly the entire heater 

surface, time = 4.8I8 sec 




Figure 5.25d Category E: the small scale protuberances are more 

distinguishable, time = 4.820 sec 
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Figure 5.25e Category E: the vapor mass with small scale interface 

protuberances engulfs vertical thermocouple wires, 
time = 4.823 sec 
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Figure 5.25f Category E: large scale protuberances are visible, 

time = 4.828 sec 
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Figure 5.25g Category E: liquid bulk momentum lifts the vapor mass off the 

heating surface, time = 4.973 sec 



169 




Side View 


Underside 

View 




Both Views 


Test Conditions 
Heater Surface Q-6 
a/g = IQ' 5 T s = 58.61 °C 


q" T = 7.8 W/cm 2 Ti = 47.50 °C 


AT sub =11.11 °C PBMTO 13 1.820 

P = 145.0 kPa Run # M-18-2 


Figure 5.26a Category F: the initial bubble at the edge of the heater surface, 

time = 0.513 sec 
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Figure 5.26b Category F: the initial bubble moves toward the heater surface 

center with two additional bubbles formed at the initial boiling 
site, time = 0.628 sec 
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Figure 5.26c Category F: more bubbles grow from both sites and move 

toward the heater surface center, time = 0.718 sec 



Figure 5.26d Category F: additional bubbles from the two boiling sites, 

time = 0.920 sec 
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Figure 5.26e Category F: one large vapor mass covers the heater surface and 

some vapor has separated from this vapor mass to form 
individual bubbles in the R113, time = 3.608 sec 
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Figure 5.27 Demonstration of transient thin gold film temperature measured 
with the Codas and the CR7X systems 
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Figure 5.28 


Comparison of predicted 
temperature for various 


and measured step rise in surface 
leveis of power input 



Delay Time, t*. (ms) 


177 



AT* SU p (Deg C) 

Figure 5.29 Measured incipient boiling delay times with step change in 
surface temperature 
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Figure 5.30 Measured incipient boiling delay times and the corresponoing 
surface temperature at boiling inception 
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Figure 5.31 Quasi-steady heating prior to incipient boiling for surface 
Q-2 with horizontal up orientation 
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Figure 5.32 Quasi-steady heating prior to incipient boiling for surface 
SH-2 with horizontal up orientation 
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Figure 5.33 Quasi-steady heating prior to incipient boiling for surfac 
SH-2 with horizontal down orientation 
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Figure 5.34 Natural convection heat transfer results 




CHAPTER 6 


DISCUSSION AND CONCLUSIONS 


In Chapter 5, measured transient surface temperatures were presented for two 
different basic thermal disturbances at the heat transfer surface: a step increase in 
imposed heat flux and a step increase in imposed temperature. In the absence of 
convection, these permit the computation of the temperature distribution within the 
liquid, using well-known mathematical techniques. The influence of the liquid 
temperature distribution at the onset of boiling on the resulting phase change dynamics 
will be demonstrated, using the calculated liquid temperature distribution, together 
with the measured heater surface temperature, both at the moment of incipient boiling, 
to explain the type of subsequent boiling propagation observed, especially the 
explosive type of boiling propagation occurring at the low heat flux levels. 
Additionally, an examination of previous concepts concerning an empirical 
heterogeneous nucleation factor will be presented. In conclusion, the underlying 
concepts of the present work will be summarized. 

6.1 Computed Superheat Temperature Distributions at Boiling 
Inception: Imposed q T " 

Figures 6.1 through 6.6 present the computed liquid superheat temperature 
distributions at the moment of incipient boiling together with the measured heater 
surface superheat temperatures in microgravity conditions obtained with the drop 
tower tests. The coincidence between the measured and computed heater surface 



184 


superheat temperatures confirms the accuracy of the calculated liquid superheat 
temperature distributions. The superheat temperature distributions, using the known 
total power flux input to the gold film, were calculated with the one dimensional semi- 
infinite media solution. For the case of a known imposed heat flux to the liquid, the 
expression for the transient one dimensional temperature rise is obtained by the 
Laplace transform solution of equation 2.19 as: 


2q", 


at 


AT =T(x) -T. = ~ \J— exp 


X2 ' 

q”. x 

' «v 

,4at. 

erfc 

i_ 

X 





(6.1) 


The initial test conditions, which were nearly identical for tests with two different 
surfaces, Q-5 and Q-6, are included in the figures in order to demonstrate how the 
results obtained may differ with different test surfaces. The tests were conducted at 
two or three levels of heat flux with three nominal subcooling levels, obtained by 
increasing the system pressure above the saturation pressure corresponding to the 
initial liquid temperature. 

For a given figure, the largest calculated superheated thermal boundary layer 
thickness at incipient boiling, 5*, is denoted for illustrative purposes in addition to the 
computed liquid temperature distributions at the time of boiling inception. 8 the total 
thermal boundary layer, or more simply the thermal boundary layer, is defined as a 
layer of liquid near the heating surface in which there is a significant increase in the 
liquid temperature above that far from the heating surface, with the liquid far from the 
heating surface remaining at the initial temperature. 5 is always > 6*. The 
superheated thermal boundary layer is the layer of liquid near the heater surface where 
the liquid temperature is greater than the saturation temperature corresponding to the 
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system pressure. Both the thermal boundary layer and the superheated thermal 
boundary layer are shown for illustrative purposes in Figure 6.7. 

Figures 6.1 and 6.2 are for the R113 initially near saturation conditions for 
surfaces Q-5 and Q-6, respectively. It may be noted that in both tests PBMT1012.800 
in Figure 6.1 with a q T " of 7.4 W/cm 2 , and PBMT0125.800 in Figure 6.2 with a q T " of 

7.6 W/cm 2 , which were conducted at the high heat flux levels, the lowest superheat 
distributions existed at incipient boiling. The superheated thermal boundary layer is 
the region within the liquid adjacent to the heater surface where the liquid temperature 
is above the saturation temperature corresponding to the initial, constant liquid 
pressure level, and it can be observed from these two Figures that the high heat flux 
levels have the smallest superheated thermal boundary layers. The boiling spread 
associated with these two tests was the "slow growth," category F type described 
earlier in Chapter 5. Tests PBMT1025.800, PBMT1026.600, PBMT1 102.800 with 
q x " of 7.5, 5.7, 7.7, W/cm 2 respectively, in Figure 6.1 and Test PBMT0207.600 with a 
q T " of 6.3 W/cm 2 in Figure 6.2 had the next larger liquid superheat distributions and 
the next larger superheated thermal boundary layers. These tests were characterized 
by the rapid unstable boiling spread pictured in Section 5. 1.1. 6, corresponding to 
category D. Test PBMT1 101.600 in Figure 6.1 with surface Q-5, was the most 
violent of all of the microgravity tests and was referred to as violent or explosive 
boiling propagation, category E, in Section 5. 1.1.6. This test resulted in a heater 
surface superheat temperature at incipient boiling nearly equal to that of tests 
PBMT1025.800 and PBMT1 102.800, but with a liquid superheat distribution 
considerably above these two tests and with the largest superheated thermal boundary 
layer of all of the microgravity tests. 

Figures 6.3 and 6.4 show the corresponding results using surfaces Q-5 and Q-6 
in R113, with a subcooling of approximately 2.8 °C. In Figure 6.3, tests 
PBMT09 14.805 and PBMT1013.805 had q x " of 7.8 and 7.2 W/cm 2 and in Figure 6.4, 
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both tests PBMT0201.805 and PBMT0223.805 had q T " of 7.7 W/cm 2 . These four 
experiments show the lowest superheat temperature distributions and the smallest 
superheated thermal boundary layers at the onset of boiling, and the subsequent 
boiling spread associated with these was of the "slow growth," category F type. On 
the other hand, the remainder of the tests in Figures 6.3 and 6.4 have large heater 
surface superheats and large superheated thermal boundary layers, and are 
representative of the rapid unstable type of boiling spread, category D. 

Figures 6.5 and 6.6 show the calculated liquid superheat temperature 
distributions and the measured heater surface temperatures at the onset of boiling for 
the microgravity pool boiling tests with the highest nominal subcooling of 11.1 °C, for 
surfaces Q-5 and Q-6, respectively. PBMT09 19.820 of Figure 6.5 with a q T " of 7.9 
W/cm 2 and of Figure 6.6 with a q T " of 7.8 W/cm 2 show the lowest superheat 
temperature distributions and the smallest superheated thermal boundary layers present 
in these two figures, and correspondingly had the "slow growth," category F boiling 
spread discussed earlier in Section 5. 1.1.6. The remainder of the tests in Figure 6.5 
demonstrate category D boiling spread. In Figures 6.1 through 6.6 for a given 
subcooling, 8* associated with boiling propagation categories D and E are always 
larger than the 8* corresponding to category F. 

6.2 Computed Superheat Temperature Distributions at Boiling 
Inception: Imposed T w 

Figures 6.8 through 6.10 present the computed liquid superheat temperature 
distributions and the measured heater surface superheat temperatures at the onset of 
boiling for tests in the laboratory pool boiling vessel in the inverted position (a/g = -1) 
only, m order to minimize convection effects and thus permit comparisons later with 
the constant heat flux microgravity tests. Taking the Laplace transform of equation 
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2.19 with a known imposed constant surface temperature up to the moment of boiling 
inception, the transient liquid temperature distribution can be computed according to: 


T(x,t)-T w 



T.-T 

1 W 



( 6 . 2 ) 


The tests conducted used three nominal subcooling levels, obtained by setting the 
system pressure to the appropriate level above the saturation pressure corresponding to 
the initial liquid temperature. The tests were all conducted with surface Q-14, and the 
specific conditions of the tests are given in the figures. The type of boiling spread 
observed in all of these experiments was always the single bubble pancake type, 
category C, previously described in Section 5. 1.1.6. q" T mMX of these figures is the 

initial step in heat flux which is imposed for 400 microseconds until the heater surface 
temperature is at the desired level, at which point q" T decreases. 

Figure 6.8 presents calculated superheat temperature distributions for imposed 
heater surface superheat temperatures of 4, 4, and 5 °C. The R113 was initially near 
saturation conditions with a T s of 48.9 °C, corresponding to the liquid pressure of 
these tests. A AT WS of 4 °C was the minimum level for which incipient boiling would 

occur. This minimum surface superheat for incipient boiling was repeatable within 1 
°C for the experiments shown here and in Figures 6.9 and 6.10. The 8* was 0.11 mm 
for experiment T107, and was the largest 5* of the tests in Figure 6.8. For all of the 
tests presented in Figure 6.8, the initial boiling always occurred at the same location, 
as observed in the high speed video pictures. 

Figure 6.9 presents calculated temperature distributions for imposed step 
increases in heater surface superheat temperatures of 3, 4, 6, and 9 °C. The test liquid 
was initially subcooled 2.3 °C at a nominal initial temperature of 49.1 °C. 3 °C was 


188 


the minimum heater surface superheat temperature at which incipient boiling just 
occurred. A 5* of 0.021 mm for test T106 was the largest superheated thermal 
boundary layer thickness for the experiments in Figure 6.9. The initial boiling always 
occurred at the same location, but was at a different location from that observed in the 
tests of Figure 6.8. 

Figure 6.10 presents the calculated temperature distributions for imposed step 
increases in heater surface temperatures of 1, 5, 9, and 13 °C. The test liquid was 
initially subcooled by 11.0 °C at a nominal initial temperature of 49.2 °C. 1 °C was 
the minimum measured heater surface superheat temperature at which incipient boiling 
would just occur. 5* was largest for test T1 14 with a value of 0.012 mm. The initial 
boiling was always from the same location, but was different from the initial boiling 
locations observed in the tests of Figures 6.8 and 6.9. 

The minimum heater surface superheat temperature required for incipient 
boiling near saturation conditions was 4 °C, while the minimum heater surface 
superheat temperatures required for incipient boiling at the nominal subcoolings of 2.3 
and 11.0 °C decreased to approximately 3 and 1 °C, respectively. The time delays 
from the onset of heating to boiling inception were much shorter than those observed 
for the imposed constant heat flux tests, being on the order of three to four 
milliseconds. It might be expected that with increasingly greater imposed steps in 
surface temperature that t* would be correspondingly smaller. This was not explicitly 
observed here because of the relatively large uncertainty in t* due to the camera 
framing rate. In the future, improved methods would be in order for more precise 
determination of t*. The superheated thermal penetration depth associated with the 
minimum surface superheat for inception, 5* min , decreased with increasing subcooling 
in the tests involving surface Q-14: near saturation conditions, 8* min was 0.027 mm; 
at a subcooling of 2.3 °C, 8* min was 0.016 mm; and at a subcooling of 1 1.0 °C, 5* min 
was 0.002 mm. Although the superheated thermal boundary layer of 0.002 mm is 
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quite thin, it still contains a large number of molecules, on the order of 7 x 10 18 for a 
heater surface area of 7.25 cm 2 . 

6.3 Comparisons of Computed Temperature Distributions: Imposed q T " and 


Figures 6.11 through 6.15 present the computed liquid temperature 
distributions at boiling inception for both the imposed constant step increase in surface 
temperature and the imposed step increase in heat flux, contrasting the relative 
magnitudes of the superheated thermal boundary layer thicknesses which persist 
following heating by these two methods. The delay times from the onset of heating to 
the appearance of the first bubble are also shown. The measurements were conducted 
in both the laboratory pool boiling vessel with the heater surface in the inverted 
position (a/g = -1) and in microgravity. 

Figure 6. 1 1 presents a comparison of the computed one dimensional superheat 
temperature distributions at boiling inception for both imposed step increase in surface 
temperature and imposed constant heat flux, for surface Q-14 with the bulk liquid 
initially saturated. The imposed surface superheat temperature levels shown in Figure 
6.1 1 are 4 and 1 1 °C. These two superheat levels were selected as representative of the 
other superheat levels conducted in the imposed surface temperature tests near 
saturation conditions. Since these tests were performed in earth gravity and the 
buoyancy force was present, the high level heat flux for the constant heat flux tests of 
7.8 W/cm 2 was chosen in order that boiling inception would take place before the 
onset of natural convection. In Figure 6.11, the 8* associated with incipient boiling 
for the imposed constant heat flux test, T207, was 0.82 mm while the largest of the 5* 
corresponding to the constant imposed surface temperature tests, T110, was only 0.03 


mm. 
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Figure 6*12 similarly shows the calculated one dimensional liquid superheat 
temperature distributions for imposed heater surface temperatures and imposed heat 
fluxes for surface Q-14 with an initial subcooling of 2.3 °C. Measured heater surface 
temperatures at the onset of boiling are also shown. Imposed heater surface superheat 
temperatures of 3 and 9 °C were considered representative of the other tests with the 
same initial conditions. The heat flux of 7.8 W/cm 2 with an incipient boiling delay 
time of 1.012 seconds, was used in order to minimize the effects of buoyancy on the 
computed temperature distribution. As can be noted in Figure 6.12, the 6* of 0.52 mm 
for T205, an imposed constant heat flux test, was considerably greater than the 5* of 
experiments T103 and T106, which were on the order of 0.02 mm. 

Figure 6.13 presents a similar comparison of the one dimensional superheat 
temperature distributions at the time of boiling inception with a nominal subcooling of 
11.0 °C, for surface Q-14 also. The representative imposed heater surface superheat 
temperatures of 1 and 13 °C with the corresponding liquid superheat distributions were 
plotted along with a test at the same initial conditions with an imposed heat flux of 7.8 
W/cm 2 . The boiling inception delay time for the constant heat flux test was 1.518 
seconds. Similar to Figures 6.11 and 6.12, the 6* for the imposed heat flux 
experiment again was greater than the 8* corresponding to the imposed heater surface 
temperature experiment. The 8* for T206 was 0.40 mm and the 8* for T1 14 was 0.01 
mm. 

Liquid superheat temperature distributions computed at boiling incipience but 
using a thin film heating surface different from surface Q-14 show a considerably 
larger difference between the superheated thermal boundary layer for the constant heat 
flux tests than for the imposed increase in surface temperature tests. Figures 6.14 and 
6.15 show comparisons of the liquid superheat distributions with both the imposed 
heater surface temperature tests and constant heat flux tests for another surface, 
designated as Q-6. Tests with the imposed heat fluxes were conducted under 
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microgravity conditions in the drop tower, and no question exists as to the absence of 
natural convection during the heating process, while tests with a step in surface 
temperature were conducted at earth gravity in the inverted position, as described 
previously. Figure 6.14 is for experiments near saturation conditions with two 
imposed heat flux levels, both conducted in microgravity conditions, of 6.3 and 7.6 
W/cm 2 . AT WS is given as zero for test T51 in Figure 6.14, but AT WJ cannot be zero as 

nucleation would not occur. This apparent contradiction of nucleation occurring at a 
superheat of 0 °C results from the heater surface temperature uncertainty, which is + 1 
°C. AT WJ must actually be greater than 0 °C, but the actual value is less than the 

uncertainty associated with the heater surface temperature measurement. The heat flux 
of 7.6 W/cm 2 has the smaller t* and resulting smaller superheated thermal boundary 
layer of the two imposed heat flux tests, but the superheated thermal boundary layer 
thickness, 0.15 mm, of the 7.6 W/cm 2 test is still greater than the calculated 
superheated thermal boundary layer thicknesses of 0.008 mm and 0.059 mm for tests 
T51 and T52, respectively. Figure 6.15 is for tests similar to those of Figure 6.14 
using Q-6, but with a nominal subcooling of 2.8 °C. The computed liquid superheat 
temperature distribution with two different heat flux levels, 4.9 and 7.7 W/cm 2 , are 
shown with the imposed increase in heater surface temperature tests. The constant 
heat flux tests were performed in microgravity, and the imposed surface temperature 
tests were performed with the heating surface inverted. Of the two imposed heat flux 
tests, the test with the higher heat flux of 7.7 W/cm 2 resulted in a lower heater surface 
superheat at the onset of boiling, and consequendy a smaller superheated thermal 
boundary layer. However, the superheated thermal boundary layer thickness of 0.2 
mm resulting from the 7.7 W/cm 2 imposed constant heat flux test is still considerably 
larger than the superheated thermal boundary layer thickness of 0.027 mm resulting 
from an imposed surface temperature. 
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In the tests using surface Q-14, 8* min decreased with increasing subcooling, but 
for the tests with an imposed surface temperature involving surface Q-6, 5* min 
increased with increasing subcooling. S* min for T50 with a subcooling of 2.61 °C was 
0.027 mm while 5* min for T5 1 near saturation conditions was 0.008 mm as presented 
in Figures 6.14 and 6.15. Table 6.1 lists computed 8* values in microns for imposed 

heat fluxes and imposed surface temperatures for surface Q-14 with similar initial 
conditions. 

6.4 Estimation of Change in Internal Energy 

It is conjectured that the type of boiling spread observed is related directly to 
the amount of thermal energy within the superheated thermal boundary layer. 
Superheated thermal boundary layers with more energy can be expected to exhibit 
more violent boiling propagation than those with less energy. With constant properties 
evaluated at an appropriate mean temperature, estimates of the increase in internal 
energy of the superheated thermal boundary layer up to the instant of incipient boiling 
can be made. 

The change in internal energy of the superheated thermal boundary layer per 
heater surface area up to the point of boiling inception, AU*, can be estimated by: 


8*(0 

AU* = p c J(T(x)-T)dx (6.3) 

0 

where S*(t) denotes the superheated thermal boundary layer thickness. The integration 
of equation 6.3 from x = 0 to x = 8*(t) for the change in internal energy of the 
superheated thermal boundary layer and from x = 0 to x = 5(t) for the change in 
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internal energy of the thermal boundary layer were performed using Simpson's Rule, 
using the appropriate analytical expressions for T(x) and values of 5*(t) and 5(t) 
corresponding to the measured onset of boiling. The values of the superheated 
penetration depth, internal energy change of the superheated boundary layer, total 
internal energy change, the ratio of the change in internal energy of the superheated 
thermal boundary layer to that of the total thermal boundary layer, and the type of 
subsequent boiling spread observed are presented in Table 6.2 for the microgravity 
tests where a constant heat flux was imposed. 

The experiments in which the increase in internal energy of the superheated 
thermal boundary layer was 19.150 kJ/m 2 or more at the moment of boiling inception 
are those in which the boiling propagation was either the rapid unstable type, category 
D, or the rapid violent type, category E. The rapid violent type of boiling propagation 
occurred only for tests in which the liquid was initially near saturation conditions. The 
tests with an increase in internal energy of the superheated thermal boundary layer of 
6.815 kJ/m 2 or less are characterized by the slow marangoni type of spreading, 
category F. Experiment PBMT1 101.600, which was the most violent of all tests, had 
the largest increase in internal energy of the superheated thermal boundary layer, 
46.68 kJ/m 2 , and also the largest ratio of increase in internal energy of the superheated 
thermal boundary layer to the increase in energy of the thermal boundary layer, 
0.9977. Two other microgravity experiments with category E type of boiling 
propagation, PBMT1025.800, and PBMT1 102.800, had a change in internal energy in 
8* of 33.22 kJ/m 2 and 31.28 kJ/m 2 , respectively, and a ratio of the increase in internal 
energy of 8 to 8* of 0.9974 and 0.9969, respectively. From the estimated change in 
internal energy, it is observed that the energy content for the most explosive type of 
boiling spread is the largest, 46.68 kJ/m 2 . The most rapid vapor bubble growth 
occurred in the microgravity experiment with a heater surface temperature of 120 °C. 
This experiment, PBMT1 101.600, had one of the lowest temperature gradients at the 
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heater surface, but had the largest superheated thermal boundary layer energy content 
of all the rmcrogravity tests. The boiling spread that occurred with relatively small 
increase in internal energy in the superheated thermal boundary layers were 
characterized by category F spreading, while the boiling spread that occurred with 
relatively large increases in internal energy in the superheated thermal boundary layer 
were of the rapid unstable type, or on occasion for very large increases in internal 
energy in the superheated thermal boundary layer, the explosive type. 

As can also be noted from a detailed examination of Table 6.2, the largest 
ratios of increase in internal energy of the superheated thermal boundary layer to the 
increase in energy of the thermal boundary occur with initially saturated liquid, as can 
be anticipated. The last two digits of the filenames in Table 6.2 provide the nominal 
liquid subcoolings in degrees Fahrenheit. 

In the case of the imposed surface temperature tests, incipient boiling began 
before the superheated thermal boundary layer thickness reached the values 
corresponding to the imposed heat flux tests. The more non-equilibrium conditions of 
a steeper temperature gradient and a smaller time scale associated with the imposed 
heater surface temperature tests may have made the liquid-vapor interface of the nuclei 
in the surface cavities more unstable than under the conditions of an imposed heater 
surface heat flux. In Table 6.1, the boiling delay time and the liquid temperature 
gradient at the surface both for 400 microseconds and the onset of boiling are 
presented. It is noted that at a time of 400 microseconds that the surface temperature 
gradient tends to be larger than that associated with the imposed heat flux tests. 
Particularly striking is the contrast in magnitude of the boiling delay times between the 
imposed heat flux tests and the imposed surface temperature tests. The surface 
temperature at incipient boiling was always greater for an imposed heat flux than for 
an imposed surface temperature for similar initial conditions and the same surface. 
Furthermore, from examination of the area under the temperature distribution curves. 
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Figures 6.11 through 6.15, the tests associated with an imposed heater surface 
temperature have a smaller increase in the internal energy of 5* as compared to the 
imposed heat flux tests, and thus category F boiling propagation is expected for the 
tests involving the imposed heater surface temperatures. 

The integration involved in the estimation of the internal energy with equation 
6.3 was found to be more convenient with an approximate polynomial representation 
of the liquid temperature distribution. If the temperature distribution is approximated 
with a quadratic polynomial, the thermal boundary layer thickness is found by the 
integral method as described in Arpaci (1966) to be: 

5(t)=*v/6at (6.4) 

and with the x coordinate axis shown in Figure 2.2a the liquid temperature distribution 
from the integral method is: 


(x - \f6at) 

_ if 

^ l 

AT. = T(x) - T. = (6.5) 

2k rz — 

i v6at 


The integration of the approximate solution, equation 6.5, is more convenient than that 
of the analytic solution, equation 6.1, with an error of approximately 2%. However, 
the integration of the exact expression, equation 6.1, was used in Table 6.2. 


196 


6.5 Bubble Shape at Onset of Rapid Unstable Spreading in Microgravity 


From Figures 5.23 through 5.25, the photos for the microgravity tests in which 
the imposed heat flux resulted in the rapid unstable and violent boiling spreading, 
categories D and E respectively, the bubble was initially smooth and symmetrical 
(This smooth and symmetrical bubble is not pictured here.), but later in time. Figures 
24a and 24b for example, the bubble ruptures at the liquid-vapor interface and then 
rapid spreading occurs. The correlation of the bubble shape at the moment of this 
rupture in dimensionless form with the dimensionless heat flux is examined here. 

For the physical situation in microgravity with no fluid motion, it is reasonably 
assumed that one dimensional conduction describes the heat transfer from the heated 
surface. For this situation the differential form of the energy equation is: 



2 

9x 


9T _1_ 
9t a 


( 6 . 6 ) 


with the boundary condition: 


q" = -k 


9T 


9x 


*=o 


(6.7) 


From equation 6.6, in terms of the physical dimensions of the situation: 

T T - 

or 8 „x 2 (6.8) 

where 5 is the diffusion length scale defined in terms of some characteristic 
temperature difference, x is the characteristic time scale. Figure 6.16 shows a small 
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bubble with D, the largest diameter parallel to the heater surface, and with H, the x 
coordinate dimension. The projected area of the bubble over the headng surface is 
much smaller than the total area of the heating surface, and the effect of the bubble 
presence, resulung in local cooling due to the phase change, on the remainder of the 
superheated thermal boundary layer is assumed to be negligible. 

From consideradon of the energy equation 6.6 and equations 6.7 and 6.8 it may 
be deduced that: 


D = fj(T, q", 8, k, H) (6.9) 

where T is a characteristic temperature, q" is the heat flux, and k is the liquid thermal 
conductivity. From the application of the Buckingham Pi theorem, three 
dimensionless groups are found: 


D_ f5 
H " 2 IH ’ kT . 


( 6 . 10 ) 


H appears on both sides of equation 6.10, and it is desirable to have the dependent 
variables on one side of the dimensionless relationship only. The two dimensionless 
groups on the RHS of equation 6.10 are multiplied resulting in two dimensionless 
groups: 


D 

H 



( 6 . 11 ) 
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T is a characteristic temperature which may be a temperature difference. With T w -T s 
as the characteristic temperature in this situation, 5 is taken as 5*, the superheated 
thermal boundary layer at the time just prior to rupturing, and: 




q"5* 1 


k(T -T ) 

W S 


( 6 . 12 ) 


Measurements of D, H, and (T w - T s ) were obtained at the onset of boiling for the 
experiments conducted in microgravity, and combined with the corresponding values 
of 5 in equation 6.12. The dimensionless shape of the bubble about to rupture was 
correlated with the heat flux imposed at the heater surface, which was made 
dimensionless in equation 6.12, and these quantities are plotted in Figure 6.17. The 
boiling spread became more violent as the dimensionless heat flux increased, which 
corresponds to the tests to the right of the dashed line. The lower limit of this 
relationship is characterized by the slow, category F boiling spread. 

6.6 Explosive Boiling Spread in Microgravity 

Following the appearance of a single bubble on the heating surface at a low 
heat flux level near saturation conditions, the explosive type of boiling spread, 
category E, was observed for the first time in the present microgravity experiments as 
described in Section 5. 1.1.6, Figures 5.25a through 5.25g and briefly referred to in 
section 6.5 in the discussion of the bubble shape and corresponding dimensionless heat 
flux. The rapidly growing vapor-liquid interface was covered with small scale 
disturbances which remained small at initial times, but these protuberances grew in 
size with time, as seen in Figures 5.25a through 5.25g. The rapidly growing vapor 
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mass with the interface disturbances pushed into the subcooled liquid located above 
the visible thermal boundary layer. Boiling inception was accompanied with a 
momentary spike in the measured pressure as indicated in Figure 5.10. 

In classical bubble growth theories, the bubble is assumed to grow from a 
smooth spherical vapor bubble in a uniformly superheated liquid, and heat transfer 
effects are considered in the liquid but thermodynamic equilibrium is assumed for the 
vapor. These theories describe three consecutive stages of bubble growth: very early 
surface tension controlled growth followed by inertia controlled growth which is, in 
turn, followed by a heat diffusion dominated stage. An example of such theories are 
presented in Prosperetti and Plesset (1978). Figure 6.18 shows two limiting conditions 
for vapor bubble growth, along with the measured growth rate of the vapor mass in 
microgravity for experiment PBMT1 101.600. The Rayleigh solution for bubble 
growth which is applicable to boiling, cavitation, and underwater explosions is shown 
as the upper limit on the growth rate. The Rayleigh solution applies to a spherical 
bubble growing in an infinite fluid, and assumes that the sphere expands radially 
outward producing a purely irrotational radial flow through the action of pressure 
forces alone. Panton (1984) gives the Rayleigh solution using the velocity potential 
and continuity equation in spherical coordinates with the Bernoulli equation. The 
boundary conditions specified that the fluid was stationary at infinity and that it 
followed the sphere motion at the sphere surface, with the resulting differential 
equation for the interface motion: 


P(R)-P- 

= R 




2 


P 


(6.13) 
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where R is the bubble surface radius and p is taken as a constant. If P(R), the 

saturation pressure corresponding to the heater surface temperature at the time of the 
appearance of the first bubble, and are assumed constant as a first approximation, 

then R(t) and R(t) are readily determined from equation 6.13. The step in pressure 
from P^ to the saturation pressure corresponding to the heater surface temperature is 

considered to provide a higher than expected growth rate. The solution for R(t) 
appears in Figure 6.18. The other analytical solution is for heat diffusion controlled 
growth of the vapor bubble, and provides the lower limit. During early stages of 
bubble growth the pressure forces are expected to control growth and at later times 
heat diffusion is expected to control growth. For the present case, neither of these 
solutions is closely followed. The effective radius used in the measurements of Figure 
6.18 is defined in Figure 6.19. Limitations in the measurements of the growth rates at 
the early time periods exist due to the relatively slow framing rate of the camera, 400 
fps, which are indicated by the data points in Figure 6.18. 

The type of boiling spread shown in Figures 5.25a through 5.25g is explosive 
and has not been previously observed for the conditions of no liquid motion during 
heating provided by microgravity. These photographs suggest that extremely large 
evaporation rates occur during early stages of the bubble growth which drive an 
instability at the bubble interface. A wrinkling of the liquid-vapor interface occurs 
which persists for some time after the vapor mass has covered the entire heating 
surface, Figures 5.25a through 5.25e. The pattern of protuberances covers the entire 
liquid-vapor interface, and the regularity of the pattern is striking. It is speculated that 
the presence of the small scale protuberances greatly increases the area of the liquid- 
vapor interface, which significantly increases the mass flux to the growing vapor 
bubble, which in turn drives the instability. It is also speculated that the wrinkles 
somehow arise from very rapid evaporation at the interface. 
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Similar small scale interface disturbances and large growth rates have been 
observed by Shepherd and Sturtevant (1981) and McCann et al (1988). In the 
experimental investigation of Shepherd and Sturtevant, a droplet of butane immersed 
in (immiscible) ethylene glycol in a bubble column apparatus at the homogeneous 
nucleation temperature was observed to vaporize explosively. At the periphery of the 
bubble, wrinkles were observed similar to those shown in Figures 5.25a through 5.25e 
and were believed to be driven by rapid instability. The wrinkles were speculated, but 
not proven, to be the result of the so-called Landau mechanism of instability, which is 
described in Landau (1944). Shepherd and Sturtevant used experimental pressure 
measurements to estimate the evaporation rate associated with the observed instability, 
and they found that the evaporation rate corresponding to this instability was two 
orders of magnitude greater than that predicted by conventional bubble-growth 
theories. Shepherd and Sturtevant believed that the large mass flux across the 
distorted liquid-vapor interface generated vorticity responsible for the Landau 
instability in the vapor by the flow transition at the interface. Vapor bubble growth in 
slightly superheated liquids has been investigated extensively, and this wrinkling of 
the interface was not seen earlier. The liquid superheat temperatures prior to the 
occurance of homogeneous nucleation, in the work of Shepherd and Sturtevant, were 
atypically large for butane, about 105 °C. 

The present experiments are different from those of Shepherd and Sturtevant in 
that the liquid in the present study was heated from a solid boundary in the presence of 
a comparatively large temperature gradient, and the heater surface temperature in the 
present study was below the homogeneous nucleation temperature for the given liquid 
pressure, in test PBMT1 101.600 (T w / T hom = 0.9). The heater surface temperature at 

incipient boiling for test PBMT1 101.600 was 120 °C, which corresponded to a heater 
surface superheat of nearly 75 °C. This heater surface superheat temperature at boiling 
inception for the present test conditions is large when compared to other observed 
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values of the heater surface superheat temperature at incipient boiling for R113 on 
glass in earth gravity, as described in You et al (1990) . 

6.7 Empirical Heterogeneous Nucleation Factor 


As discussed previously in Chapter 2, the form of the homogeneous nucleation 
rate used here was given by Fisher (1948): 

_ nkT -AG* 

J= ~ exp (“kf“) (6.14) 

In the presence of a solid heating surface, AG* can be expected to be considerably 
lower than that for homogeneous nucleation. This lowering of the thermal activation 
energy has been incorporated into the description of heterogeneous nucleation in terms 
of a dimensionless parameter defined as: 


AG 


het 


F= 


AG* 


(6.15) 


F is a positive number varying between zero and unity: zero implies no nucleation 
while unity corresponds to homogeneous nucleation. Kottowski (1973), Cole (1974), 
and Blander and Katz (1975) defined F in terms of a specific cavity geometry. In most 
situations, however, the actual cavity geometry and surface topography are seldom 
known, since typical boiling surface cavities are extremely small (on the order of 10 -7 
to 10 9 m). 

Since the cavity geometry is difficult to determine at present, Gerum et al 
(1979) and Ngheim et al (1981) attempted to define F empirically, as discussed 
previously in Chapter 2. Gerum calculated F values based on measured heater surface 
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temperatures during quasi-steady heating, and included boiling crisis data in his work. 
Ngheim et al also defined F empirically in terms of the measured heater surface 
temperature during transient heating at the instant of the appearance of the first bubble. 
Since the heterogeneous nucleation appeared to take place in the fluid immediately 
adjacent to the heater surface, the fluid temperature was taken to be the same as the 
heater surface temperature. The validity of using the measured heater surface 
temperature and measured system pressure to determine F is re-examined here. Since 
the present study is primarily concerned with transient incipient boiling, the work of 
Ngheim et al (1981) is discussed first, followed by a discussion applicable to both the 
work of Ngheim et al (1981) and Gerum et al (1979). 

In the work of Ngheim et al (1981), the dimensionless heterogeneous 
nucleation factor, F, was calculated based on the measured heating surface temperature 
at the instant of nucleation, and was plotted against the dimensionless temperature 
difference: 


T -T 

W S 

0= (6.16) 

T. -T 

hom s 


for various transient boiling nucleation data. 

F may be written from equations 2.4, 2.5, and 2.16 as: 
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nkT t* 
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t*, as defined earlier, is the time delay from the onset of heating of the liquid to the 
appearance of the first vapor bubble. In the development of equation 6.17, J was taken 
to be 1/t* as a reasonable first approximation. P v was taken approximately as the 

saturation pressure with a flat interface, corresponding to the heater surface 
temperature at nucleation. <j was evaluated at the heater surface temperature at 
nucleadon. 

t* may be neglected as a parameter in equation 6.17 since it appears in the 
logarithmic term. This can be demonstrated by expanding the logarithmic term as: 

nkT t* nkT 

W W 

In ( — h — ) = ln > + ln l * ( 6 . 1 8 ) 

The first term of the RHS of equation 6.18 is the natural log of a quantity with a 
typical magnitude of about 10 30 seconds, but t* of the second term of the RHS of 
equation 6.18 is much less than 10 30 seconds for most boiling situations. (As a matter 
of reference, 10 8 seconds is about 100 years.) 

To quantify this explicitly, consider R1 13 near saturation conditions initially at 
a temperature of 310 K, and heated to a temperature of 320 K at the surface. With the 
appropriate properties substituted into equation 6.17, F becomes: 

F = 7 x IQ 8 (93 + ln t*) (6.19) 


t* may be varied by as much as 10 5 seconds with a resulting change of only 12% in F. 
Similar behavior also occurs with liquid nitrogen and water. Hence, F, assumed to be 
a measure of the liquid-heater surface character, is a strong function of the heater 
surface temperature at nucleation, but a weak function of the delay time to nucleation. 
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and the nucleation delay time cannot necessarily be used as a parameter to characterize 
the heater surface. Thus t* may be neglected when: 


In t* 
nkT 


« 1 


( 6 . 20 ) 


The condition of equation 6.20 is expected to hold for typical engineering boiling 
applications for reasons given in connection with equation 6.18. 

Both the work of Ngheim et al (1981) and of Gerum et al (1979) rely on the 
concept that F, which results from the measured heater surface temperature and the 
liquid pressure, are representative of the heater surface-fluid character. Because of the 
exponential nature of J in equation 6.14, a narrow range of F may be calculated for a 
large range in J (over several orders of magnitude). To demonstrate the consequences 
of this behavior, J was set equal to 1 sec -1 cm' 3 (a reasonable value for the present 
experimental study) and curves for F were calculated for R113 at three different liquid 
pressure levels of 1, 2, and 3 atm. Figure 6.20 presents F, as defined in equation 6.17, 
as a function of 9, defined in equation 6.16, for R113 for these three pressures. T hom 

for each of the three pressure levels and the corresponding curves which show 
calculated data points are also given in the figure. A distinctly different curve results 
for each of the pressure levels, which implies that the relationship for F is property 
dependent, and not universal. The F curves, by definition, approach unity as 0 
approaches unity and approach zero as 9 approaches zero. F may indeed be a function 
of the solid surface and liquid-surface character and is not described by nucleation 
delay times, but perhaps from heater superheat temperature at nucleation, or some 
related function, together with equation 6.17. 
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In the work of Ngheim et al (1981) and Gerum et al (1979) and similar works, 
it was thought that there might be a general relationship between F and the 
dimensionless heater surface temperature, which was independent of the fluid type 
because F data calculated from measured heater surface temperatures for water and 
R113 heated with the same heater surface fell approximately on a single curve. In 
light of the present discussion, it is fortuitous that the calculated F values fell on a 
single curve. Considerable work remains for the more precise description of the 
nucleation characteristics of solid surfaces. 

6.8 Conclusions 

The main points of the present study are now summarized: 

1. Known temperature distributions in the heated liquid at the moment of boiling 
incipience, obtained in microgravity or for short times with the heater surface inverted 
in standard earth gravity, were shown to characterize the subsequent boiling 
propagation across the heated surface. Large superheated thermal boundary layers 
formed in the conditions of microgravity and low heat flux levels on a "smooth" 
surface are associated with a particularly explosive type of boiling spread previously 
unobserved, which is believed to be the result of an interfacial instability driven by a 
large mass flux across the wrinkled liquid-vapor interface. 

2. Buoyancy was shown to greatly affect the phase change dynamics associated with 
incipient boiling. Buoyancy, which produced fluid motion, affected the liquid 
temperature distribution during heating and the subsequent boiling propagation. 
Certain types of boiling spread, such as the explosive type of spreading, were observed 
only in microgravity, and not in earth gravity for either the horizontal up or the 
horizontal down heater surface orientations. Buoyancy induces convective motion 
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which reduces the energy of the superheated thermal boundary layer that would 
otherwise be available for boiling propagation. The buoyancy force carried vapor 
bubbles away from the heater surface when the heater surface was in the horizontal up 
orientation in earth gravity, but held vapor bubbles at the heater surface for the 
horizontal down orientation. 

3. The method of heating the liquid determines the temperature distribution in the 
liquid, which affects the resulting phase change dynamics significantly. The tests 
conducted with an imposed heat flux resulted in larger superheated thermal boundary 
layers than did the tests conducted with imposed step changes in surface temperature. 
As a result, the boiling spread after boiling incipience was markedly different. The 
temperature gradient in the fluid at the heater surface was greater with the imposed 
step change in surface temperature than in the fluid with the imposed heat flux, but the 
phase change dynamics with the imposed heat flux were more violent. In addition, the 
onset of boiling with the imposed step change in surface temperature occurred at 
significantly shorter delay times and lower heater surface temperatures than those for 
imposed heat fluxes. 

6.9 Recommendations for Future Study 

Because of the physical limitations of the drop tower, the microgravity tests in 
the present work were of short duration (about 5 seconds). For these tests, boiling 
propagation categories D and E at the low heat flux levels, on the order of 5 W/cm 2 , 
were particularly striking. It would be desirable to perform microgravity incipient 
boiling tests at lower heat flux levels, on the order of 0.5 to 2 W/cm 2 , for longer times 
than was possible with the drop tower experiments. As a result of the knowledge 
gained from the present experiments, low level heat flux experiments are being 
planned for a future space shuttle flight. 
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Figures 6.1 - 6.20 follow. 
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• PBMT 10 12.800, q" T = 7.4 w/an 2 , i* = 0.52 see 

+ PBMT1025.800, q” T = 7.5 W/cm 2 , t* = 2.467 sec 

a PBMT1 102.800, q" T = 7.7 W/an 2 , t* = 2.264 sec 

o PBMT1026.600, q" T = 5.7 w/an 2 , t* = 3.293 sec 

M PBMT1 101.600, q”T = 5.4 W/cm 2 , l* = 4.813 sec 
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Figure 6.1 Computed liquid superheat in microgravity from imposed heat 
flux at boiling inception, saturation conditions (surface Q-5) 
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Figure 6.2 Computed liquid superheat in microgravity from imposed heat 
flux at boiling inception* saturation conditions (surface Q-6) 


+++ 


♦ PBMT10 13.805, q" T = 7.2 W/cm 2 , t* = 0.045 Sec 

+ PBMT09 14.805, q"y = 7.8 W/cm 2 , t* = 0.300 Sec 

a PBMT1213.805, q" T = 7.6 W/cm 2 , t* = 1.728 Sec 

• PBMT1 128.805, q" T = 7.7 W/cm 2 . t* = 1.918 Sec 

a PBMT11 15.605, q" T = 5.8 W/cm 2 , t* = 3.083 Sec 

0 PBMT10 19.605, q" T = 5.3 W/cm 2 , t* = 4.083 Sec 

11 Measured Heater Surface Temperature 
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Figure 6.3 Computed liquid superheat in microgravity from imposed heat 
flux at boiling inception. 2.8 °C subcooling (surface Q-5) 


212 
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Figure 6.4 Computed liquid superheat in microgravity from imposed heat 
flux at boiling inception. 2.8 °C subcooling (surface Q-6) 
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Figure 6.5 Computed liquid superheat in microgravity from imposed heat 
flux at boiling inception, 11.1 °C subcooling (surface Q-5) 
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Figure 6.6 Computed liquid superheat in microgravity from imposed heat 
flux at boiling inception. 11.1 °C subcooling (surface Q-6) 
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heated surface 


Figure 6.7 The concepts of superheated thermal boundary layer and 
thermal boundary layer 
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Figure 6.8 
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Figure 6.9 Computed liquid superheat from imposed surface temperature at 
boiling inception, 2.3 °C subcooling (surface Q-14) 
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0.00 0.04 0.08 


X Distance From Heater Surface (mm) 


Figure 6.10 Computed liquid superheat from imposed surface temperature at 
boiling inception, 11.0 °C subcooling (surface Q-14) 



Liquid Superheat T - Ts , (Deg C) 


219 



X Distance From Heater Surface (mm) 


Figure 6.11 Comparison of computed imposed constant heat flux and 

constant surface temperature liquid temperature distributions, 
saturation conditions (surface Q-14) 
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Figure 6.12 


? 0 °nTanT 0 ™r°fL C °"! PUIed ' mp ° se<i con »ant beat flux and 
f.or'A ? temperature liquid temperature distributions. 
v_ subcooling (surface Q-14) 



Liquid Superheat T - Ts , (Deg C) 
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X Distance From Heater Surface (mm) 


Figure 6.13 Comparison of computed imposed constant heat flux and 

constant surface temperature liquid temperature distributions, 
11.0 °C subcooling (surface Q- 14) 
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* T51, q" T ^ = l.i W/cm 2 . AT wi = 0 Deg C 

+ T52 t q” T na, =14.9W/cm 2 , AT W * = 7 Deg C 

o PBMT0207.600, q" T = 6.3 W/cm 2 , t* = 1.685 Sec 

a PBMTO 125.800. q" T = 7.6 W/cm 2 . t* = 0.088 Sec 

□ Measured Surface Temperature 


Constant Heat Flux 
a/g = 10-5 

Boiling 
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Constant Surface Temperature 
a/g = -1 


Boiling 
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5* PBMT0207.600 


Nominal T s = 48.7 Deg C 


Surface Q-6 
Saturation Conditions 
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1.0 


2.0 


X Distance From Heater Surface (mm) 


Figure 6.14 Comparison of computed imposed constant heat flux and 

constant surface temperature liquid temperature distributions, 
saturation conditions (surface Q-6) 



Liquid Superheat T - Ts , (Deg C) 
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X Distance From Heater Surface (mm) 


Figure 6.15 Comparison of computed imposed constant heat flux and 

constant surface temperature liquid temperature distributions. 
2.8 °C subcooling (surface Q-6) 
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Figure 6.16 Small bubble neglecting the liquid temperature distribution 
immediately prior to boiling propagation 
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Figure 6.18 Comparison of measured bubble growth from test 

PBMT1 101.600 with the Rayleigh solution and thermal diffusion 
solution 
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Figure 6.19 Definition of effective bubble radius 



23i 



Figure 6.20 heterogeneous nucleation factor for three different 
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The uncertainty associated with the tabulated heat flux values are 3% of the given value. 



A. 2 Constant Heat Flux and Constant Surface Temperature: Pool Boiling 

Constant Heat Flux 
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B. CALIBRATION 

B.l Gold Film Heater Surface 

A massive copper fixture was used to calibrate the tempered gold film heater 
surfaces and is shown in Figure B.l. The thick copper insured uniform temperature 
conditions for calibration. The copper fixture, which was sealed with an o-ring, was 
immersed in a constant temperature water bath. The Colora constant temperature bath 
(model NB-35 640) provided the temperatures for calibrations over the temperature 
range 20 to 96 °C. Legs on the fixture positioned the surface about 5 cm from the 
bottom of the container. A calibrated chromel-constantan thermocouple positioned 
near the heater surface within the sealed fixture indicated the temperature of the gold 
film. The voltage of the thermocouple, which used a water-crushed ice reference 
junction, was measured with the CR7X data acquisition system. Once the surface 
reached a steady-state temperature as indicated by the thermocouple, the gold film 
resistance was ready to be determined. 

Since resistance is calculated from voltage measurements, rather than being 
measured directly, voltage measurements across the heater surface and across a 
calibrated shunt were used to calculate the heater surface resistance at an arbitrary 
constant temperature bath level. The voltage measurement circuit for calibration of 
the gold surfaces is presented in Figure B.2. A six volt battery provided the d.c. 
current source. The variable resistor was adjusted to a nominal value of 400 O so as to 
produce a small current which would, in turn, result in a negligible amount of ohmic 
heating in the gold film (0.01 watts), so that the film resistance would not be affected. 
A National Bureau of Standards calibrated standard shunt (Leeds and Northrup Model 
No. 4221-B) used for the determination of the current had a resistance of 0.1000008 
Q with an uncertainty of 0.001% at the current levels corresponding to this surface 
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calibration, as stated by the calibration certificate. The CR7X determined the voltages 
across the thin film and the shunt. A typical gold film calibration is shown in Figure 
B.3, and the coefficient of determination, described in Montgomery (1983), for this 
linear fit of the resistance and temperature is unity. 

B.2 Thermocouples 


Sections of a chromel-constantan thermocouple spool were calibrated 
individually. Copper extension wires were used between the thermocouple wires and 
the water-crushed ice reference junction, both during calibration and during use. Each 
set of thermocouple wire and extension wires were numbered for identification 
purposes. A Colora constant temperature bath provided the variable temperature for 
comparison with a platinum resistance thermometer, while a hypsometer at 
atmospheric pressure provided the steam point. Using a Mueller Temperature Bridge 
(Leeds and Northrup Model 8067) with a null detector, the resistance of the platinum 
element was measured to an accuracy of + 0.0001 12, which corresponds to an 
uncertainty in the temperature of the platinum resistance thermometer of + 0.001 °C. 
The output of the thermocouples were measured using a Leeds and Northrup K-5 

potentiometer with a built-in null detector. The sensitivity in measuring this voltage is 
within +0.1 (iV. 

B.3 Pressure Transducers 

A Ruska Pressure Dead Weight Calibration System (2400 Series dead weight 
gage) was utilized in the calibration of the Heise pressure transducers used in the 
laboratory boiling vessels and the Setra pressure transducer for the drop vessel. The 
uncertainty in the pressure level provided by the Ruska system was ± 0.0017 kPa (+ 
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0.0003 psi) for the pressure range of the present study: 96 to 172 kPa (14 to 25 psi). 
A schematic of the Ruska system and measurement circuitry for the pressure 
transducer calibration appears in Figure B.4. A representative calibration curve is 
shown in Figure B.5. The Ruska system was also used for calibrating a Heise 
Bourdon tube gage, which in turn was used for in-place periodic checks on the 
pressure transducer calibration. In Figure B.4, a precision pressure regulator was 
connected to the pressure transducer being calibrated, which was in turn connected to 
the dead weight system. Since the mass of the weights were known accurately along 
with an accurate determination of g, the pressure transducer signal measured by the 
multimeter or K-3 potentiometer could be correlated with the known pressure from the 
dead weight tester. 


CURRENT IN 
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AW 

VARIABLE RESISTOR 
SET TU 400 Ohms 


STANDARD RESISTOR* 
0.1000008 Ohms 



Figure B.2 The voitage measurement circuit for gold film calibration 
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Figure B.4 


Schematic for pressure transducer calibration 
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C. CR7X AND CODAS DATA ACQUISITION SYSTEMS 


Four tests were performed to determine the performance of the CR7X under 
various measuring conditions. The smallest voltages involved in this study were those 
of the chromel-constantan thermocouples. The accuracy desired for temperature 
measurement was 0.06 °C (0. 1 °F) which corresponds to about 3 JJ.V. It was found, in 
the first test, that an ice reference junction for thermocouples provided much more 
satisfactory temperature measurements than with the internal reference junction 
compensation of the CR7X. A set of negative and positive terminals were shorted 
with a copper wire to establish a condition of zero voltage in a second test. Voltage 
readings were taken employing both the fast and slow sampling rates of the instrument 
every two seconds at eight different ranges, from the 1.5 mV range to the 5 V range. 
It was found that two intermediate ranges, 15 mV and 50 mV, were the most desirable 
in terms of accuracy in temperature measurement in the 20 to 60 °C ambient 
temperature range. 

In the third test, the accuracy of the CR7X with various input voltages for 
different ranges and sampling rates was verified. The accuracy at each range and 
speed was measured using a K-3 potentiometer (Leeds and Northrup Model 7553-5), 
which had a maximum sensitivity of 0.5 |iV, and a null detector (Leeds and Northrup 
Model 9834-1), which had a maximum sensitivity of 0.1 microvolts per division, as 
the voltage standard. The desired maximum deviation of 3 jiV was not met for all of 
the ranges. The ranges of 500 mV to 5000 mV in the slow sampling rate mode did not 
meet the 3 |iV criterion for any voltage level. The 1.5 mV and the 5 mV range in the 
slow sampling rate mode were unable to measure voltages exceeding 5 mV. Hence, 
their effectiveness was diminished. For thermocouple measurements, the 15 mV and 
the 50 mV ranges were used at the slow sampling rate which had an uncertainty of ± 
0.3 (iV. The maximum uncertainty on the fast sampling rate mode was ± 0.6 (J.V for 
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the 15, 50, and 150 mV ranges at the 49 °C level, which corresponds to a maximum 
uncertainty in temperature of about ± 0.01 °C. 

The dynamics of the CR7X were evaluated in a fourth test against a 7700 series 
eight channel Sanborn recorder. The Sanborn (with a Model 8803 high gain d.c. 
amplifier) and the CR7X were used individually and sequentially to measure the 
transient shunt voltage near the voltage level expected in the transient pool boiling 
measurements. For this test, the high gain amplifier of the Sanborn recorder had a 
sensitivity range of 1000 jiV per division to 1 ^iV per division. The Sanborn pre-amp 
followed from 10 to 90% of a step change in 5 ms. The CR7X and Sanborn voltage 
measurements were in excellent agreement, as shown in Figure C.l. 

Tests were performed with both the CR7X and the Codas data acquisition 
system in the 50 to 150 mV range, the range of the divided voltage difference across 
the heater surface, in order to estimate the uncertainty in the Codas system voltage 
measurement. In this voltage range, there were differences of about one millivolt in 
the voltage readings between the Codas and the CR7X with the same d.c. source. As 
a result, the uncertainty in the Codas voltage measurements in this voltage range were 
estimated to be + 1 mV. 
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Figure C. 1 Comparison of (he CR7X and the Sanborn dynamic voltage 
measurement 
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D. UNCERTAINTY ANALYSIS 


The uncertainties associated with the calculation of the gold film temperature, 
thermocouple temperature, system pressure, and total heat flux are estimated following 
the procedure of Kline et al (1953). According to this procedure, the uncertainty u r in 
computing r is given by: 


u 


2 ar 2 ar 2,11/2 

sy +( ^ + - + £v 

i 2 n 


(D.l) 


where r is a function of z lt .... 2 ^ and u Zi , u^, u^ are the uncertainties in each of 
these values, respectively. 


D. 1 Gold Film Temperature 


The CR7X was used for calibration of the gold film heater surfaces and for the 
imposed constant heat flux tests. The uncertainty in the heater surface resistance, 
which ultimately results in the heater surface temperature, was less for calibration of 
the heater surfaces than for actual power tests. From equation D.l the uncertainty in 
the gold Film resistance is expressed as: 
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Values of the quantities in equation from a representative calibration are: 
uv = 0-6 p.V uv = 0.6 |iV 

T w T sh 

V w = 120 mV V sh = 50 mV 

UR jh = 0.00 IQ 

R sh = 3.760 a 

These values substituted into equation D.2 result in an uncertainty of 0.03%. For 
power runs, the above values are different and result in a greater uncertainty. 
Representative values for a power run are: 

uv = 0.05 mV uv = 0.6 jiV 

T w T sh 

V w = 120 mV V sh = 50 mV 

u R = 0.000002 n 

AV sh 

R sh = 0.016314 Q 


with a resulting uncertainty in the heater surface resistance of 0.1 1%. 

Calibration data with the largest scatter would be expected to provide a 
conservative estimate of the uncertainty in the slope. From the calibration data of 
surface Q-2, the corresponding temperature and resistance data with the largest 
deviation from the Fit line are 54.85 °C at 2.8801 Q and 83.59 °C at 2.0172 £1 An 
estimate for the uncertainty in the slope of the linear relationship between the surface 
temperature and resistance can be found by using equation D.l together with the 
definition of this slope and the above values for temperatures and resistances, 
respectively. The uncertainty in the slope for these values is ± 0.23 °C/Q . 

The uncertainty in the surface temperature measurement can now be estimated 
with a known estimate of the uncertainty in the heater surface resistance and of the 
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slope in equation 4.1. Carrying out the indicated operations of equation D.l on the 
equation for the surface heater temperature, equation 4.1, the following expression for 
uncertainty in the temperature measurement results: 
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(D.3) 


Representative values for the quantities in equation are: 

u Tc = 0.05 °C URw = 0.003 Q 

u m = 0.23 °C/Q URc = 0.003 Q 

These values result in an estimated uncertainty in the surface temperature 
measurement of about ± 1 °C using the CR7X data acquisition system. The 
calculations were repeated in the same manner for the Codas system, and the 
uncertainty for representative values of the above quantities was found to be + 6 °C. 

D.2 Thermocouple Temperature 

As discussed in section B.2, the thermocouple spool was calibrated against the 
platinum resistance thermometer which had an uncertainty in temperature of + 0.001 
°C, and a calibration equation for the spool was determined from a fourth order 
polynomial least square fit of the measured thermocouple voltage and the determined 
platinum resistance thermometer temperature. In terms of the uncertainty in the CR7X 
voltage measurement, which was ± 0.6 JJ.V, the uncertainty in the temperature 
determination is estimated to be + 0.01 °C. 
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D.3 System Pressure 

There was little discrepancy, on the order of + 0.005 psi between the Ruska 
pressure determination, with an uncertainty in pressure of ± 0.0003 psi and the 
calculated pressure resulting from the calibration equations. The uncertainty in 
pressure resulting from the use of these equations was then taken as ± 0.005 psi. From 
equation D.l: 


Up 

~F = 

where P is the calculated pressure from the calibration equation and V is the pressure 
transducer voltage signal. For a typical run the following values were obtained: 

uy = 1 mV P = 16 psi 

dP -2 

— = 9.8608 - 2.7352 x 10 V V = 1 .6 V 

oV 

The uncertainty in pressure at this level is estimated from equation to be 0.06%. 

D.4 Total Heat Flux 

The expression for q" T was given earlier in equation 4.2. Using equation D.l 


'8P 


dV 


V 

~P~ 


(D.4) 


with equation 4.2, an expression for the uncertainty in q" T is: 
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In the experiments using 

uy T = 0.3 

V T = 5 V 
uy = 0.3 jlV 
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V sh = 0.08 V 
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the following were typical values: 

u Rjh = 0.000002 Q 

R sh = 0.016314 Q 
u^ = 0.08 cm 2 

A = 7.25 cm 2 


(D.5) 


With equation D.5 and these values, the estimated uncertainty in the calculated q" T is 

1 %. 
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E. ELECTRONIC CIRCUITRY 

E. 1 Power Supply for Step in Heat Flux and Quasi-Steady Heating 

As shown in Figure E.l, the power controller receives the set point voltage 
from either the CR7X or a manual input command, and amplifies the signal by two. 
The feedback signal (FBK1) from heater card number 1, presented in Figure E.2, is 
then subtracted from the amplified set point voltage and extraneous noise is filtered. 
For a positive difference, a proportional amount of current is sent into the base of 
transistor number 1. For a negative or zero difference, a diode prevents current from 
entering the base of transistor number 1. The base current then draws an amplified 
amount of current from the collector which is connected to bases (Bl, B2, B3, and B4) 
of the four transistors on heater card number 1. Thus, as the base currents to the heater 
card transistors arc increased, more power is sent to the heater. 

E.2 Power Supply for Step in Heater Surface Temperature 

This power supply produces a rapidly varying, programmed schedule of power 
to the test heater to result in a constant temperature of the test heater. The schedule 
begins with the maximum power to the heater, and the power is progressively 
decreased. 

Discrete changes in the heating power permit the use of switches. N-channel 
MOSFET switches are used. The heater is on a common bus which connects the 
sources of all of the MOSFETS. The power supply voltage is connected to a separate 
loading resistor for each switch, the other side of the load resistor is connected to the 
drain of the MOSFET. The series combination of the heater and load resistances 
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determines the power to the heater when the appropriate switch is on. Each gate is 
optically coupled to its particular switching signal. 

The voltage-time schedule is divided into 24 different time intervals; as short 
as 20 micro-seconds for early intervals and as long as 1.3 seconds for the final interval. 
A micro-processor programmed in assembly times the signals to the switches. 

Since the power supply starting transient requires about 50 milliseconds, an 
output resistor with a MOSFET permits the power supply to operate near maximum 
power before the actual power switching transient is imposed. Details of the electronic 
circuitry are presented in Figure E.3. 

E.3 Pressure Controller 

For large deviations from the set point pressure level, the pressure controller 
operates as an on/off controller by opening and closing the solenoid actuated fill/vent 
valves in response to large pressure deviations. More precise control is achieved 
through regulation of the proportional valve while a needle valve admits a constant 
flow rate of air into the bellows and the fill/vent valves remain closed. The 
proportional valve is operated within its linear range of operation so that as voltage to 
the proportional valve is increased or decreased, the flow rate of air exiting the bellows 
in direct proportion to the voltage. The electronic circuits are presented in Figures E.4 
through E.6. 

E.4 Drop Vessel Temperature Controller 

Before a drop test, the drop vessel was brought to a uniform temperature by 
means of heating films attached to the external surfaces of the drop vessel. The 
heating films were powered by means of the drop vessel temperature controller. The 
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electronic circuits for the drop vessel temperature controller are presented in Figures 
E.7 through E.10. 



Figure E.l Main healer controller 
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Figure E.2 Heater controller 
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Figure E.3 Heater surface power controller for imposed step in 
temperature 












Figure E.5 Pressure controller: 


part two 
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+15V 10 K 



Figure E8 Drop vessel external heater power controller: part two 
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F. ANALYSIS OF IMPOSED STEP IN TEMPERATURE 


One dimensional conduction to two semi-infinite media, in this case R113 and 
quartz, from the common plane heat source, which produces a step in the interface 
temperature at the boundary, are described analytically. Die problem is considered in 
two domains as shown in Figure F.l. The following governing equation, initial 
condition, and boundary conditions are for the R1 13: 

90 920 

Tt =a ^ (F-i> 

where 9 = T - T w 

0(x,O) =T m -T w (F.2) 

6(0, t) = 0 (F.3) 

0(°°,t) = T„-T w (F.4) 

The solution to F.l and the heat flux necessary at the boundary to maintain the 
constant temperature are by the Laplace transform method as described in Arpaci 
(1966): 
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The temperature distribution in the quartz is also given by equation F.5, but with the 
properties for quartz substituted for the R 1 1 3 properties and the origin for the quartz 
domain as indicated in Figure F.l. The heat flux to the quartz is given by equation F.6 
with the appropriate quartz properties. The sum of the heat flux to the quartz and the 
heat flux to the R 1 1 3 results in equation 3.1, the total power input to the heater surface 
divided by the surface area to impose a step increase in surface temperature: 


T -T 
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G. ANALYSIS OF IMPOSED STEP IN HEAT FLUX 


The transient temperature distribution resulting from conduction in two 
connecting semi-infinite media with a uniform plane heat source q T " at the interface 

may be found analytically for each of the two media. The solutions for the R113 and 
quartz substrate are, respectively: 

x . 

2 


(G.l) 
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H. THREE DIMENSIONAL CONDUCTION NUMERICAL ANALYSIS 

A three dimensional numerical solution for the transient temperature 
distribution in the quartz and the R113 was obtained by solving the conduction 
equation with temperature dependent properties and power generation in the common 
plane. The Crank-Nicholson finite difference method, which is unconditionally stable, 
was used to solve the equations. The numerical grid had equivalent node spacing in all 
three coordinate directions, and the spacing was adjustable. Since the problem was 
symmetric, the temperature in one fourth of the quartz and R113 soluuon domain 
provided the three dimensional temperature distribution. Temperature solutions in 
several vertical planes of the quartz and R113, including the heater surface, were 
printed at the desired time. Boundary conditions on the bottom of the quartz surface 
included those of convection and an insulated boundary. The edges of the heater 
surface and the R1 13 above were assumed to be insulated. 
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